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1.  OVERVIEW 
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1 . 1  INTRODUCTION 


Earth’s  neutral  atmosphere  affects  many  aspects  of  Ail  Force 
operations.  These  include  satellite  tracking,  lifetime  and  reentry 
prediction,  design  and  operation  of  aerospace  systems  and  the 
specification  of  absolute  density  for  point  analysis.  Requirements 
for  knowledge  of  absolute  density  and  aerodynamic  drag  variations 
within  5-10%  are  not  met  by  current  models.  Evaluations  of  models 
above  150  km  have  shown  that  despite  improvements  in  our 
morphological  description  of  the  atmosphere,  model  accuracies  for 
drag  variability  have  remained  at  about  +  15%  (one  sigma)  during 
the  past  two  decades.  This  plateau  has  been  dubbed  "the  15% 
barrier" . 

The  objective  of  this  workshop  was  to  assess  our  capability 
to  significantly  improve  models  of  density,  winds,  temperature  and 
composition.  Emphasis  was  focussed  on  the  altitude  region  80-250 
km  where  there  continues  to  be  a  lack  of  coordinated  neutral 
density  programs.  The  meeting  format  was  similar  to  that  initiated 
at  NOAA's  Satellite  Drag  Workshop  in  1982  and  followed  at  the  NASA 
MSFC  Drag  Workshop  in  1985.  There  was  an  extremely  strong  response 
to  this  meeting  from  the  scientific  and  user  community.  The  103 
registered  participants,  listed  at  the  end  of  this  report,  included 
representatives  from  22  Government  agencies,  7  Universities,  6 
research  centers  and  7  private  organizations. 

Section  1.2  gives  the  meeting  agenda.  Papers  submitted  by 
workshop  speakers  are  in  Sections  2-8.  Some  of  the  titles  differ 
from  those  in  the  agenda  and  a  few  papers  were  not  received. 
Through  presentation  and  discussion,  the  three-day  meeting  provided 
extensive  coverage  of  the  following  areas:  Air  Force  relevance, 
atmospheric  physics,  empirical  and  theoretical  models,  measurement 
techniques,  drag  coefficients  and  solar  influences.  A  meeting 
summary  is  given  in  Section  9.  It  was  concluded  that  significant 
improvements  in  our  80-250  km  modeling  capability  are  achievable. 
A  team  approach  to  this  multi-faceted  problem  was  recommended.  The 
major  requirement  is  to  obtain  a  comprehensive  data  base  of 
accurate  neutral  atmosphere  measurements  and  associated  forcings 
to  extensively  test  and  improve  formulations  in  thermospheric 
general  circulation  models. 


(Received  for  publication  23  January  1990) 
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AIR  FORCE  GEOPHYSICS  LABORATORY 
IONOSPHERIC  PHYSICS  DIVISION 


WORKSHOP 


ATMOSPHERIC  DENSITY  AND  AERODYNAMIC  DRAG 
MODELS  FOR  AIR  FORCE  OPERATIONS 


20-22  OCTOBER  198  7 

HOLZMAN  SEMINAR  ROOM 
BUILDING  1100A 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731-5000 


MAJOR  OBJECTIVE  IS  TO  ASSESS  CAPABILITY  FOR 
SIGNIFICANT  IMPROVEMENT  IN  MODELS  OF  DENSITY. 
WINDS,  TEMPERATURE  AND  COMPOSITION 
-  FOCUS  ON  ALTITUDE  REGION  80-250  km 


FINAL  AGENDA,  SCHEDULE  AND  AREA  INFORMATION 
ARE  ATTACHED 


ORGANIZER:  FRANK  A.  MARCOS 
IONOSPHERIC  PHYSICS  DIVISION,  AFGL 
TELEPHONE.  (617)  377-3037 
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FINAL  AGENDA 


WORKSHOP  ON  ATMOSPHERIC  DENSITY  AND  AERODYNAMIC  DRAG 
MODELS  FOR  AIR  FORCc  OPERATIONS 

AIR  FORCE  GEOPHYSICS  LABORATORY 
20-22  OCT  87 

**  TUESDAY,  20  OCTOBER  1987  ** 

I.  INTRODUCTION 

TIME:  0900-1000 

-  CHAIR:  F.  A.  Marcos  (AFGL ) 

1 .  Welcoming  Address 

2.  Meeting  Overview  and  Objectives 

II.  AIR  FORCE  REQUIREMENTS 
TIME:  1000-1145 

-  CHAIR:  (Capt.  T.  Schott,  HQAWS) 

1.  Air  Weather  Service's  Neutral 
Density  Model 

2.  Point  Density  Specification 

3.  National  Aerospace  Plane 

4.  Issues  on  Atmospheric  Drag 
Models  in  Operational  Systems 

5.  I onosphere- Thermosphere  System 

LUNCH:  1145-1315 
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Col .  J.  R.  K idd 
Commander,  AFGL 
H.  Carlson  (AFGL/LIS) 


Capt.  T.  Schott  (HQAWS) 

C.  Schaaf  (AFGL/LYA ) 

Lt.  S.  Johnson  (ASD/WE) 
J.  Liu  ( SPACECMD ) 

H.  Carlson  (AFGL/LIS) 


III.  PHYSICS  OF  DENSITY  AND  DRAG  VARIABILITY 
TIME:  1315  1450 
-  CHAIR:  R.  G.  Roble  ( NCAR ) 

1.  Mesosphere  -  Lower  Therm.  Coupling 

2.  Gravity  Waves  in  Mesosphere-Therm. 

3.  Tidal  Effects  on  Therniospher ic 
Circulation,  Temperature  and 
Density  Structure 

4.  Thermosphere  Dynamics  and  the 
Influence  on  Atmospheric  Density: 

The  !.E  Experience 

5  .  Session  Summary 

BREAK:  1450-1510 


J.  Forbes  (Boston  Univ.) 

R.  Walterscheid  (Aerospace  Corp) 
C.  Fesen  (Univ.  Colol 


T  Killeen  (Univ.  Mich) 


R.  Roble  (NCAR) 


IV.  NEW  DENSITY  AND  DRAG  MEASUREMENTS/REMOTE  SENSING 

A.  GROUND-BASED 
TIME:  1510-1655 


-  CHAiR:  P.  B.  Hays  (Univ.  of  Michigan) 


1 .  AFGL  L idar  Faci 1 ity 

2.  Lidar  Density  Measts 

3.  Radar  Winds/Gravity  Waves 

4.  Incoherent  Scatter  Density 

5.  High  Lat.  Lower  Therm.  Dynamics 

EVENING:  L  idar  Facil ity  Tour 
**  WEDNESDAY,  21  OCTOBER  1987  ** 


D.  Sipler  (AFGL/L 1 0 ) 

C.  R.  Philbrick  ( AFGL/L ID) 

R.  Vincent  (Univ.  Adelaide) 

W.  Oliver  (MIT  Haystack  Obs.) 
R.  Johnson  (SRI  Int‘  1 ) 

(D.  Sipler,  AFGL/L I D ) 


IV.  NEW  DENSITY  AND  DRAG  MEASTS/REMOTE  SENSING  (Cont'd) 
TIME:  0830-1000 


-  CHAIR:  P.  B.  Hays  (Univ.  of  Mich) 

B.  SATELLITE  REMOTE  SENSING 

1 .  Density  80-120  Km 

2.  Determination  of  Atomic  Oxygen 
Density  and  Temperature  of  the 
Thermosphere  by  Remote  Sensing 

3.  RAIDS 

C.  COUPLED  ENERGETICS  AND  DYNAMICS 
ATMOSPHERIC  RESEARCH 

-  Session  Summary 
BREAK:  1000-1020 


P.  Hays  (Univ.  Mich) 

R .  Sharna  (AFGL/LSI ) 

R.  McCoy  (NRL) 

W.  Sharp  (Univ.  Mich.) 

P.  Hays  (Univ.  Mich) 
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V.  EMPIRICAL  MODELS 


TIME:  1020-1140,  1310-1340 
-  CHAIR:  T.  L.  Killeen  (U.  of  Michigan) 

1.  MSIS  86  Empirical  Model,  Status 
and  Future  Directions 

2.  Model  Evaluations  from  AFGL 
Accelerometer  Data 

3.  CIRA  86,  80-120  KM 

LUNCH:  1125-1310 

4.  GRAM  Model 

5.  Session  Summary 


VI.  GENERAL  CIRCULATION  MODELS 
TIME:  1340-1505,  1525-1615 
-  CHAIR:  J.  M.  Forbes  (Boston  Univ.) 

1.  NCAR  TGCM 

2.  Operational  TGCM's 

3.  Global  Thermospheric  Specifications 
and  Forecasting 

4.  U.  Coll,  Lond.  TGCM 

BREAK:  1505-1525 

5 .  Future  TGCM  Needs 

6.  Coupled  Thermo-Ionosphere  TGCM 

7.  Session  Summary 

**  THURSDAY,  22  OCTOBER  1987  ** 


A.  Hedin  (NASA/GSFC ) 

F.  Marcos  (AFGL/LIS) 

G.  Groves  (Univ.  Coll.  Lond.) 


D.  Johnson/G.  Fichtl  ( NASA/MSFC ) 
T.  Killeen  (Univ.  Mich) 


R.  Roble  (NCAR) 

T.  Killeen  (Univ.  Mich) 
Capt.  C.  Tschan  (HQ  AWS) 

D.  Rees  (Univ.  Coll.  Lond) 


R.  Schunk  (Utah  St.  Univ.) 
D.  Anderson  (AFGL/LIS) 

J.  Forbes  (Boston  Univ.) 


VII.  DRAG  COEFFICIENTS 
TIME:  0800-0935 
-  CHAIR:  Blanchard  (NASA/LaRC) 

1.  Aerodynamic  Drag  Calculations  and  the 
Effects  of  Surface  Reflection 
Uncertainties 

2.  Direct  Simulation  Monte  Carlo  Drag 
Predictions 

3.  In-Orbit  Measurements  of  Atom  Scattering 
for  Drag  Calculations 

4.  Satellite  Drag  Coefficient  Variability 

BREAK:  0935-1150 


G.  Bird  (Univ.  of  Sydney) 

J.  Moss  (NASA/LaRC) 

J.  Gregory/G.  Karr  (NASA/MSFC) 

K.  Tobiska  (Univ.  of  Colo.) 
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VIII.  NEW  DENSITY  AND  DRAG  MEASUREMENTS/DIRECT  SENSING 
TIME:  0955-1150 

-  CHAIR:  H.  C.  Carlson  (AFGL/LIS) 

1.  Atmospheric  Density 

Specification  Satellite  (ADS) 

2.  ADS  Mass  Spec.,  Absolute  Density 

3.  ADS  Mass  Spec.,  Winds 

4.  ADS  FPI  and  Feasibility 

5.  Shuttle  Accel.  Density  Results 

6.  Session  Summary 

LUNCH:  1150-1300 


IX.  SOLAR  AND  flAGNETOSPHERIC  INPUTS 
TIME:  1300-1500 

-  CHAIR:  N.  C.  Maynard  (AFGL/PHG ) 

1 .  Solar  Flux 

2.  Solar  Flares  and  UV  Inputs 

3.  Gd-based  Solar  Coronal  Observations 
as  a  Possible  Input  to  Atmospheric 
Models 

4.  Future  Solar  and  Solar  Wind 
Measurements 

a.  GOES  Solar  X-ray  Imager 

b.  SAMEX 

c.  SIMPL 

d.  EUV  Calibration 

5.  Geomagnetic  Predictions 

6.  Magnetospheric  Convection 
and  DMSP 

7.  Particle  and  Joule  Energy 
Dissipation  from  DMSP 

8.  Evans  Index 

BREAK:  1500-1520 


X.  SUMMARY/CONCLUSIONS 
TIME:  1520-1620 
-  CHAIR:  F.  A.  Marcos  (AFGL ) 

1.  Absolute  Density  Chairmen 

2.  Winds,  Composition,  Temperature,  Drag  Chairmen 

ADJOURN:  1620 


R.  Donnelly  ( NOAA ) 

S.  Kei  1  (AFGL/PHS) 

R.  Altrock  (AFGL/PHS) 


W.  Wagner  (NOAA) 

S.  Kei 1  (AFGL/PHS) 

M.  Shea  (AFGL/PHP) 

W.  Parkinson  (SAO) 

J.  Joselyn  (NOAA) 

'N.  Maynard/F.  Rich  (AFGL/PHG) 

F.  Rich/S.  Gussenhoven/ 
(AFGL/PHP) 

D.  Evans  NOAA) 


F.  Marcos  (AFGL/LIS) 

J.  Ballenthin  (AFGL/LID) 

D.  Kayser  (Aerospace  Corp) 
T.  Killeen  (llniv.  Mich) 

R.  Blanchard  (NASA/LaRC) 

H.  Carlson  (AFGL/LIS) 
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MODELING  LIMITATIONS 
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AWS  SOLAR  OBSERVATORIES 
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SATELLITE-BASED  OBSERVATIONS 
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INTEGRATED  SPACE  ENVIRONMENTAL 
MODELING  SYSTEM  (ISEMSI 
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POINT  DENSITY  SPECIFICATION 

Crystal  L.B.  Sihaaj 
AFGL/LYA 


The  IMPROVED  POINT  ANALYSIS  MODEL  (1PAM)  is  a  new  Air  Force  Globa! 
Weather  Centra!  operational  model  that  describes  the  atmosphere  above  and  around  a 
given  point  on  the  earth’s  surface  for  a  specific  time.  This  new  model,  which  will  be 
operational  in  the  fall  of  1988,  will  represent  a  significant  improvement  in  the  quality  of 
point  analyses  and  an  expansion  in  the  delivery  of  customer  required  parameters. 
However,  there  are  still  a  number  of  areas  where  enhancements  will  be  needed.  AFGL, 
under  guidance  from  AWS,  is  undertaking  a  variety  of  research  efforts  aimed  at  further 
improving  IP  AM.  As  techniques,  algorithms,  or  models  are  developed,  emphasis  will  then 
be  placed  on  successfully  transitioning  the  technology  from  AFGL  researchers  to  AWS 
personnel  who  can  then  implement  the  ideas  at  Global  Weather  Central. 

IP AM  will  provide  its  customers  with  a  listing  of  the  atmospheric  conditions  at  and 
around  a  requested  point.  A  pseudo-surface  observation  will  be  deduced  from  any  and  all 
available  meteorological  surface  information.  A  cloud  depiction  will  be  provided, 
extracted  from  the  Real-Time  Nephanalysis  (Fye,  1978).  Information  will  be  produced 
describing  the  atmosphere  and  choosing  the  appropriate  default  inputs  needed  to  run 
LOWTRAN  (Kneizys,  et  ai,  1983  --  LOWTRAN  is  a  separate  model  that  calculates 
atmospheric  transmittance  and  radiance  over  a  given  path).  Finally  a  vertical  profile  of 
the  atmosphere  will  be  provided.  Below  100,000  ft,  this  profile  is  based  on  RAOB, 
ROCOB,  satellite  soundings,  aircraft  reports,  and  cloud-tracked  winds  and  is  produced 
using  optimum  interpolation  techniques.  Above  100,000  ft,  the  profile  merely  depends  on 
Groves  Model  (Groves,  1985)  linked  to  MSIS83  (Hedin,  1983)  in  a  fashion  similar  to  that 
described  by  Groves  (1987). 

There  is  a  requirement  to  improve  the  specification  of  temperature,  pressure,  and 
relative  humidity  to  120  km  and  density  to  200  km.  Both  vertical  and  line-of-sight  profiles 
are  needed.  Suggestions  and  techiques  discussed  at  this  workshop  will  be  used  to  enhance 
and  expand  the  capabilities  of  IP  AM. 


Fye,  F.K.,  1978;  The  AFGWC  automated  cloud  analysis  model.  AFGWC  Tech.  Memo. 
78-002,  97p. 

Groves,  G.V.,  1985;  A  global  reference  atmosphere  from  18  to  80  km.  Air  Force  Surveys 
in  Geophysics,  no.  448,  AFGL-TR-85-0129,  121p  ,  AD  AJ.62499. 

_ ,  1987;  Modeling  of  atmospheric  structure:  70-  130  km.  AFGL-TR-87-0226, 

AD  A201077. 

Hedin,  A.,  1983;  A  revised  thermospheric  model  based  on  mass  spectrometer  aud 
incoherent  scatter  data.  J.  Geophys.  Res.,  88,  p.  10170. 

Kneizys,  F.X.,  Shettle,  E.P.,  Gallery,  W.O.,  Chetwynd,  J.H.,  Abreu  L.W..  Selby,  LF  A  , 
Clough,  S.A.,  and  Fenn,  R.W.,  1983;  Atmospheric  transmittance/radiance 
computer  code  LOWTRAN  6.  Environmental  Research  Papers,  no.  8-10. 
AFGL-TR-83-0187,  200p  ,  AD  A137786. 
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POINT  PAPER 

on 

X-30  National  Aero-Space  Plane  (NASP)  Environmental  Support 

Lt.  Steven  Johnson 
DET  1,  2WS,  Wright  Patterson  AFB 

Presented  at  the  Atmospheric  Density  and  Aerodynamic  Drag  Models  for  Air  Force 
Operations  Workshop,  Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  20  Oct  87. 

DISCUSSION 

Program  organization  and  goal. 

Joint  DOD  (DARPA,  Air  Force,  Navy,  and  SDIO)  and  NASA  program. 

Joint  Program  Office  (JPO)  located  at  Wright-Patterson  AFB. 

-  As  of  1  Nov  Program  Manager  will  be  Dr  Robert  Barthelemy. 

-  1 Lt  Steve  Johnson  (Det  1,  2  WS)  provides  support  to  the  JPO. 

Program  goal  is  to  demonstrate  the  technologies  required  for  hypersonic 
cruise  and  single-stage-to-orbit  (SSTO)  capabilities. 

-  X-30  is  a  research  vehicle  not  being  designed  with  any  specific 

mission  or  with  mass  production  in  mind. 

Flight  envelope. 

Limited  time  in  the  troposphere. 

Stratosphere  is  the  primary  flight  regime. 

-  Hypersonic  cruise:  Mach  8-12  at  altitudes  100,000  to  150,000  ft. 

-  Acceleration  to  high  Mach  numbers  for  orbital  insertion  here. 

Above  50  km  NASP  will  fly  an  orbital  insertion  trajectory. 

Above  80  km  NASP  is  essentially  spacecraft. 

-  Powered  by  rocket  motor,  speeds  near  escape  velocity. 

-  Ascent  trajectory  very  different  from  Space  Shuttle. 

-  NASP  accelerates  within  the  atmosphere. 

-  Descent  trajectory  similar  to  Space  Shuttle. 

Ultimately  to  obtain  a  low  Earth  orbit. 
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Desired  environmental  data. 


The  most  current  reference  atmosphere  above  90  km. 

-  Standard  deviation  data  desired  for  flight  control  system  design. 

-  Applicability  of  data  under  varying  solar  conditions. 

Identify  spectrum  of  gravity  waves  to  be  encountered. 

-  An  as  yet  undefined  crucial  magnitude  to  vehicle  design  fully 

understood . 

-  Wave  sources  known  and  forecastable. 

A  characterization  of  turbulent  extremes  near  the  turbopau^c  is  needed. 

-  Effects  of  breaking  gravity  waves  quantified. 

Effects  of  solar  tides. 

-  Changes  in  density/drag  with  diurnal  and  semidiurnal  tides. 

-  Areal  extent  at  varying  latitudes. 

Test  flights  are  set  to  begin  in  1993. 

-  Sunspot  cycle  places  mid  1990's  in  period  of  increased  solar  activity. 

Conclusions. 

Phase  2  work  now  beginning. 

-  Engine  and  airframe  contracts  have  been  awarded. 

-  Emphasis  shifting  from  conceptual  design  to  technology  development. 

More  accurate  data  desired. 

-  Means  and  perturbations  both  important  to  vehicle  design. 

-  Atmospheric  models  applicable  over  a  wide  range  of  conditions. 

Test  flights  scheduled  for  the  mid  1990's. 

iLt  Steven  Johnson/Det  1 ,  2  WS/AUTOVON  785-2208 
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BRIEFER:  LT  STEVE  JOHNSON 


DOD /NASA  NATIONAL  AEROSPACE  PLANE 
PHASE  2  ORGANIZATION 
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SOURCES  KNOWN  AND  FORECASTABLE 


DESIRED  ENVIRONMENTAL  DATA 


PERIOD  OF  INCREASED  SOLAR  ACTIVITY 
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TEST  FLIGHTS  IN  THE  MID-1990S 


2.4 

ATMOSPHERIC  DENSITY  MODELS 
IN  OPERATIONAL  SYSTEMS 

BY 

JOSEPH  J  F  LIU 

Space  Division/DOA 

HQ  AFSPACECOM 


VGs  PREPARED  FOR 

Workshop  on  Atmospheric  Density  and  Aerodynamic  Drag  Models  for 
A i r  Force  Operations 

20-22  October  1987 

AFGL 
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OPS  REQUIREMENTS  ON  SATELLITE  DYNAMICS 

-  Accuracy  of  Neutral  Density 

-  Parameter  Predictability 
Short -term/ Long -term 

-  Altitude  Range 

-  Computational  Efficiency 

-  Improvements  Over  Existing  Model 

-  Compatibility  Issues 

-  Additional  Features : 

Wind,  Cmhrge  "drag7. 
New/additional  Parameters,  etc- 
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ATMOSPHERIC  DENSITY  MODELS  IN  USE 


Aiialytic  Models 

-  Exponential  Models 

-  Basic  Model 

-  Modified  and  Improved  Models 

-  Power  Function  Models 

-  Basic  Model 

-  Modified  and  Improved  Models 
Em? i rical  Models 

-  DENSEL 

-  Harris-Priester 

-  Jacchia  Models  and  Derivatives 

-  Jacchia  64,  70,  71,  77 

-  Jacchia-Bass  11 

-  Walker  64 

-  Lafontaine's  Analytic  Model 

-  Earlier  Model 

-  Mass  Spectrometer  and  Incoherent 
Scatter  Data  Models 
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OPS  PROBLEMS  RELATING  TO  DENSITY  MODELS 


tlo  Significant  Improvements  on  Neutral  Density 
No  Reliable  Parameter  Prediction  Capability 
No  Improvement  on  Efficiency 
Proliferation  &  Compatibility  Problems 
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PROJECTS  REQUIRE  BETTER  ATM  MODEL 


-  TIP 

-  COMBO 

-  Laser  Clearinghouse 

-  Space  Shuttle 

-  Sensor  Calibration 

-  Military  Applications 

-  Others 
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SPADOC  EPHEMERIS  MODELS  USI1J6  ATM  DENSITY  MODEL 


-  General  Perturbations  Theories.-  SGPA,  HAN DE 

-  Semi -analvt ic  Theory:  SALT 

-  Special  Perturbations  Theory:  SP 

IJote:  across-the-board  effects  on  all  ops  applications 
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OPS  £  SCIENCES  COOPERATION 


Initiate  OPS  £  Sciences  Working  Group 

-  Time:  week  of  15  Feb  1988 

-  Location:  AFSPACECOM 

Identify  Mutual  Needs 

-  Emphasize  Team  work  2  Compatibility 
Set  OPS  Requirements  and  RsD  Goals 

-  Short-term  £  Long-term 
Assess  Impacts  £  Priorities 
Consolidate  Resources 
Establish  Feasible  Action  Items 
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NUMERICAL  EXAMPLE  1 


Orbital  Characteristics 

SEMI “MAJOR  AXIS  =  7000  KM,  ECCENTRICITY  =  -0667 
INCLINATION  =  50  DEG,  EPOCH  ALTITUDE  =650  KM 
BALLISTIC  COEF  =  *00213  m**2/kg,  epoch:  1970  136 


Days  Since 
Epoch 


Propagation  Errors  (km) 

10%  15% 


NUMERICAL  EXAMPLE  2 


Orbital  Characteristics 


SEMI -MAJOR  AXIS  =  6718  KH, 
INCH  HAT  I  ON  =  50  DEG, 
BALLISTIC  COEF  *  -0018, 


ECCENTRICITY  =  -0293 
EPOCH  ALTITUDE  =  466  KM 

epoch:  1970  220 


Days  Since  Propagation  Errors  (km) 


5  % 

10% 

15% 

20% 

1 

11 

22 

33 

43 

2 

46 

91 

137 

182 

5 

280 

559 

839 

1120 

10 

1261 

2516 

3750 

5000 

* 


2-37 


NUMERICAL  EXAMPLE  3 


Orbital  Characteristics 

SEMI-MAJOR  AXIS  =  6600  KM,  ECCENTRICITY  =  -0012 
INCLINATION  =  32  DEG,  EPOCH  ALTITUDE  =  217  KM 
BALLISTIC  COE F  =  *0286,  EPOCH:  1978  258 


Hrs  Since  Propagation  Errors  (km) 

Epoch 


5% 

10% 

15% 

20% 

12 

20 

A0 

60 

80 

2  A 

108 

215 

320 

A30 

32 

260 

530 

80  0 

1100 

38 

550 

1000 

im 

1750 
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RECOMMENDED  REQUIREMENTS 


AlTITUjE 

Average  Neutral 
Density  Error 

<  200  KM 

5 1 

500-700  km 

10% 

>  700  KM 

15% 

2-39 


SLIDE  1 


SPACE  OPERATIONS  WORKSHOP  III 
AFSPACECOM  DO  PANEL:  TOPIC  4  -  ASTRO DYNAMIC  STANDARDS 


PLAN  FOR  ASTRODYNAMIC  STANDARDS 

1.  Assess  the  Need  for  Astrodynamic  Standards 

2.  Identify  the  Areas  that  Require  Standards 

3.  Quantify  Operational  Impact 

4.  Prioritize  Standard  Elements 

5.  Develop  Recommendations 
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SLIDE  2 


SPACE  OPERATIONS  WORKSHOP  III 
AFSPACECOM  DO  PANEL:  TOPIC  4  -  ASTRODYNAMIC  STANDARDS 

JUSTIFICATIONS  FOR  ASTRODYNAMIC  STANDARDS 

1.  Ensure  the  Integrity  of  Operational  Systems 

2.  Optimize  the  Communication  between  Centers  and  Users 

3.  Provide  Baseline  for  Future  Systems 

4.  Enhance  Compatibility  and  Inter-operability 
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3.  PHYSICS  OF  DENSITY  AND  DRAG  VARIABILITY 


3  -  1 


3.1 


MESOSPHERE-THERMOSPHERE 

COUPLING 


Jeffrey  M.  Forbes 

Department  of  Electrical,  Computer,  and  Systems 

Engineering 

Boston  University,  Boston,  MA  02215 


Workshop  on 

Atmospheric  Density  and  Aerodynamic  Drag 
Models  for  Air  Force  Operations 
Air  Force  Geophysics  Laboratory 
20  -  22  October  1987 
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1.  MESOPAU SE-REG ION  PROCESSES 


The  mesopause  is  the  minimum  in  the  earth's  atmospheric  temperature 
vertical  profile  near  90  km  (Figure  1).  and  as  such  represents  the  bound¬ 
ary  between  the  mesosphere  (a  region  of  negative  temperature  gradient, 
ca.  50  -  90  km)  and  the  thermosphere  (a  region  of  positive  temperature 
gradient,  ca.  90  -  500  km).  A  schematic  summary  of  the  important  in¬ 
teractions  which  take  place  in  the  height  region  around  the  mesopause  is 
provided  in  Figure  2.  The  primary  sources  of  heat  input  into  the  lower  ther¬ 
mosphere  include  absorption  of  EUV  and  UV  radiation,  Joule  dissipation 
of  ionospheric  currents,  particle  precipitation,  and  dissipation  of  wave  and 
turbulent  energy.  The  primary  heat  sink  between  100  and  140  km  is  down¬ 
ward  heat  conduction.  Below  100  km  the  densities  and  collision  frequencies 
of  polyatomic  molecules  C02,  O 3.  H20 .  and  NO  are  sufficiently  great  that 
heat  loss  by  infrared  radiation  plays  a  significant  role  in  determining  the 
thermal  structure.  Differential  heating  in  latitude  tends  to  set  up  a  large  - 
scale  (geostrophic)  circulation,  as  it  does  in  the  stratosphere.  However,  up¬ 
ward  propagating  gravity  waves  can  achieve  convective  instability  at  these 
heights,  leading  to  a  momentum  flux  divergence  which  effectively  causes 
a  force  imbalance,  sometimes  called  a  ’wave  drag’,  in  the  mean  eastward 
momentum  equation.  The  force  imbalance  induces  a  mean  meridional  cir¬ 
culation,  and  thereby  a  Coriolis  torque  opposing  the  drag  and  adiabatic 
heating/ cooling,  giving  rise  to  a  reversed  meridional  temperature  gradient 
in  the  region  of  the  mesopause  (Figure  1  inset). 

Gravity  waves  are  thought  to  originate  in  such  tropospheric  sources  as 
storms  and  topographic  excitation.  In  addition  to  depositing  momentum, 
dissipating  gravity  waves  induce  turbulent  mixing  which  strongly  influences 
the  transport  and  distributions  of  long-lived  chemical  species  such  as  NO. 
Molecular  diffusion,  as  its  efficiency  increases  with  decreasing  density,  usu¬ 
ally  becomes  faster  than  turbulent  diffusion  rates  above  about  100  km,  the 
level  sometimes  referred  to  as  the  ’turbopause’.  It  is  well  known  that  the 
level  of  the  turbopause  is  a  major  factor  controlling  the  composition  and 
hence  the  total  mass  density  in  the  thermosphere. 

Similarly,  tidal  oscillations  of  24-hour  (‘diurnal’)  and  12-hour  (’semidi¬ 
urnal)  periods  are  generated  by  H20  and  O3  insolation  absorption  in  the 
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troposphere  and  stratosphere  and  propagate  through  the  mesopause  into 
the  lower  thermosphere  (Figure  3).  Tides  represent  a  major  source  of  tem¬ 
perature  and  wind  variability  in  the  mesopause  region,  and  are  also  capable 
of  generating  turbulence  and  depositing  mean  momentum  and  heat  above 
85  km.  In  addition,  interactions  with  the  magnetoplasma  at  higher  levels 
induce  important  redistributions  of  ionization  both  directly  and  through 
the  electric  fields  generated  via  dynamo  action.  These  effects  feed  back 
via  the  ion  drag  mechanism  to  affect  the  neutral  dynamics  and  density 
variations. 

Within  this  workshop  session  the  subject  of  gravity  waves  is  addressed 
by  Dr.  R.L.  Walterscheid,  and  that  of  the  penetration  of  tidal  waves  above 
100  km  and  their  consequences  above  about  120  km,  by  Dr.  C.G.  Fesen. 
This  paper,  therefore,  will  concentrate  on  tidal  oscillations  between  80  and 
120  km  and  their  contributions  to  density  variability  within  this  height 
regime. 


2.  THE  STATUS  OF  ATMOSPHERIC  TIDES 
RESEARCH 

Atmospheric  tides  represent  one  obvious  mechanism  for  dynamically  cou¬ 
pling  the  lower  atmosphere  with  the  upper  atmosphere  and  ionosphere. 
Semidiurnal  and  diurnal  tidal  oscillations  are  excited  by  i/20  and  O 3  in¬ 
solation  absorption  in  the  troposphere,  stratosphere,  and  mesosphere,  and 
propagate  upwards  into  the  lower  and  upper  thermosphere.  In  the  meso¬ 
sphere  mode  coupling  occurs  via  interactions  with  the  mean  wind  field  to 
modify  the  mixture  of  modes  driven  directly  by  thermal  excitation.  The 
shorter  wavelength  modes  dissipate  and  deposit  mean  momentum  and  heat 
at  the  base  of  the  thermosphere.  The  longer  wavelength  modes  penetrate 
to  upper  levels  and  combine  with  oscillations  excited  in-situ  to  produce 
the  total  observed  diurnal  variation  of  winds,  temperatures,  densities,  and 
pressures  at  any  given  height  and  latitude. 

Quantitatively,  the  degree  to  which  tides  of  low-er  atmospheric  origin 
affect  the  diurnal  and  latitudinal  structure  of  the  thermosphere  as  a  func¬ 
tion  of  height,  season,  and  solar  cycle  is  still  not  well  known.  Certainly,  the 
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effects  are  expected  to  become  more  prominent  as  the  EUV  excitation  di¬ 
minishes  with  declining  solar  activity.  Insofar  as  modelling  is  concerned,  a 
major  uncertainty  iti  thermospheric  general  circulation  models  (TGCMs)  is 
specification  of  tidal-period  oscillations  at  the  lower  boundary  near  95  km. 
This  can  be  greatly  aided  by  global  observations  available  from  MST,  me¬ 
teor  wind,  and  medium-frequency  radars  (Figures  4  and  5),  combined  with 
spectral  model  results  for  the  height  regime  below  100  km.  The  spectral 
models  as  they  now  exist  must  be  upgraded  to  incorporate  more  realistic 
specifications  of  heating  rates,  turbulent  diffusivities,  and  mean  winds  on  at 
least  a  monthly  basis.  Modelling  and  observational  efforts  must  be  directed 
both  along  the  lines  of  ’climatology’  and  ’weather’  ;  that  is,  while  it  is  im¬ 
portant  to  understand  seasonal  and  global  trends,  other  efforts  must  focus 
on  transient  phenomena  such  as  equinox  transitions  or  stratwarm  effects. 

In  addition,  further  work  is  required  to  understand  the  physics  of  gravity 
wave/tidal  interactions,  related  implications  with  regard  to  the  interpreta¬ 
tion  of  observations,  and  how  such  effects  can  be  parameterized  in  models 
which  simulate  tides.  For  instance,  w'ork  by  D.  Fritts,  R.  Vincent,  and  R. 
Walterscheid  show's  that  tides  can  induce  a  daily  variation  in  the  momen¬ 
tum,  heat,  and  turbulence  deposition  by  dissipating  gravity  wraves.  The 
accompanying  wind,  temperature,  and  density  fluctuations,  which  at  this 
point  in  time  can  only  be  speculated  on,  would  be  expected  to  occur  at 
higher  levels  as  wrell  as  in-sit u. 


3.  DENSITY  VARIATIONS  DUE  TO  TIDES 


Numerical  models  typically  solve  a  reduced  set  of  equations,  either  cast  in 
terms  of  the  perturbation  geopotential  .  or  the  three  velocity  components 
and  temperature,  depending  on  how  diffusion  is  parameterized.  Assuming 
the  tidal  temperature  perturbations  bT  are  small  compared  with  the  mean 
background  temperature  T0,  one  can  derive  an  approximate  expression  for 
the  relative  variation  in  total  mass  density: 


bp 

Po 


bT  \ 

T,  ( 


e.ip{ 


}  I 


where 
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Using  the  numerical  simulations  by  Forbes  (J.  Geophys.  Res.,  87,  p. 
5222  -  5252,  1982),  tidal  variations  in  density  as  illustrated  in  Figures  6  and 
7  can  be  computed.  The  estimated  variations  of  5  -  20  percent  between  80 
and  120  km  ought  to  be  important  for  aerospace  operations  in  this  height 
regime. 


4.  THE  LOWER  THERMOSPHERE  COU¬ 
PLING  STUDY 

The  subject  of  mesosphere-thermosphere  coupling  is  the  focus  for  the  Lower 
Thermosphere  Coupling  Study  (LTCS)  of  the  CEDAR  (  Coupled  Energetics 
and  Dynamics  of  the  Atmospheric  Regions  )  and  WITS  (  World  Ionosphere 
-  Thermosphere  Study)  Programs.  The  LTCS  (Figure  8),  which  began  in 
1986,  is  a  coordinated  investigation  of  the  lower  thermosphere  (  80  -  180  km) 
combining  observational  and  numerical  modelling  efforts  with  the  ultimate 
goal  of  better  understanding  the  dynamic  and  electrodynamic  processes 
coupling  the  mesosphere,  lower  thermosphere,  and  upper  thermosphere  re¬ 
gions.  Primary  emphasis  is  placed  on  neutral  winds  and  temperatures,  and 
electric  fields  and  currents,  as  these  are  among  the  parameters  of  main  sci¬ 
entific  interest,  and  there  is  a  need  to  keep  the  study  as  narrowly  focused 
as  possible.  Nevertheless,  the  better  understanding  of  physical  phenomena 
and  improved  modelling  capabilities  likely  to  emerge  from  the  study  will 
enhance  the  capability  to  specify  total  mass  density  as  well. 
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TIDAL  LOWER  BOUNDARY  CONDITION 
(  PERTURBATION  GEOPOTENTIAL  )  FROM  WIND  DATA 


Momentum  equation  from  classical  tidal  theory: 


du 

dt 


1 

acosO  d\ 


where 

u  =  eastward  velocity 
v  =  northward  velocity 
6  =  latitude 
A  =  longitude 
£  =  time 

/  =  coriolis  parameter  =  20  sin  0 
9  =  acceleration  due  to  gravity 
a  =  radius  of  Earth 
0  =  27r/24  hours _1 


introducing  , 


tun'  —  fv'  = 


is 


acosB 


& 


For  migrating  tide  u;  =  sO  ;  therefore: 

,,  „  .  f  ,  ,2sin0  ,1 

$  =  — oOcos  0  <  i/  +  j - v'  > 

[see  also  Portnyagin  and  Kajdalov,  1986  ] 


(2) 


(3) 
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From  Classical  Tidal  Theory: 


$'  =  £>„©„((?) 
n 

(4) 

u>  = 

n 

(5) 

(6) 

n 


where 

0n(0)  =  Hough  function  (orthogonal) 

Un(0)  =  eastward  velocity  expansion  function 
Vn(6)  —  northward  velocity  expansion  function 

The  <f>n  s  can  be  determined  from  radar  wind  data  by  least- 
squares  fitting 

(a)  <&'  [  Equations  (3)  and  (4)  ] 

l*>}  u‘  and/or  v‘  directly  [  Equations  (5)  and/or  (6)  ] 


DEGREES  LATITUDE 


degrees  latitude 


PERCENT  DENSITY  DEVIATIONS 

EQUINOX  z  =  120  km 


Figure  -  Estimated  percent  density  derivations  from  mean 
equinox  due  to  atmospheric  tides. 
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at 


THE 

LOWER  THERMOSPHERE  COUPLING  STUDY  (LTCS) 

of  the 

CEDAR  and  WITS  Programs 


PURPOSE  AND  SCOPE 

The  Lower  Thermosphere  Coupling  Study  (LTCS)  is  a  co¬ 
ordinated  investigation  of  the  lower  thermosphere  (80  -  180 
km)  combining  observational  and  numerical  modelling  efforts 
with  the  ultimate  goal  of  better  understanding  the  dynamic 
and  electrodynamic  processes  coupling  the  mesosphere,  lower 
thermosphere,  and  upper  thermosphere  regions.  Primary  em¬ 
phasis  will  be  placed  on  neutral  winds  and  temperatures,  and 
electric  fields  and  currents. 
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ATMOSPHERIC  GRAVITY  WAVES 
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DENS 
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CCELERATION  DUE  TO  GRAVITY 


Period/hr 

100  24  10  2  0  5  01 


Seasonally  averaged  spectral  densities  for  horizontal 
wind  fluctuations  observed  at  a  height  of  85  km  at  Adelaide 
and  Townsville.  The  spectra  are  the  average  of  the  zonal  and 
meridional  spectra  and  the  bars  indicate  the  range  of  the 
indivial  spectra.  A  Show  spectral  densities  of  wind  motions 
observed  at  Poker  Flat  (after  Carter  and  Balsley,  1982). 
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THERMOSPHERE 
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Thule  (  77n.  69h  i  primrose  lrke  (  sun.  now  i 


HI  ROTA  (1984) 


Wave  Distribution 
4489  Wave  Samples 
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Magnetic  Local  Time(hr) 

The  wave  occurrence  histogram  shows  the  distribution 
of  waves  of  all  amplitudes  as  a  function  of  magnetic  local  time. 
The  greatest  number  of  waves  are  found  in  the  midnight  sector, 
although  secondary  maxima  near  noon  indicate  the  possible 
transpolar  propagation  of  waves  assumed  to  be  generated  near 
midnight. 
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[_  Averoged  Mognetic 
3  hr  Activity  Index 
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Geomognetic 
lot.  E.  long. 
Goose  Boy  64.6  12.1 

Lorirtg  58.4  2.1 
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This  figure  shows  a  sum. nary  of  wave  occurrences  dur¬ 
ing  a  magnetic  storm  in  February  1974.  The  top  and  bottom 
panels  show  the  N,  number  density  at  330  km  and  wave  occur¬ 
rence  and  amplitude  at  200  km,  respectively.  The  center  panel 
shows  an  averaged  magnetic  activity  index  for  two  high-latitude 
stations.  The  vertical  arrows  show  the  closest  approach  to  the 
geomagnetic  pole.  The  largest  waves  were  seen  near  noon  in 
regions  of  lower  IN,!  and  appear  to  be  generated  near  midnight, 
as  is  indicated  by  the  increases  in  magnetic  activity  at  these 
times. 
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SUMMARY 
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3.3 


Tidal  Effects  on  Thermospheric  Circulation,  Densities,  and  Temperatures 

C.  G.  Fesen 

Before  discussing  the  influence  of  tides  on  the  thermosphere,  a  brief  review  of  atmospheric 
tides  is  presented.  The  review  attempts  to  answer  the  most  basic  questions  usually  encountered 
in  tidal  research,  and  are  listed  in  Figure  1.  These  questions  include:  What  are  tides?  Why  are 
they  important?  Where  do  they  come  from?  The  incorporation  of  tidal  effects  in  thermospheric 
models  follows,  along  with  a  brief  summary  of  model  results.  Finally,  suggestions  are  made 
regarding  what  can  be  done  to  improve  the  thermospheric  tidal  research  now  being  done. 

Tides  are  defined  to  be  global  scale  oscillations  of  atmospheric  quantities  with  periods 
that  are  integral  multiples  of  a  solar  (or  lunar)  day.  Observations  which  show  variations  with 
a  G  hr,  8  hr,  12  hr,  or  24  hr  period  are  exhibiting  the  effects  of  a  tidal  perturbation. 

Figure  2  shows  an  example  of  tidal  effects  seen  in  neutral  wind  observations  from  Arecibo. 
The  top  two  panels  are  plots  of  the  meridional  winds,  and  the  bottom  panel  the  vertical  winds, 
as  a  function  of  height  and  local  time  for  consecutive  days.  The  shading  indicates  southward 
or  downward  winds.  There  is  an  obvious  oscillation  in  the  winds,  with  about  a  12  hr  period, 
or  a  semidiurnal  tide.  Higher  in  the  atmosphere,  this  changes  to  a  diurnal  tide. 

Similar  behavior  would  be  found  in  observations  of  the  temperatures  and  the  densities. 
Further,  these  oscillations  are  seen  at  other  observat ories  around  the  globe,  differing  in  the 
degree  of  variability,  but  still  exhibiting  time-dependent  behavior. 

Figure  2  therefore  illustrates  why  tidal  effects  are  important  in  thermospheric  studies: 
tides  are  very  evident,  very  strong  in  the  upper  atmosphere.  In  fact,  in  the  lower  thermosphere, 
tides  dominate  the  thermospheric  response;  i.e.,  in  analysis  of  wind,  temperature,  and  density 
observations,  the  semidiurnal  component  typically  overwhelms  the  mean  component. 

Since  tides  are  time-dependent,  they  arise  due  to  time-dependent  forcings.  There  are 
three  sources  in  the  thermosphere,  as  listed  in  Figure  3.  The  sources  due  to  solar  heating 
and  ion-neutral  momentum  coupling  are  relatively  easy  to  include  in  numerical  models  of  the 
thermosphere.  The  source  due  to  the  upward  propagation  of  energy  from  the  lower  atmosphere 
is  more  difficult  to  parametrize;  this  is  the  only  tidal  source  discussed  in  detail  here. 

The  solar  energy  absorbed  in  the  troposphere  end  stratosphere  generates  a  tidal  response 
in  these  atmospheric  regions.  Since  tides  are  a  form  of  internal  gravity  wave,  they  can  and  do 
propagate  vertically,  growing  exponentially  in  amplitude  with  altitude,  and  carrying  energy 
with  them.  Figure  4  shows  a  model  calculation  from  Walterscheid,  in  which  he  calculated  the 
semidiurnal  zonal  wind  from  the  ground  up  to  100  km,  resulting  from  the  absorption  of  solar 
energy  in  the  lower  atmosphere.  The  three  different  curves  are  for  three  different  model  runs; 
the  most  accurate  result  is  the  solid  curve.  A  wavelike  behavior  is  evident,  with  amplitudes 
growing  with  height  which  implies  that  the  phases  increase  downwards.  Obviously,  by  the  time 
the  wave  reaches  the  thermosphere  it  is  quite  large.  This  is  part  of  the  explanation  for  the 
strong  semidiurnal  wave  seen  in  the  lower  thermosphere  in  the  Arecibo  observations  in  figure 
2. 
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Thermospheric  models  must  be  able  to  include  the  effects  of  the  tides  that  are  generated 
in  the  lower  atmosphere  and  subsequently  penetrate  into  the  thermosphere.  Models  which 
extend  from  the  ground  through  the  thermosphere  are  not  possible  at  present.  A  simpler, 
do-able  scheme  is  based  on  the  following  approach:  if  a  model  such  as  that  shown  in  figure  4  is 
correctly  formulated,  it  should  be  able  to  reproduce  the  model  results  at  upper  levels  by  taking 
the  model  solutions  at  some  lower  level,  using  them  as  new  lower  boundary  conditions,  and 
running  the  model  from  this  new  lower  boundary.  If  correctly  formulated,  the  model  should 
reproduce  the  solutions  from  the  original  ’’full"  model  for  those  altitudes  above  the  new  lower 
boundary. 

This  is  the  approach  used  in  incorporating  lower  atmosphere  tidal  effects  in  the  NCAR 
TGCM.  A  distribution  of  the  model  variables  is  prescribed  at  the  lower  boundary  which 
is  wavelike.  The  wavelike  distribution  represents  the  waves  that  were  generated  in  the  lower 
atmosphere,  and  travelled  up  to  the  thermosphere.  An  alternate  description  is  that  the  wavelike 
distribution  is  analogous  to  shaking  the  lower  boundary,  like  shaking  a  rug:  this  will  cause 
oscillations,  which  will  travel  upwards  as  described  by  the  model  equations. 

The  wave  description  at  the  lower  boundary  of  the  NCAR  TGCM  is  in  terms  of  Hough 
modes.  The  Hough  modes  represent  part  of  the  solution  to  a  simplified  set  of  tidal  equations 
which  are  approximately  valid  near  and  below  100  km.  The  Hough  modes  describe  the  lat¬ 
itudinal  distribution  of  a  particular  wave,  and  are  functions  of  Legendre  polynomials  which 
are  well  known  mathematically.  Figure  5  shows  the  semidiurnal  Hough  modes.  The  modes 
are  identified  as  the  (2,2),  (2,3),  (2,4)  and  (2,5).  The  first  subscript  is  a  zonal  wavenumber, 
denoting  the  number  of  waves  that  fit  around  a  latitude  circle.  Two  waves  represent  a  semid¬ 
iurnal  tide.  The  second  index  is  the  meridional  wavenumber,  which  measures  the  degree  of 
waviness  or  structure  in  the  meridional  direction.  Subtraction  of  these  two  numbers  (the  zonal 
and  meridional  wavenumbers)  gives  the  number  of  nodes,  or  zero-crossings  of  t lie  equator,  for 
each  of  these  modes.  Thus,  the  2,2  has  no  nodes  and  is  symmetric  about  the  equator,  the  2,3 
has  one  node  at  the  equator,  and  is  asymmetric,  etc. 

Figure  5  also  lists  the  vertical  wavelengtlis  of  these  modes,  which  increase  with  meridional 
wavenumber.  Thus,  the  higher  order  modes  do  not  penetrate  far  into  the  thermosphere,  since 
they  are  more  likely  to  suffer  interference.  Consequently,  the  lower  thermosphere  exhibits  more 
wave-like  structure  and  variability,  which  gradually  damps  out  higher  in  the  thermosphere, 
until  only  the  lowest  order  modes  such  as  the  2,2  and  the  2,3  are  contributing  to  the  response. 

The  Hough  modes,  then,  are  the  functions  used  to  describe  the  latitudinal  structure  of  the 
model  variables —  pressure,  temperature,  winds,  and  density  —  at  the  lower  boundary  of  the 
NCAR  TGCM  which  is  at  97  km.  An  oscillatory  time  behavior  is  also  folded  into  the  lower 
boundary  description.  The  model  solutions  are  checked  by  comparison  with  data  for  similar 
conditions. 

Figure  G  presents  a  brief  summary  of  the  steps  followed  in  producing  a  model  run  of  the 
TGCM  which  includes  the  effects  of  upward  propagating  tides.  The  dashed  line  in  the  figure 
from  step  5  to  step  3  indicates  that  if  the  model  winds  and  temperatures  do  not  agree  well 
with  observational  data,  a  new  lower  boundary  condition  is  formulated,  and  the  model  run  is 
repeated  until  suitable  agreement  is  obtained  between  the  results  and  the  data. 
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To  illustrate  the  importance  of  including  upward  propagating  tides  in  theoretical  models, 
Figure  7  presents  TGCM  calculations  for  solstice  solar  cycle  minimum  which  are  identical 
except  for  the  inclusion  of  the  upward  propagating  tides.  The  top  (bottom)  panels  represent  the 
upper  (lower)  thermosphere,  near  approximately  300  (120)  km.  The  solid  contours  represent 
perturbation  temperatures,  and  the  arrows  wind  vectors,  with  the  magnitude  of  the  maximum 
wind  arrows  for  each  panel  shown  in  the  center  of  the  figure.  The  left  (right)  panels  are  the 
results  without  (with)  upward  propagating  tides. 

The  lower  thermosphere  shows  very  obvious  differences  in  the  two  results.  There  is  a 
definite  semidiurnal  character  to  the  temperatures,  reflecting  the  fact  that  the  semidiurnal 
tides  predominate  at  altitudes  of  100-150  km  or  so.  The  winds  obviously  will  change  to  reflect 
this,  as  will  the  densities.  There  are  significant  meridional  flows  at  low  latitudes,  unlike  the 
results  without  tides.  The  maximum  wind  speed  is  seen  to  have  doubled  with  tides.  The 
amplitude  of  the  temperature  variation,  from  maximum  to  minimum,  is  also  much  larger  with 
the  propagating  tides. 

The  upper  thermosphere  also  shows  significant  differences,  especially  near  midnight.  Note 
that  the  wind  speeds  again  are  stronger  than  they  were  without  the  tides,  and  that  the  ampli¬ 
tude  of  the  temperature  variation  is  larger,  which  will  translate  into  larger  amplitude  density 
variations. 

Figure  8  shows  that  the  inclusion  of  tides  in  calculations  is  important  and  necessary  to  pro¬ 
duce  realistic  simulation  of  observations.  The  neutral  wind  observations  from  Arecibo,  shown 
in  Figure  2,  are  reproduced,  along  with  model  calculations  of  the  wind  at  Arecibo's  location, 
without  and  with  the  tides.  The  improvement  is  quite  obvious.  The  lower  thermosphere  oscil¬ 
latory  behavior  is  not  reproduced  at  all  without  accounting  for  the  upward  propagating  energy 
in  the  model. 

Figure  9  shows  the  neutral  mass  densities  associated  with  the  wind  and  temperature 
results  of  figure  7,  with  the  same  format.  Again,  the  lower  thermosphere  calculations  show 
a  tremendous  difference  in  the  densities  calculated  with  and  without  upward  propagating 
tides.  The  densities  without  tides  show  a  diurnal  variation  that  is  not  realistic.  The  densities 
calculated  with  tides  show  strong  latitude  and  local  time  gradients. 

The  upper  thermosphere  shows  a  simpler  structure,  predominantly  diurnal,  but  with 
stronger  local  time  variations  at  most  latitudes  in  the  tides  case.  It  is  obvious  that  the  am¬ 
plitude  of  the  density  variation  is  obviously  much  greater  with  tides  than  without  throughout 
the  thermosphere. 

Seasonal  differences  are  illustrated  in  figure  10,  which  shows  solstice  and  equinox  tidal 
calculations.  Obviously,  the  equinox  distributions  are  more  symmetric,  since  the  forcing  is 
symmetric:  the  sun  is  on  the  equator.  The  solstice  results  are  asymmetric  because  the  forcing 
is  asymmetric:  the  sun  is  in  the  southern  hemisphere.  Again,  the  tides  are  seen  to  be  very 
important  in  the  lower  thermosphere.  There  are  much  stronger  gradients  and  more  structure  at 
solstice  than  equinox,  reflecting  the  fact  that  the  solstice  tides  have  contributions  from  all  the 
semidiurnal  modes,  both  symmetric  and  asymmetric.  The  amplitude  variation  of  the  density 
is  larger  during  solstice  than  equinox.  In  the  upper  thermosphere,  the  local  time  variation  at 
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equinox  is  relatively  flat  during  the  day  at  most  latitudes,  unlike  the  solstice  case. 

So  far,  then,  the  tides  have  shown  altitudinal,  latitudinal,  and  seasonal  dependences. 
Figure  11  shows  the  solar  cycle  dependence  of  the  tides.  Illustrated  in  the  top  panels  are  the 
calculated  temperatures  near  Arecibo  for  solar  cycle  minimum  and  maximum  as  a  function 
of  altitude  and  time;  the  bottom  panels  show  the  temperatures  and  winds  at  fixed  pressure 
levels  near  120  and  300  km.  It  is  evident  that  the  tides  remain  important  to  higher  altitudes 
near  solar  cycle  minimum  as  compared  to  solar  maximum,  which  is  due  to  the  fact  that  the 
density  increases  at  a  fixed  pressure  level  with  solar  activity.  This  means  that  the  atmosphere 
is  thicker,  denser,  and  consequently  more  viscous  during  solar  cycle  maximum,  so  that  the 
tides  coming  up  from  the  lower  atmosphere  are  damped  more  rapidly  than  during  low  solar 
activity. 

The  results  thus  far  have  been  mainly  theoretical;  figure  1?  shows  some  comparisons 
between  the  model  results  and  observations  to  illustrate  how  well  the  model  reproduces  the 
observed  thermospheric  response.  The  data  in  figure  12  are  on  the  right-hand  side,  and  were 
taken  by  the  Atmosphere  Explorer  satellite  'E’  (AE-E).  The  data  are  from  solstice  periods 
from  197C  through  197S.  The  satellite  was  in  a  circular  orbit  near  275  km  at  this  time,  and  was 
restricted  to  latitudes  between  ~  ±  20°.  The  satellite  observations  of  temperature  and  number 
densities  during  this  period  were  averaged  together  to  determine  average  temperatures,  and 
average  atomic  oxygen  and  N2  number  densities.  These  averages  were  then  decomposed  into 
mean,  diurnal,  and  semidiurnal  components.  This  work  was  done  by  USAF  Lt.  Mark  Storz, 
as  part  of  his  master’s  thesis  at  the  University  of  Michigan  in  the  summer  of  19S7. 

Only  semidiurnal  results  are  plotted  in  figure  12,  with  the  corresponding  results  from  the 
TGCM  shown  on  the  left  of  the  figure.  The  highs  and  lows  show  the  maxima  and  minima  of 
the  semidiurnal  temperature  and  density  variations;  the  amplitudes  are  double  either  of  these 
numbers.  The  solid  lines  running  through  the  highs  represent  the  phases,  or  times  of  maxima, 
which  are  read  off  the  bottom  axis. 

The  agreement  is  quite  encouraging.  Granted,  at  first  glance,  the  model  results  are  skewed 
toward  the  north,  while  the  data  show  symmetry  about  the  equator.  But  the  numbers  and  the 
morphology  show  decent  agreement  for  a  first  cut.  The  density  amplitudes  are  quite  similar, 
while  the  phase  agreement  is  fair,  with  the  model  phases  about  2  hours  early.  The  temperature 
amplitudes  are  50%  too  high  in  the  model,  while  the  phase  is  4  hours  early.  Why,  then,  is  the 
agreement  claimed  to  be  ’’encouraging”?  This  represents  the  first  attempt  to  match  t Ire  model 
output  to  this  type  of  extensive  data  coverage.  There  was  no  tuning  of  the  model  to  these 
conditions,  beyond  matching  the  season,  solar  cycle,  and  geomagnetic  conditions.  And.  due 
to  certain  parametrizations  in  the  model,  the  low  latitude  region  is  a  difficult  region  to  model 
accurately.  Further,  data  coverage  is  very  limited  in  this  region.  The  comparison,  then,  seems 
to  show  that  the  model  is  on  the  right  track,  and  progress  can  be  made  in  understanding  and 
alleviating  the  differences  between  the  two.  For  example,  the  model  run  had  the  amplitudes  of 
the  2,2  and  2,3  modes  the  same  at  the  lower  boundary.  As  a  reminder,  the  2,2  and  2.3  latitude 
structure  is  repeated  at  the  bottom  of  the  figure.  It  seems  quite  obvious  that  there  is  too 
much  2,3  in  the  model  results;  the  observations  show  a  2,2-type  behavior.  This  illustrates  how 
valuable  observations  are  to  these  theoretical  tidal  studies:  the  observations  very  much  guide 
our  modelling  effort.  After  attempting  to  better  match  these  solstice  observations,  equinox 
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observations  will  be  modelled,  followed  by  a  similar  study  for  geomagnet  icaily  active  conditions. 
Figure  13  summarizes  the  main  results  of  the  tidal  studies  presented  here,  which  are: 

-  tidal  effects  are  an  integral  part  of  the  thermospheric  circulation,  densit’es,  and  temper¬ 
atures;  this  requires  that 

-  theoretical  models  must  account  for  the  tides 

-  the  tides  vary  with 

-  season  (up  to  aM  over  25%  in  amplitude  for  densities,  depending  on  altitude) 

-  solar  cycle  (50-100%  and  more  in  density  amplitudes  at  high  altitudes) 

-  latitude  (see  below) 

-  the  tides  are  also  very  variable  from  day  to  day;  variations  may  be  on  the  order  of  30-40% 
in  the  density  amplitudes,  depending  on  altitude; 

-  tidal  influence  on  the  neutral  densities  are  shown  below,  for  three  different  heights. 
These  are  calculated  percentage  differences  in  the  results  with  tides  and  without  tides  at 
the  pole  and  at  the  equator: 

Z  %p  (pole  -  equator,  with  and  without  tides) 


120  km 

100-1000% 

180  km 

10-25% 

300  km 

10-40% 

-  to  improve  modelling  and  predictive  capabilities,  there  is  a  great  nc  -d  for  more  obser¬ 
vations,  meaning  more  frequent  observations,  and  from  more  sites;  these  should  extend 
from  the  lower  to  the  upper  thermosphere,  and  be  made  at  several  latitudes  in  order  to 
isolate  the  important  tidal  modes. 

-  the  major  obstacle  to  our  future  progress  in  these  tidal  studies  is  directly  attributable  to 
the  lack  of  data:  it’s  impossible  to  tell  how  good  (or  bad)  the  models  are  working  without 
substantial  amounts  of  data  for  the  conditions  being  simulated. 
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Tidal  Effects  on  Thermospheric 
Circulation,  Temperature,  and  Density 

Structure 


•  What  are  tides? 

•  Why  are  they  important? 

•  Where  do  they  come  from? 

•  How  do  we  put  them  in  models? 

•  What  happens  when  we  do  get  them  in  the 
models? 

•  How  are  we  doing? 

•  What  else  needs  to  be  done? 


Figure  1. 
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HEIGHT  (to**) 


Sources  of  Thermospheric  Tides 


1 .  In-situ  absorption  of  solar  UV  and  EUV  energy 

2.  Upward  propagating  tides  from  the  stratosphere 
and  mesosphere 

3.  Ion-neutral  momentum  coupling 


Figure  3. 
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APPROXIMATE  ALTTTUOT  (Vn) 


SEMIDIURNAL  HOUGH  FUNCTIONS 
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Incorporating  Tides  in  the  TGCM: 


1 .  Specify  geophysical  conditions 

2.  Run  the  model  to  equilibrium  for  solar  forcing 
only 

>3.  Put  in  the  semidiurnal  tides  at  the  lower 
boundary  (97km) 

4.  Run  the  model  to  equilibrium 

5.  Compare  T,  VN.S  to  observations: 
presently,  Arecibo  (18°N),  Saint  Santin  and 

Millstone  Hill  (both  ~45°N) 


Figure  6. 
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Figure  12. 


Summary 


•  Tidal  effects  are  an  integral  part  of  the 
thermospheric  circulation,  densities,  and 
temperatures 

t 

Theoretical  models  must  account  for  the  tides 

•  The  tides  vary  with 

•  Season  (up  to  and  over  25%  in  amplitude, 
depending  on  altitude) 

•  Solar  Cycle  (  >50%  at  higher  altitudes) 

•  Latitude,  Altitude 

•  The  tides  are  also  very  variable  from  day  to  day: 
Variations  may  be  on  the  order  of  30-40%  in  the 
density  amplitude 

•  Tidal  influence  on  the  neutral  densities 
(fixed  season,  solar  cycle) 

Z  %  p  (pole/equator) 

120km  100-1000% 

180  -1 0-25% 

300  1 0-40% 

•  To  improve  modelling  and  prediction,  need 
observations: 

at  several  levels,  at  several  latitudes 


Figure  13. 
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Thermosphere  Dynamics  and  iLs  influence _ on  atmospheric  density:  the 

Dynamics  Explorer  experience 
T.  L.  Kiileen, 

Space  Physi-s  Research  Laboratory. 

University  of  Michigan. 

Improvements  in  the  accuracy  with  which  atmospheric  densities  (100-250  km)  are  currently 
modeled  may  be  expected  to  occur  due  to  two  separate,  but  related,  efforts.  The  fust  of  these  is  in  the  am  ol 
empirical  or  semi-empiricai  models.  Such  models  are  based  on  fits  to  a  large  body  of  data  using  analytical 
expressions.  These  expressions  are  typically  modified  through  consideration  of  physical  principles 
associated  with  the  domain  under  consideration.  For  example,  the  assumption  of  diffusive  equilibrium  is 
often  used  for  empirical  modeling  of  thermospheric  composition.  This  assumption  provides  a  constraint  on 
the  analytical  expressions  chosen.  Empirical  models  can  be  improved  by  increased  resolution  (i.e.,  greater 
number  of  harmonics),  increased  experimental  coverage  of  the  fitted  parameters,  and/or  improvements  in  u.e 
analytical  expressions  used  to  convey  the  information.  To  improve  accuracy  of  empirical  models  in  the 
lower  thermosphere  in  the  next  several  years  will  require  full  use  of  the  existing  oata  set  as  well  as  efforts  to 
provide  additional  systematic  coverage  of  composition  and  densities  in  the  region,  i  measurements  are 
required  to  limit  the  use  of  extrapolations  from  higher  altitudes. 

The  second  effort  is  in  i!ie  area  of  first-principles  theoretical  models,  such  as  tht.mospheric  general 
circulauon  models  (TGCMs).  The  use  of  realistic,  time-dependent  boundary  conditions  and  inputs  of  energy 
and  momentum  within  the  framework  of  a  TGCM  can  lead  to  reasonably  accurate  simulations  of  the 
thermospheric  response,  involving  all  thermospheric  state  variables  (winds,  temperatures,  densities, 
composition,  etc.).  Improvements  i.ere  will  come  about  due  to  increased  sohpisticauon  of  the  models,  1 

improved  realism  of  the  input  parameters  used  and,  most  importantly,  improved  understanding  of  the 
physical,  chemical,  radiative  and  dynamical  processes  responsible  for  the  modeled  structure. 

This  talk  addresses  the  second  of  the  two  approaches  outlined  above  and  makes  the  fundamental 
point  the'  experience  has  shown  that  TGCMs  are  powerful  tools  '.hat  need  comprehensive  experimental 
constraints  in  ord  r  to  be  tested  against  reality  and  (incrementally)  made  more  realisuc.  An  analogy  is  di  iwn 
between  anticipated  future  improvements  in  our  understanding  of  lower  thermosphere  structi  using 
TGCMs  and  the  recent  rapid  progress  made  in  modeling  the  dynamics,  composiuon  and  thermal  structure  of 
the  upper  thermosphere.  For  tne  case  of  the  upper  thermosphere,  data  from  the  Dynamics  Explorer-2  (DPT) 
spacecraft  have  bc^n  used  to  constrain  and  test  the  NCAR  TGCM  and  the  experiment-model  comparisons 
have  led  directly  to  improved  model  accuracy  PE-2  provided  direct  measurements  of  all  tn  -miosphcnc  state 
variables  at  F-region  altitudes  at  a  formative  umi  m  the  development  of  the  NC.AR-T  GCM.  Comparisons 
and  critical  evaluation  of  discrepancies  between  (lies  r  and  experiment  led  to  a  quantitative  understanding  cf 
the  major  upper  thermospheric  forcing  processes  and  sig.  ificantly  more  realistic  predictive  capabilio. 
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For  future  progress  to  be  made  in  understanding  the  structure  and  variability  of  the  thermosphere  at 
lower  altitudes,  similar  comprehensive  coverage  of  ALL  lower  thermospheric  state  variables  is  needed.  Such 
coverage  will  be  used  to  test  and  improve  the  TGCM  in  this  relatively  poorly  understood  region.  The  new 
data  are  needed  essentially  to  improve  our  understanding  of  the  basic  physical  controlling  mechanisms, 
w  ithout  which  understanding  no  rapid  theoretical  progress  can  be  made. 
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Thermosphere  Dynamics  and  its  influence 
on  atmospheric  density:  the  Dynamics 
Explorer  experience 

Our  understanding  of  F-region  thermosphere 
dynamics,  composition,  density  perturbations,  thermal 
structure  and  thermospheric  forcing  mechanisms  was 
improved  considerably  due  to  the  comprehensive 
coverage  of  essential  variables  by  DE-2  and  DE-1. 

Knowledge  of  magnetospheric,  solar,  and  lower 
atmosphere  inputs  to  the  thermosphere  -  as  well  as 
all  thermospheric  state  variables  -  essential. 

Interative  comparison  with  general  circulation 
models,  followed  by  detailed  examination  of  the 
physics  within  the  models,  has  led  to  a  much  superior 
predictive  capability  at  F-region  altitudes. 

Need  to  do  the  same  for  the  lower  thermosphere  in 
order  to  understand  and  thereafter  model  the  density 
perturbations.  Need  comprehensive  experimental 
coverage  as  well  as  sophisticated  theoretical  tools. 
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Changes  in  the  thermospheric  mass  mixing  ratio  along  a 
constant  pressure  gradient  with  respect  to  time  are  given  by 

Dickinson  et  al.,  [1984]  as: 
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Tlie  terms  on  the  right  hand  side  of  this  equation  are  from  left  to 

right: 


molecular  diffusion  (basic  restoring  force): 
eddy  diffusion: 

horizontal  advection  (effect  of  horizontal  transport): 
vertical  advection  (upwelling): 
photodissociation: 
recombination: 
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3.5 


THERMOSPHERE  DYNAMICS:  RECENT  CONTRIBUTIONS  FROM  DYNAMICS  EXPLORER-2 


T.  L.  KILLEEN*  and  R.  G.  ROBLE+ 

Instrumentation  Down  oo  the  Dynamics  Explorer-2  (DE-2)  spacecraft  enabled  the 
dynamics  of  both  the  neutral  and  ionized  components  of  the  Earth's  upper 
atmosphere  to  be  monitored  over  (he  lifetime  of  the  spacecraft  from  August  1981 
to  February  1983.  The  direct  measurements  of  global  thermospheric  vector  neutral 
winds  and  iou  drifts  were  supplemented  by  observations  of  neutral  and  ionic 
constituent  abundances  and  temperatures,  precipitating  particle  fluxes  and  pitch¬ 
ing!*  di^ribuUO”,  electric  and  magnetic  fields,  and  global-scale  auroral 
luminosity  distributions.  Analysis  of  this  unique  and  comprehensive  global  data 
base  has  led  to  advances  in  our  understanding  of  the  manner  in  which  the  Earth's 
thermosphere  responds  dynamically  to  the  insertion  of  energy  and  momentum 
from  the  magnetosphere.  In  this  paper,  we  review  the  scientific  progress  achieved 
over  the  past  five  years  through  the  efforts  of  members  of  the  Dynamics  Explorer 
science  team  and  other  interested  scientists  who  have  used  the  DE  data  base  to 
investigate  the  dynamics  and  composition  of  the  thermosphere. 

INTRODUCTION 

The  Dynamics  Explorer  (DE)  dual-spacecraft  mission  represents  an  important  element  of  the  National 
Aeronautics  and  Space  Administration  (NASA)  program  dedicated  to  improving  our  understanding  of  the 
processes  that  control  man's  near  space  environment.  The  DE  mission  developed  out  of  scientific 
discussions  in  the  late  1960's  and  early  1970's  during  which  it  became  dear  that  a  full  understanding  of  the 
near-space  environment  would  require  comprehensive  experimental  and  theoretical  studies  of  the  physical, 
chemical,  dynamic  and  radiative  processes  that  couple  the  Earth's  magnetosphere,  ionosphere,  and  upper 
atmosphere  -  and  not  just  independent  studies  of  the  individual  regions.  A  mission  plan  was  subsequently 
designed,  involving  two  companion  spacecraft  in  co-planar  orbits:  a  high-altitude  spacecraft  (DE  1) 
instrumented  to  probe  the  Earth's  inner  magnetosphere  and  plasmasphere  and  a  low-altitude  spacecraft  (DE 
2)  instrumented  to  probe  the  ionosphere  and  neutral  upper  atmosphere  simultaneously  with  the  DE-1 
measurements.  The  principal  mission  objective  was  to  investigate  the  strong  interactive  processes  coupling 
the  hot,  tenuous,  coovecdng  plasmas  of  the  magnetosphere  with  the  cooler,  denser  plasmas  and  neutral 
gases  of  the  Earth's  upper  atmospheric,  ionosphere,  and  plasmasphere.  To  pursue  this  task,  secondary 
mission  objectives  were  defined  and  several  of  these  relate  specifically  to  the  investigation  of  the  dynamical 
response  of  the  upper  atmot  V  to  the  input  of  energy  and  momentum  from  the  magnetosphere.  A 
question  of  particular  interc.  <  ur  r-ves  the  processes  of  momentum  and  energy  transfer  between  the 
con  veering  ionospheric  kns  (dr*  y  high-latitude  electric  fields  of  magnetosphere:  origin)  and  the  neutral 

gas  particles  in  the  thermosphere.  An  overview  of  the  main  scientific  goals  for  the  DE  program  has  been 
given  by  Hoffman  and  Schmeriing  (1981)*. 

In  this  paper,  we  present  a  synopsis  of  the  progress  made  during  the  first  five  yean  of  the  DE  mission,  with 
emphasis  placed  on  the  new  insights  obtained  through  use  of  the  unique  data  base  of  thermospheric  neutral 
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wind  observations.  Review  articles  touching  on  early  DE  results  have  been  presented  previously  2*3*4.  In 
the  next  section,  we  briefly  describe  the  DE-2  instrumentation  that  has  permitted  the  vector  neutral-wind 
observations  and,  in  the  subsequent  sections,  we  describe  aspects  of  the  scientific  work  accomplished  using 
the  DE  data  base,  including  the  very  important  parallel  theoretical  model  development  and  the  ongoing  data- 
theory  comparative  studies  that  have  been  a  hallmark  of  the  program.  Finally,  we  summarize  the  state  of 
the  field  after  the  first  five  years  of  analysis  of  data  from  the  Dynamics  Explorer  mission 

INSTRUMENTATION 

Two  instruments  on  DE  2,  each  of  which  measures  a  single  orthogonal  component  of  the  horizontal  vector 
neutral  wind,  when  combined,  provide  the  horizontal  neutral-wind  vector.  These  instruments  are  the  Fabry- 
Pcrot  interferometer,  FPI-\  and  the  Wind  and  Temperature  Spectrometer,  WATS^.  In  this  section  we 
briefly  describe  the  measurement  technique  for  FPI  and  WATS;  a  more  detailed  discussion  of  these  and  the 
other  DE  instruments  can  be  found  in  the  articles  published  in  the  special  issue  of  Space  Science 
Instrumentation  (VoL  5, 1981)  devoted  to  the  DE  miss  km  objectives  and  instrumental  details. 

ThcFabry-PerQilmcffsrpmcigf  qhDE.2. 

The  Fabry-Perot  interferometer  on  DF  2  measures  the  meridional  (in-orbit-plane)  component  of  the  neutral- 
wind  vector  and  the  thermospheric  kinetic  temperature  as  a  function  of  altitude  along  and  below  the  satellite 
track.  A  full  description  of  the  instrument  has  been  given  previously^.  Fig.  1  illustrates  schematically  the 
measurement  viewing  geometry  for  the  optical  remote-sensing  device.  FPI  measures  altitude  profiles  of  the 
neutral  wind  and  the  nnitral  kinetic  temperature  ahead  of  and  below  the  spacecraft  by  remotely  sensing  the 
Doppler  shift  and  Doppler  line-broadening,  respectively,  for  the  0(*D)  thermospheric  emission  line  profile 
at  630  nm.  The  instrument  field  of  view  is  directed  along  and  below  the  spacecraft  velocity  vector  and 
horizon  scans  of  the  Earth's  limb  are  performed  in  the  orbital  plane  using  an  internal  scanning  mirror.  The 
most  critical  component  of  the  FPI  is  a  plane  Fabry-Perot  etaloo  which  provides  the  required  ultra-high 
spectral  resolution  necessary  to  quantify  the  Doppler  characteristics  of  the  measured  emission  line  profiles. 
The  measured  Doppler  shift  is  converted  into  a  line-of-sight  neutral  wind  measurement  by  subtracting  that 
component  of  the  shift  due  to  the  known  spacecraft  velocity.  The  ambient  kinetic  temperature  of  the 
atmosphere  is  derived  by  deconvolving  the  calibrated  instrumental  broadening  function  from  the  measured 
signal,  leading  to  the  characterization  of  the  Gaussian  half-width  of  the  emission  line^. 

IbsLWiiKl  and  Tcmpgature  Spcctromflg 

The  Wind  and  Temperature  Spectrometer,  WATS,  has  been  described  in  detail  previously^.  The  instrument 
is  a  modified  closed-source  mass  spectrometer,  with  scanning  baffles  providing  modulation  of  the  detected 
signal  from  which  winds  and  temperatures  can  be  derived.  It  is  located  on  the  spacecraft  base-plate  such  that 
the  instrument  aperture  faces  towards  the  ram  direction  and  neutral  gas  particles  are  swept  into  the  mass 
spectrometer  chamber  by  the  motion  of  the  spacecraft.  The  basic  measurement  technique  is  illustrated 
schematically  in  Fig.  2.  Each  scanning  baffle  is  commanded  to  move  across  the  instrument  aperture  at  a 
precisely-determined  speed  during  measurements  of  the  abundances  of  the  major  thermospheric  constituents, 
such  as  O  and  N2-  The  exact  baffle  angle  at  which  the  incoming  (ram)  beam  of  neutral  particles  is 
intercepted  maximally  is  a  function  of  the  cross-track  wind  component  in  the  plane  of  the  baffle  motion 
(Le.,  cither  in  the  vertical  or  in  the  zonal  direction,  depending  on  which  baffle  is  exercised).  The  wind 
component  can  be  simply  obtained  by  solving  the  vector  triangle  shown  in  Fig.  2,  where  the  solution 
obviously  requires  very  accurate  knowledge  of  both  the  speed  and  direction  of  motion  of  the  spacecraft.  The 
spacecraft  velocity  vector  is  considered  to  be  known  sufficiently  well  in  the  general  case  to  provide  WATS 
wind  accuracies  within  ~±15m/sec,  Measurements  with  this  precision  can  be  made  whenever  the  spacecraft 
is  sufficiently  low  in  altitude  to  provide  statistically  significant  neutral  particle  counting  rates.  Typically, 
this  factor  constrains  the  WATS  wind  measurements  to  altitudes  below  about  600  km. 

The  neutral  kinetic  temperature  is  also  obtained  from  an  analysis  of  the  modulation  of  the  WATS  count  rate 
js  a  function  of  baffle  location.  The  shape  of  the  drop-off  in  count  rate  during  the  interception  of  the 
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incoming  beam  of  neutral  particles  by  ibe  scanning  baffle  is  a  function  of  the  transverse  velocity 
distribution  of  the  particles  and,  therefore,  of  the  ambient  neutral  kinetic  temperature.  The  temperature  can 
oe  recovered  from  the  measured  count  rate  profile  using  a  least  squares  fitting  technique^. 

Neutral  wind  vector  measurements 

The  EFT  and  WATS  instruments  briefly  described  above  are  complementary  in  that  they  provide  independent 
measurements  of  the  two  components  of  the  horizontal  neutral  wind.  A  technique  has  been  developed®  that 
merges  the  data  from  both  instruments  to  provide  a  description  of  the  horizontal  vector  of  the  neutral  wind 
along  the  orbital  track  of  DE  2.  As  an  example  of  the  utility  of  the  horizontal  vector  wind  measurement 
from  DE  2,  Fig.  3  shows  wind  vectors  from  FPI  and  WATS  obtained  during  a  single  orbit  of  DE  2  where 
the  vector  winds  arc  plotted  for  (a)  the  Southern  Hemisphere  polar  cap  and  foi  (b)  the  Northern  Hemisphere 
polar  cap.  The  data  in  both  hemispheres  provide  ar.  indication  of  the  characteristics  of  the  circulation 
pattern,  showing  high-speed  anti-sunward  flow  (-500  m/sec)  in  the  central  polar  cap  bounded  by  sunward 
flow  in  both  the  dawn  and  dusk  sectors.  This  pattern  of  thermospheric  flow  is  similar  to  the  classical 
ionospheric  convection  pattern  driven  by  magnetospheric  electric  Gelds,  with  two  contra-rotating  vortices  at 
high-latitudes.  Vector  wind  measurements  from  FPI  and  WATS  such  as  these  have  been  used  in  many 
experimental  and  comparative  theoretical  studies,  reviewed  in  more  detail  below. 

VECTOR  WIND  MEASUREMENTS  AND  MEAN  HIGH-LATITUDE  THERMOSPHERIC 
CIRCULATION 

The  availability  of  global-scale  thermospheric  neutral  wind  vector  measurements  from  DE,  such  as  those 
described  above,  has  led  to  a  number  of  studies  over  the  past  few  years  that  have  concentrated  on  providing 
the  first  large-seal*,  morphological  description  of  the  high- latitude  thermospheric  circulation  system.  We 
next  review  the  contributions  made  in  this  area  through  use  of  DE  data.  We  note  at  the  outset  that  the  wind 
'mapping*  studies  are  for  solar  maximum  conditions  existing  during  the  DE-2  lifetime  when  the  effects  of 
the  energy  and  momentum  coupling  between  the  magnetosphere,  ionosphere  and  thermosphere  might  be 
expected  to  be  at  a  maximum.  It  would  be  incorrect  and  misleading  to  assume  that  the  DE  measurements 
provide  a  correct  description  of  thermospheric  circulation  for  other  parts  of  the  solar  cycle. 

The  Universal  Tunc fffcfi at  high-latitodcs 

One  of  the  first  critical  tests  performed  using  the  FPI/WATS  vector  wind  measurement  capability  was  to 
check  for  'ground-truth*  by  comparing  DE-2  wind  observations  with  simultaneously-determined  and 
common-volume  ground-based  measurements  of  winds.  This  study  was  carried  out?  using  the  coordinated 
measurements  of  R.  W.  Smith  and  colleagues  from  the  Ulster  College  ground-based  Fabry-Perot 
interferometer  (UFPI)  located  in  Svalbard,  Norway  for  overflights  and  near  overflights  of  the  DE-2 
spacecraft  The  UFPI  and  DE-2  observations  were  found  to  be  in  satisfactory  agreement  for  the  set  of 
simultaneous  near-common- volume  measurements  available,  indicating  that  the  ground-based  and  satellite 
techniques  were  both  providing  realistic  and  consistent  determinations  of  the  thermospheric  flow. 

Having  established  the  compatibility  of  the  ground-based  and  spacecraft  observation  techniques,  these 
authors  went  on  to  generate  snapshots  of  the  thermospheric  circulation  at  high-latitudes  by  using 
simultaneous  data  from  the  UFPI  and  DE-2  for  periods  when  the  instruments  were  separated  spatially.  These 
snapshots  illustrated  the  important  role  played  by  the  ion-drag  momentum  force  in  establishing  the  twin- 
cell  neutral  thermospheric  circulation  pattern  by  demonstrating  basic  similarities  between  the  measured 
high-latitude  thermospheric  neutral  wind  system  and  the  classical  twin-cell  ionospheric  convection  pattern. 
Since  the  latter  is  controlled  by  magnetospheric  electric  fields  mapping  into  the  ionosphere  along 
geomagnetic  field  lines,  it  was  natural  to  expect  that  the  component  of  the  neutral  wind  system  driven  by 
ion  drag  or  Joule  heating  at  high-latitudes  might  also  display  some  degree  of  geomagnetic  control  To  study 
this  question,  additional  studies  were  conducted  using  averaged  neutral  wind  data  from  many  DE-2 
orbital  passes.  Firstly,  a  statistical  analysis  of  neutral  wind  data  from  58  orbital  passes  of  DE  2  over  the 
Southern  Hemisphere  polar  region  was  performed The  data,  which  all  came  from  a  single  LT  plane,  were 
averaged  into  bins  of  5  degrees  in  laiinvV-  and  two  hours  in  UT  to  illustrate  the  basic  diurnal  coherence  in 
the  neutral  wind  pattern  as  seen  from  a  polar  orbit.  It  was  also  shown  how  the  high-latitude  wind  field  could 
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be  described  experimentally  in  terms  of  a  superposition  of  a  "background*  wind  field  driven  by  aolar  heating 
with  a  wind  field  driven  by  high-latitude  momentum  **»d  energy  sources. 

la  a  second  and  more  comprehensive  study,  this  approach  was  extended^  to  investigate  tbe  mean  UT- 
depcodeot  neutral  winds  in  both  summer  and  winter  hemispheres,  showing  dearly  the  basic  geomagnetic 
control  of  the  higb-Iatimde  thermospheric  wind  pattern.  The  mean  thermospheric  winds  were  also  compared 
with  tbe  predictions  of  the  two  thermospheric  general  drcnlatioe  models  (TGCMs),  namely  tbe  University 
College  London  (UCL-TGCM)12  and  the  National  Center  for  Atmospheric  Research  (NCAR- 
TGCM)13.14  The  similarities  and  differences  between  the  experimental  observations  and  tbe  theoretical 
model  predictions  indicated  that  the  models  provided  reasonable  predictions  for  the  mean  UT-d'pendert  upper 
thermospheric  winds  and  that  they  contained  a  reasonably  valid  description  of  the  dhirnally-modulated  ion 
drag  momentum  force.  This  type  of  global-scale  experiment-theory  comparison  has  been  very  important  for 
essentially  validating  tbe  physical  mechanisms  incorporated  within  the  TGCM  model  descriptions,  for 
identifying  physical  processes  missing  from  the  theoretical  formulations,  and  for  leading  to  improvements 
in  the  input  prescriptions  and  parameterizations  used  in  tbe  models. 


H  »v«<  ill  npi  ia  STiTuTi  [, 


In  a  later  study^,  all  available  F-regioo  neutral  wind  data  from  seven  ground-based  Fabry-Perot 
interferometers  (GBFPls)  as  well  as  from  several  hundred  DB-2  orbital  passes  taken  during  December  1981 
were  collated  in  order  to  provide  tbe  fullest  possible  experimental  description  of  tbe  mean  thermospheric 
circulation.  Since  weather  conditions  at  the  ground  stations  limited  the  number  of  occasions  when 
simultaneous  observations  were  possible,  the  data  were  averaged  by  UT  to  provide  12  thermospheric 
“climatological  maps'  for  every  two  boms  of  UT  during  a  "mean  day"  in  December  1981.  Fig.  4  presents 
three  of  the  experimentally-derived  maps  Cor  the  UT  bins  centered  on  07:00.  09jOO  and  11  fX)  hrs  UT  as  well 
as  the  theoretical  predictions  from  the  NCAR-TGCM  in  a  common  format.  Here;  the  averaged  PE  winds  are 
shown  as  arrows  and  the  ground-based  winds  are  given  by  the  meteorological  symbols.  It  was  concluded 
firstly,  that  DE  measurements  and  data  from  multiple  GBFPI  stations  could  be  used  to  provide  consistent 
determinations  of  the  global-scale  thermospheric  circulation  pattern  and  secondly,  that  the  mean  diurnal 
wind  pattern  in  the  northern  hemisphere  high-latitude  region  at  winter  solstice  could  be  well  modeled using 
appropriate  model  inputs  for  the  mean  thermospheric  forcings.  This  latter  conclusion  was  important  since  it 
implied  that  the  TGCMs  could  be  used  with  confidence  to  study  quantitatively  the  forcing  processes 
responsible  for  the  mem  circulation  and,  furthermore,  could  be  used  to  predict  the  mean  circolatioo  in  areas 
where  measurements  are  not  available.  Tbe  successful  theoretical  prediction  of  the  mean  circulation 
characteristics  is  also  useful  in  that  it  provides  a  "baseline"  or  "reference"  pattern  for  use  in  studies  of  storm¬ 
time  perturbations  where  large  departures  from  the  mean  can  oocur. 

It  was  noted  in  this  and  other  studies  that  the  anti-cyclonic  vortex  in  the  thermospheric  circulation  pattern, 
associated  with  the  das k  cell  of  ionospheric  convection,  is  almost  invariably  stronger  than  the  cyclonic 
vortex  associated  with  the  dawn  ion  convection  cefL  This  large  asymmetry  was  ascribed^  to  the  effects  of 
the  Coriolis  force  which  tends  to  destroy  cyclonic  vorticity  and  maintain  anticyclonic  vorticity  in  both 
hemispheres,  in  accord  with  an  earlier  theoretical  study  *6.  The  dominance  of  the  dusk  cell  in  the  neutral 
circulation  pattern  is  also  due,  at  least  in  part,  to  tbe  tendency  of  the  dusk  cell  of  iooospheric  convection  to 
be  associated  with  a  greater  potential  drop  than  the  dawn  ceQ^. 

In  an  analogous  investigation  to  tbe  wind  studies  mentioned  above,  neutral  composition  data  from  the 
Neutral  Atmosphere  Composition  Spectrometer  (NACS)  on  DE-2  was  used  18  to  investigate  the 
morphology  of  thermospheric  composition  variations  for  N2  and  He  in  tbe  polar  regions  during 
geomagnetically  quiet  conditions.  Averaged  measurements  from  the  NACS  instrument  were  compared  with 
the  MSIS-83  empirical  model^  of  neutral  composition.  It  was  found  that  tbe  density  for  the  heavier 
species  (N2)  was  enhanced  in  the  general  vicinity  of  tbe  geomagnetic  poles,  with  tbe  enhancement  more 
closely  confined  to  the  polar  regions  in  winter.  The  lighter  species  (He)  was  depleted  in  the  same  regions, 
with  maximum  density  response  in  the  morning  hours  (magnetic)  on  avenge.  The  distribution  of  densities 
about  the  magnetic  poles  was  found  to  be  dependent  on  UT.  These  results,  when  combined  with  the 
morphological  wind  and  temperature  studies  using  DE  data,  show  bow  tbe  mean  thermospheric 


composition,  temperature  structure,  sod  dynamics  all  have  clear  signatures  that  may  be  related  to 
magnaospheric  forcing. 

Interplanetary  Mayn^ie  Reid  dependency  of  high-latitude  thermospheric  winds 

To  follow  oo  from  the  morphological  and  mapping  studies  mentioned  above,  attempts  have  been  made  to 
identify  additional  physical  processes  and  mechanisms  that  can  provide  order  to  the  large  DE  data  base  of 
neutral  wind  observations.  Since  the  neutral  winds  at  high  latitudes  were  found  to  partially  mimic  the  twin¬ 
cell  ioo -convection  pattern,  and  since  previous  work2*'  had  indicated  that  the  geometry  of  the  ion- 
coovcctioo  pattern  is  a  strong  function  of  the  direction  and  magnitude  of  the  Interplanetary  Magnetic  Field, 
IMF,  several  studies  were  carried  out  to  test  the  sensitivity  of  the  measured  neutral  wind  circulation  to  the 
magnitude  and  orientation  on  the  IMF. 

Firstly,  the  dependency  of  the  neutral  circulation  on  the  east-west  (By)  component  of  the  IMF  was 

investigated24'2^.  Both  single-orbit  (lata  as  well  as  averaged  data  from  many  DE-2  orbital  passes  over  the 
Northern  Hemisphere  polar  region  were  used  to  illustrate  changes  in  the  configuration  of  the  neutral 
circulation  pattern  in  geomagnetic  coordinates  that  were  clearly  ordered  by  the  IMF  By  component.  To 
accomplish  this  work  measurements  of  the  IMF  from  the  IS  EE- 3  spacecraft  were  used  to  select  from  and 
order  the  DE  neutral  wind  vector  data.  From  a  study  of  individual  DB  orbital  passes,  a  "W" -shaped  envelope 
to  the  anti-sunward  polar  cap  neutral  flow  was  reported24  that  resembled  similar  By-dependent  signatures  in 
ionospheric  convection.  In  a  lata  study2**,  *  significantly  larger  data  set  was  used  to  characterize  the  By- 
depcadencc  of  the  neutral  circalatioo  for  both  southern  (summer)  and  northern  (winter)  hemispheres.  Results 
from  this  study  are  shown  in  Fig.  5,  where  data  from  many  hundreds  of  DE  2  orbital  passes  over  six 
months  of  the  mission  have  been  separated  according  to  the  sign  of  By  and  averaged  into  bins  of 

g<Tuna£Tv-fir  latitiuV.  wnA  gmm»gTWk-.  local  tiny,  for  the,  two  tvmkphffre^  These  datm  were  nhtliiv*!  daring 

the  nwmth«  November  1981  to  January  1982  and  November  1982  to  January  1983.  The  criteria  used  to 
separate  the  data  sets  by  By  involved  the  pre-existence  of  a  definite  positive  or  negative  value  for  By  as 
measured  by  ISEE-3  for  1  hour  prior  to  the  DE  orbital  pass.  Passes  occurring  either  during  very  high  or 
very  low  levels  of  geomagnetic  activity  were  excluded  from  the  study  and  allowance  was  made  for  the 
propagation  time  to  the  Earth  of  the  IMF  measured  at  ISEE-3  altitudes.  As  can  be  seen  from  fig.  5,  there 
are  significant  differences  in  the  mean  thermospheric  circulation  patterns  for  By  positive  and  negative.  To 
some  extent  the  effects  are  mirrored  between  hemispheres;  that  is.  a  By  positive  (negative)  signature  in  the 
northern  hemisphere  resembles  a  By  negative  (positive)  signature  in  the  southern  hemisphere.  In  all  cases, 
the  average  circulation  patterns  shown  illustrate  the  dominance  of  the  dude  anticyclonic  neutral  vortex 
mejituwri  above.  The.rparial  magnitude  of  this  vortex,  however,  is  greater  far  By  positive  conditions  in  the 
northern  hemisphere  and  for  By  negative  conditions  in  the  southern  hemisphere.  A  distinct  rotation  in  the 
direction  of  the  polar  op  anti-sunward  flow  that  was  dependent  on  the  sign  of  By  was  also  noted. 

A  clear  dependence  of  the  neutral  wind  field  as  measured  from  DB  2  oo  the  north-south  (B^)  component  of 
the  IMF  was  also  demonstrated  in  a  separate  study27.  For  conditions  of  strong  and  steady  northward 
(positive)  Bz.  neutral  and  ion  winds  were  presented  which  showed  that  sunward  flow  of  neutral  gas  in  the 
center  of  the  polar  cap  cm  be  established  and  maintained  by  ion  drag  for  positive  Bz. 

DE  neutral  wind  measurements  have,  therefore,  demonstrated  clear  morphological  dependencies  on  the 
orientation  of  the  IMF  and  these  dependencies  are  similar  in  nature  to  those  previously  observed  for 
ionospheric  convection.  The  identification  and  quantification  of  both  large  and  more  subtle  dependencies  of 
the  neutral  wind  circulation  pattern  oo  the  orientation  of  the  IMF  has  been  one  of  the  most  remarkable 
results  from  the  DE  mission  in  that  it  clearly  illustrates  the  strong  coupling  that  occurs  throughout  the 
Earth's  near-space  environment  all  the  way  from  the  solar  wind,  through  the  magnetosphere  and  ionosphere 
and  into  the  neutral  ionosphere  itself.  Such  dependencies  would  perhaps  not  have  appeared  intuitively 
obvious  prior  to  the  DE  mission.  Ongoing  but  as  yet  unpublished  studies  have  shown  that  the 
thermospheric  composition  and  density  structure  at  high  latitudes  is  also  strongly  dependent  on  the 
orientation  of  the  IMF. 

Geomagnetic  dependency  of  mean  circulation 
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The  burning  and  averaging  studies  reviewed  above  tend  to  mask  the  high  degree  of  variability  that  occurs  in 
high- latitude  thermosphere  dynamics  as  a  function  of  geomagnetic  activity.  This  variability  was  first 
highlighted  using  DE  observations  by  Spencer  a  aL  (1982)^  who  presented  WATS  zonal  and  vertical  wind 
measurements  in  the  polar  region  that  demonstrated  the  existence  of  strong  vertical  winds,  having  velocities 
of  100- 150  m/sec  at  times,  correlated  with  strong  (many  hundreds  of  m/sec)  zonal  winds.  They  also  reported 
high-Laritude  temperature  perturba/'oos  of  several  hund»*«i  degrees  and  associated  these  with  gravity  waves. 

In  order  to  quantify  the  geomagnetic  dependency  of  high-latitude  neutral  dynamics,  a  statistical  study  using 
6  months  of  DE-2  data  was  recently  conducted^  specifically  to  determine  the  relationship  between  the 
mean  circulation  pattern  in  the  summer  and  winter  high-latitude  regions  and  the  level  of  geomagnetic 
activity,  as  represented  by  the  Kp  and  AE  indices.  Data  were  sorted  and  binned  according  to  these 
geomagnetic  indices  and  the  results  from  this  study,  shown  in  Figs.  6  and  7,  illustrate  bow  the  anti- 
sunward  neutral  winds  in  the  central  polar  cap  increase  in  speed  as  a  function  of  enhanced  geomagnetic 
activity  levels.  The  values  shown  in  Tables  1  and  II  are  least  squares  fits  to  the  specified  regression 
relationships  between  the  various  avenged  high-latitude  wind  features  and  the  geomagnetic  indices.  This 
study  may  be  considered  to  be  a  first  step  towards  a  satisfactory  empirical  description  of  the  high-latitude 
neutral  wind  field,  enabling  experimentally-based  wind  values  to  be  used  for  different  geophysical  conditions 
in  relevant  studies  and  providing  a  challenge  to  theoretical  modelers  to  match  quantitatively  the  observed 
dependencies. 

ION-NEUTRAL  COUPLING  AT  HIGH-LATITUDES 

The  various  DE  studies  of  high-latitude  circulation  reviewed  above  have  indicated  the  important  role  played 
by  the  high  latitude  sources  of  momentum  (ion  drag)  and  energy  (Joule  beating)  in  the  establishment  of  the 
mean  thermospheric  wind  system.  A  primary  conclusion  is  that  above  60°  of  iaritude  the  neutral  winds  in 
the  upper  thermosphere  are  best  described  in  a  geomagnetic  coordinate  system.  That  is,  the  thermosphere 
wind  system  responds  primarily  to  magnetospheric  driving  functions  that  map  into  the  auroral  region 
surrounding  the  magnetic  poles  of  the  Earth.  In  order  to  study  these  processes  further  a  large  set  of 
atmospheric,  ionospheric  and  magnetospheric  parameters  must  be  considered.  For  example,  the  ion-drag  and 
Joule  hearing  driving  functions  are  dependent  on  the  magnitude  of  the  ion-neutral  collision  frequency,  as 
well  as  on  the  differential  velocity  of  the  ions  and  neutrals.  Thus,  in  addition  to  the  Deutral  winds 
themselves,  the  effects  of  particle  precipitation,  convection  electric  fields,  ionospheric  recombination,  and 
neutral  composition  and  temperature  changes  must  all  be  taken  into  account  to  provide  the  framework  for 
auartitative  studies  of  the  dynamical  influence  of  magnetospheric  processes  on  the  upper  atmosphere.  The 
comprehensive  set  of  instrumentation  on  DE  2  has  enabled  valuable  initial  studies  of  this  type  to  be  carried 
out  Wc  next  review  these  studies. 

Data  from  six  instruments  on  DE  2  were  collected  to  investigate  the  detailed  ionospheric  and  neutral 
atmospheric  structure  far  a  set  of  orbital  passes  over  the  southern  hemisphere  polar  region™  These  data 
enabled  the  in-si  tu  ico  acatral  vector  difference  velocity  and  the  ion-neutral  collision  frequency  to  be 
obtained  simultaneously  for  the  first  tune,  leading  to  a  determination  of  the  insonrimcocc  local  rates  of 
transfer  of  momentum  and  energy  between  the  species.  The  simultaneously-measured  km  drifts  and  neutral 
winds  indicated  that  the  spatial  morphology  and  temporal  behavior  of  the  momentum  and  energy  coupling 
between  species  is  complex,  with  the  neutrals  responding  to  magnetosphericaUy-controlled  changes  in 
ionospheric  convection  on  time  scales  of  many  (lCTs  to  100’s)  minutes,  thereby  filtering  out  much  of  the 
high-frequency  structure  in  the  ion  convection.  Fig.  8  illustrates  this  last  point  well,  showing  two 
examples  of  simultaneous  zonal  neutral  winds  and  ion  drifts  measured  along  the  DE-2  orbital  track  as  the 
spacecraft  crossed  the  northern  hemisphere  polar  region.  In  this  figure,  the  ion  and  neutral  winds  have  been 
plotted  with  the  same  spatial  and  temporal  resolution.  It  can  be  seen  that  the  neutral  wind  signatures 
resemble  tliose  of  the  ion  drifts,  with,  however,  a  reduction  in  magnitude  (note  the  factor  of  3  difference  in 
scales  for  the  different  species)  and  with  much  less  evident  spatial  structure.  Thus,  because  of  the  long  (-1 
hour)  ioc-neutral  momentum  coupling  time  constant,  the  neutral  thermosphere  can  be  seen  to  act  as  an 
integrator,  partially  mimicking  the  iooospheric  convection  pattern,  but  displaying  reduced  amplitudes  and 
generally  less  spatial  structure. 

This  study  of  ion-neutral  momentum  transfer  time  constants  was  later  extended-*  1,  providing  statistical 
averages  for  sola,  maximum  (using  DE-2  measurements)  and  for  solar  minimum  (using  Atmosphere 


Explorer -C  measurements).  These  averaged  time  constants  were  compared  with  cakulatxxu  using  the 
MSIS-83  ^  and  the  Chiri32  semi-empirical  models  for  thermospheric  and  ionospheric  densities, 
respectively.  The  comparisons  showed  that  the  models  provided  values  for  the  time  constants  ax  solar 
maximum  that  were  in  good  agreement  with  the  measurements,  but  that,  at  solar  minimum,  the  modeled 
lime  constants  were  too  small,  indicating  that  the  ion-neutral  momentum  coupling  is  significantly  less 
important  than  that  predicted  for  low  solar  activity.  The  discrepancy  was  attributed  to  the  overestimation  of 
polar  ion  densities  at  solar  minimum  by  the  Chui  model.  This  study  pointed  lo  the  large  dynamic  range  of 
the  coupling  between  the  magnetosphere  and  thermosphere  which  occurs  as  a  function  of  time  within  the 
solar  cycle,  with  the  DE-2  mission  providing  information  during  a  period  when  thermospheric  effects  of  the 
coupling  were  at  a  maximum. 

Simultaneous  measurements  have  been  presented34  of  global-scale  auroral  luminosity  distributions  and 
vector  neutral  winds  over  the  northern  (winter)  polar  cap  using  data  from  the  spin-scan  auroral  imager, 
SA133  on  DE  1  and  the  FPI  and  WATS  instruments  on  DE  2,  respectively  (see  Fig.  9).  The  purpose  of  this 
work  was  to  illustrate  the  spatial  relationship  between  large-scale  morphological  features  of  the  F-region 
neutral  wind  field  in  the  winter  polar  region  and  the  location  and  spatial  extent  of  the  aurora.  A  definite 
correlation  was  found  to  exist  between  reversals  and  boundaries  in  the  neutral  wind  field  and  the  location  of 
the  visible  auroral  ovaL  The  dusk  sector  reversal  to  sunward  flow  in  the  neutral  gas  was  again  found  to  be 
more  pronounced  than  that  occurring  in  the  dawn  sector.  Ion  drift  and  composition  data  were  used  to 
determine  the  ion  drag  and  Coriolis  forces  on  the  neutral  gas  and  it  was  concluded  that  (he  dawn-dusk 
asymmetry  in  the  thermospheric  neutral  wind  response  was  related  to  the  differing  effects  of  the  Coriolis 
force  in  the  dawn  and  dusk  sectors.  Examples  of  such  simultaneous  data  sets  indicated  that  the  neutral  wind 
boundaries  follow  the  subs torm -dependent  expansion  and  contraction  of  the  auroral  oval.  The  superposition 
of  vector  neutral  wind  measurements  from  DE  2  and  auroral  imagery  from  DE  1  demonstrates  clearly  the 
power  of  the  DE  mission  to  determine  parameters  associated  with  the  input  of  energy,  momentum  and 
charged  particles  to  the  upper  atmosphere  as  well  as  its  dynamical  response  to  such  inputs. 

LOW-LATITUDE  THERMOSPHERIC  WINDS 

Measurements  of  thermospheric  winds  at  middle  and  low  latitudes  from  DE  2  have  been  studied  in  an  effort 
to  understand  better  the  mean  circulation  as  well  as  the  propagation  characteristics  of  wave?  and 
perturbations  originating  at  high  latitudes.  A  major  contribution  of  the  DE  2  data  has  been  in  the  area  of 
determining  experimentally  the  degree  of  atmospheric  superrotation  using  a  direct  technique.  Previous  work 
in  this  area  had  suggested  rather  large  values  for  the  superrotation  of  the  equatorial  thermosphere  based  on 
indirect  inferences  from  satellite  drag  data  (c.g.  Ref.  35).  Such  large  superrotation  velocities  at  low  latitudes 
presented  theoretical  difficulties^.  Using  WATS  zonal  wind  observations,  average  profiles  of  the  equatorial 
zonal  wind  for  different  levels  of  solar  dux  and  for  different  altitude  levels  were  obtained3^  and  the  results 
demonstrated  an  atmospheric  superrotation  of  only  a  few  percent  in  good  agreement  with  model 
calculations3®"4*).  Fig.  10  shows  the  local  time  dependency  of  the  WATS  zonal  wind  measurements  at  the 
equator  in  the  altitude  range  300-500  km.  As  can  be  seen,  while  there  is  a  significant  LT  variation  driven 
primarily  by  the  diurnal  variation  of  soiat  EUV  and  UV  beat  sources  as  well  er  momentum  sources 
associated  with  the  diurnal  tide3^  ,  the  average  value  is  very  near  to  zero.  The  table,  shown  in  Fig.  10. 
presents  the  results  for  die  mean  (LT-averaged)  zonal  wind  as  a  function  of  latitude  for  two  altitude  bins.  In 
both  cases,  the  latitudinally-averaged  (equator  -  40®  latitude)  zonal  wind  has  a  magnitude  of  less  than  2 
m/sec,  implying  that  the  equatorial  thermosphere,  on  the  average,  simply  uwotaies  with  the  Earth. 

THEORETICAL  MODELING 

At  the  time  of  the  Dynamics  Explorer  launch  in  August  1981.  two  thermospheric  general  circulation 
models  (the  UCL-TGCM12  and  NCAR-TGCM13)  had  been  developed  sufficiently  for  use  in  analyzing  data 
and  for  the  interpretation  of  physical  processes  responsible  for  the  observed  thermospheric  structure  and 
dynamics.  Both  models  solve  the  primitive  equations  of  dynamic  metcrology  modified  for  thermospheric 
dynamics  and  include  coupling  between  composition  and  dynamics. 

The  initial  DE  wind  vector  measurements  clearly  showed  a  high-latitude  wind  system  that  was  strongly 
controlled  by  the  momentum  ucurce  associated  with  magnetosphere  convection.  There  were  strong  anti- 
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sun*  •»-  Mods  in  the  magnetic  polar  cap  and  strong  sunward  winds  eqa&torwsrd  of  the  polar  cap  in  the 
amoral  zone.  Bod)  TGCMs  were  used  to  investigate  the  observed  wind  pattern  and  the  effectiveness  of  the 
magnctospheric  cocvccuoo -driven  ion-drag  momentum  source  in  controlling  the  wind  system  at  high 
latitudes.  The  UCL-TGCM  results  were  compared  with  observations  for  seven  polar  passes  of  the  DC-2 
spacecraft  during  November  198l4^.  It  was  found  thanbe  major  features  of  the  experimental  winds,  such  as 
the  mean  day-to-night  circulation  caused  by  solar  EUV  and  UV  heating,  augmented  by  magnctospheric 
processes  at  high  latitude,  and  the  sharp  boundaries  and  flow  reversals  imposed  on  thermospheric  winds  by 
ion-drag  momentum  transfer  from  the  magnetosphere,  were  qualitatively  predicted  by  a  version  of  the  model 
that  used  a  semi-empirical  model  of  polar  electric  fields.  A  comparison  of  winds  predicted  by  the  UCL- 
TGCM  using  two  different  ionospheric  models  of  electron  density  with  measurements  made  by  the  DE-2 
satellite  is  shown  in  Fig.  11.  The  dawn-to-dusk  cuts  show  antisun  ward  winds  in  the  polar  cap  flanked  by 
sunward  winds  in  the  vicinity  of  the  auroral  oval.  There  is  good  general  agreement  between  the  model 
predictions  and  observations  for  these  polar  passes. 

DE-2  d  tia  from  the  Oaober-November  1981  time  period  was  compared  with  predictions  of  the  NCAR- 
TGCM  •;’4  \  and  the  results  of  the  comparison  indicated  an  effective  control  of  the  thermospheric  wind 
system  by  magnctospheric  convection.  Roble  et  al.  (1984)43  showed  that  the  high-latitude  thermospheric 
temper,  ore  and  nruxal  composition  were  well-ordered  in  geomagnetic  coordinates  and  were  able  to  obtain 
resoruci-  zg-  v  incni  between  the  TGCM  calculations  and  DE-2  measurements  of  temperature  and 
composition. 

These  TGCM  results,  as  well  as  those  from  other  theoretical  and  numerical^  "*8  studies  have  all  shown 
that  the  neutral  winds  at  F-rcgion  altitudes  are  very  effectively  controlled  by  the  ion-drag  momentum  source 
associated  with  magnctospheric  convection.  The  twin-cell  pattern  of  magnctospheric  convection  drives  a 
largely  rotational,  non  divergent,  double-vortex  wind  system  at  F-region  heights.  The  magnitude  of  the 
neutral  wind  is  related  to  the  strength  of  the  forcing  which,  in  him,  depends  upon  both  the  ion-drift  speeds 
and  the  ion  and  neutral  number  densities.  The  TGCMs  have  shown  that  the  neutral  wind  pattern  effectively 
follows,  but  generally  lags,  the  ion-drift  pattern  and  has  variations  similar  to  those  associated  with  the  ion- 
drift  pattern.  In  general,  the  TGCM  simulations  have  indicated  that,  if  the  imposed  ion-drift  pattern  has 
sharp  boundaries  and  is  relatively  steady  for  time  periods  on  the  order  of  1/2  to  1  hour,  then  the  neutral  wind 
pattern  will  develop  similar  features. 

In  a  detailed  recent  study  using  the  NCAR-TGCM4^,  an  analysis  has  been  made  of  data  from  of  a  number 
of  orbits  of  the  DE-2  spacecraft  during  21-22  November  1981,  which  was  a  period  of  relatively  steady 
geomagnetic  forcing.  This  study  demonstrated  a  reasonably  good  agreement  between  the  DE-2 
measurements  of  auroral  particle  precipitation  and  ion-drift  velocity  vectors  and  the  TGCM  auroral 
parameterizations  for  these  parameters.  Large-scale  structures  in  the  thermospheric  neutral  wind  field  were 
reasonably  well  reproduced,  but  the  DE-2  measurements  of  neutral  winds  showed  considerably  more  small- 
scale  structure  than  the  model  could  simulate,  because  of  limitations  imposed  through  use  of  the  five  degree 
latitude-by-longitude  TGCM  grid.  Fig.  12  shows  a  comparison  of  the  TGCM-predicted  wind  and 
temperature  structure  of  the  polar  cap  along  with  DE-2  measurement*  of  the  cross-track  wind  component 
along  the  satellite  track  for  six  DE-2  orbits.  The  results  show  the  influence  of  magnctospheric  forcing  on 
the  high-latitude  circulation  pattern,  with  antisunward  winds  in  the  magnetic  polar  cap  and  sunward  winds 
on  the  dusk  side  of  the  auroral  oval  The  TGCM-predicted  wind  velocities,  reversal  boundaries,  and  the 
transition  from  the  mid-latitude  solar  EUV  and  UV -driven  wind  system  to  that  controlled  by  magnctospheric 
forcing  at  high  latitudes  are  all  in  reasonable  agreement  with  DE-2  measurements  for  this  period  of 
relatively  steady  geomagnetic  forcing.  This  study  represents  the  first  time  that  an  attempt  has  been  made  to 
simulate  numerically  the  time-dependent  thermospheric  response  to  real  magnetospheric  and  solar  forcings. 
The  encouraging  results  point  to  the  need  for  additional  future  work  to  be  conducted  along  these  lines  to 
investigate  the  perturbed  thermosphere  for  periods  when  DE-2  measurements  provide  both  the  necessary 
information  on  inputs  to  the  thermosphere  as  welt  as  a  detailed  record  of  the  large-scale  thermospheric 
response  to  the  time-dependent  forcing  mechanisms. 

Further  studies  investigating  the  circulation,  temperature,  and  compositional  structure  of  the  winter 
hemisphere  polar  cap  have  been  carried  oul^*^.  The  DE-2  measurements  show  considerable  seasonal 
differences  in  all  neutral  atmospheric  parameters.  In  the  summer  hemisphere,  solar  EUV  and  UV 
photoionization,  photodissociation,  and  neutral  heating  generally  dominate  auroral  inputs  for  all  except  the 


most  disturbed  periods.  In  the  winter  hemisphere,  the  diminution  of  these  solar  radiative  inputs,  because  of 
the  relatively  large  solar  zenith  angles,  makes  auroral  inputs  more  important  in  the  maintenance  of  the 
circulation,  temperature  and  compositional  structure  of  the  winter  polar  cap.  Rees  et  aL  (1985)4^  have 
shown  that  the  region  of  the  summer  pole,  under  similar  geomagnetic  and  solar  activity  conditions  at  a 
level  of  about  300km,  and  at  a  similar  level  of  geomagnetic  and  solar  activity  conditions,  is  about  300  K 
wanner  than  that  of  the  winter  pole.  The  density  of  atomic  oxygen  is  strongly  depleted  (and  nitrogen 
enhanced)  around  the  summer  pole  compared  to  the  winter  pole.  They  also  point  out  that  only  part  of  the 
differences  in  temperature  and  compositional  structure  can  be  related  to  the  seasonal  variation  of  solar 
insolation  and  that  both  polar  regions  display  structural  variations  (with  latitude  and  UT)  that  are 
characteristic  of  strong  magnetospheric  forcing.  The  magnetospheric  control  is  strong  in  the  summer 
hemisphere  because  of  the  enhanced  ion  densities  caused  by  solar  photo  ionization  and  weaker  in  the  winter 
hemisphere  where  this  coupling  is  somewhat  smaller. 

Killeen  and  Roble  (l^)4^  have  developed  a  diagnostic  package  to  examine  the  balance  of  forces  within  the 
NCAR-TGCM.  Their  analysis  showed  that  at  F-region  altitudes  the  largest  forces  in  the  “steady  state"  are 
due  to  ion  drag  induced  by  sunward  drifting  ions  on  the  dawnside  and  deskside  of  the  auroral  ovaVpolar  cap 
boundary.  At  F-rcgion  altitudes,  the  basic  balance  of  forces  is  between  ion  drag  and  pressure,  whereas  at  E- 
region  altitudes,  Coriolis  and  advectioo  forces  also  have  significant  contributions.  Later,  these  authors 
extended  their  diagnostic  package  to  include  neutral  parcel  trajectories  in  the  high-latitude  E-  and  F- 
rcgions^.  They  found  that  trajectories  analyzed  for  'steady  state'  diumally-reproducible  runs  of  the  TGCM 
show  the  following  features:  1)  the  magnetospheric  convection  throat  region  plays  an  important  regulatory 
role  for  the  overall  flow  pattern  at  high  la  tinges,  2)  the  trajectories  of  typical  F-region  neutral  gas  parcels 
pass  through  the  throat  region  at  least  once  and  sometimes  twice  in  a  24-hour  period  as  a  result  of  the 
convergent  geometry  of  the  ion  convection  pattern,  3)  parcels  transiting  the  dayside  cusp  region  are  heated 
rapidly  in  the  region  of  soft  particle  precipitation  and  this  energy  is  adverted  over  the  polar  cap  and  deposited 
several  hours  'downstream'  along  the  parcel  trajectory  beating  the  polar  cap,  4)  the  'postmidnight*  surge  in 
the  meridional  neutral  wind  if  the  result  of  parcel  acceleration  in  the  convection  throat  region  and  subsequent 
deflection  towards  the  corotation  direction,  and  5)  E-region  parcels  move  less  rapidly  than  their  F-region 
counterparts  and  can  dwell  for  appreciable  periods  of  time  in  the  auroral  oval;  only  a  few  geometrically- 
selected  parcels  escape  to  lower  latitudes. 

Finally,  a  spectral  analysis  capability  to  be  used  with  the  NCAR-TGCM  as  part  of  its  diagnostic  package 
has  been  developed^2.  This  capability  allows  the  gtidded  global  TGCM  output  to  be  decomposed  according 
to  a  series  of  vector  spherical  harmonic  functions,  enabling  the  power  in  the  various  wave  modes 
transmitted  though  the  TGCM  grid  to  be  intercom  pared  quantitatively.  The  expansion  of  the  wind  field  in 
vector  spherical  harmonics,  combined  with  a  fourier  expansion  in  time  and  a  polynomial  expansion  in 
altitude,  has  enabled  the  time-dependent,  3-D,  diurnal  neutral  wind,  ion  drift,  neutral  temperature  and 
composition  fields  to  be  reproduced  using  a  simple  and  fast  computer  subroutine.  This  code,  termed  the 
VSH  (sector  spherical  harmonic)  model,  is  currently  being  used  as  the  basis  for  an  effort  to  produce  a 
Dynamics  Explorer  semi -empirical  model  of  global  thermosphere  dynamics.  A  comparison  of  the  global 
output  irom  the  VSH  model  (currently,  simply  a  replication  of  the  NCAR-TGCM  predictions)  with  neutral 
wind  measurements  made  on  single  orbital  passes  of  DE  2  were  reported-*2.  The  agreement  between  the 
TGCM  predictions  and  the  experimental  data  was  very  good  for  both  high-latitudc  regions,  giving 
confidence  in  the  TGCM  itself  as  well  as  in  the  future  use  of  the  VSH  representation  to  yield  a  realistic  and 
useful  semi-empirical  model  of  thermosphere  dynamics.  Such  a  model  will  represent  one  of  the  more 
significant  legacies  of  the  DE  mission. 

CONCLUDING  REMARKS 

Considerable  progress  has  been  sustained  ov<t  the  past  5  years  though  use  of  the  Dynamics  Explorer  data 
base  to  improve  our  understanding  of  the  dynamics  of  the  Earth’s  upper  atmosphere.  The  combination  of 
comprehensive  experimental  coverage  of  the  key  parameters  involved  in  magnetosphere-ionosphere- 
thermosphere  coupling  together  with  improved  theoretical  techniques  has  led  to  important  advances  in  both 
the  morphological  description  of  thermosphere  dynamics  and  the  identification  and  quantification  of  the 
main  forcing  mechanisms. 
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Ia  the  area  of  uKxpbological  studies,  the  DE  mission  has  provided  information  on  the  mean  circulation  at 
F-rcgioo  altitudes  for  both  high-  and  low-latitudes.  The  high-latitude  thermospheric  wind  pattern  has  been 
shown  tc  have  definite  dependencies  on  UT.  LT,  season,  IMP  orientation,  and  geomagnetic  activity  level 
and  progress  has  been  made  in  quantifying  these  dependencies  for  the  solar  maximum  conditions  appropriate 
to  the  DE-2  mission  lifetime.  The  extent  to  which  the  equatorial  thermosphere  co-rotatrs  or  super-rotates 
with  the  Earth  has  been  investigated,  leading  to  a  satisfactory  agreement  between  theory  and  experiment. 
This  experimental  result,  indicating  that  very  little  supperrotation  occurs  in  the  icnestnil  thermosphere,  has 
fundamental  implications  for  the  study  of  the  dynamical  aerooomy  of  other  planets  wlere  superrotation  is 
known  to  be  an  important  feature.  DE  measurements  have  been  used  to  investigate  th:  detailed  coupling 
between  the  ionosphere  and  the  neutral  thermosphere.  The  strong  influence  of  energy  and  momentum 
mapping  down  from  the  magnetosphere  into  the  ionosphere  and  influencing  thermospheric  circulation  and 
composition  has  been  demonstrated,  and  the  inputs  to  the  theoretical  models  improved  through  comparative 
studies.  Complex  feedback  processes  acting  to  modulate  thermosphere-ionosphere  coupling  have  been 
identified  and  studied.  The  role  of  ion-drag  momentum  forcing  and  the  Coriolis  force  in  establishing  the 
thermospheric  dreuiatioo  has  been  elur  dated. 

Perhaps  the  greatest  benefit  of  the  DE-2  measurements  has  been  to  provide  the  evolving  theoretical  models 
of  the  thermosphere  with  rigorous  and  global  experimental  constraints.  The  close  and  interative  interaction 
between  experimentalists  and  theorists  working  with  DE  data  has  led  to  great  improvements  in  the 
theoretical  understanding  of  the  region,  leading  to  a  situation  where  a  predictive  capability  has  become  a  real 
possibility  and  a  goal  for  future  work.  This  interaction  is  of  particular  importance  for  die  improvement '  f 
predictions  of  atmospheric  density  fluctuations  from  first-principal  calculations.  Because  of  the  strong 
coupling  between  composition,  atmospheric  density  perturbations  and  dynamics,  a  full  and  quantitative 
understanding  of  the  thermospheric  wind  systems  and  their  influence  on  other  atmospheric  parameters  is 
required  to  constrain  and  test  sophisticated  models  such  as  the  TGCMs. 
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TABLES 


Table  L 

Relationship  between  various  characteristic  features  of  the  neutral  thermospheric  winds  in  the  polar  regions 
and  Kp  for  the  specified  regression  line.  Taken  from  RTef.  29. 

labkH 

Same  as  Table  I  except  that  the  relationships  are  for  the  Auroral  Electrojet  index  AE. 
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FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3 


Figure  4. 


Figure  5. 


Figure  6. 


Typical  measurement  geometry  for  the  DE-FPI  (bottom  schematic).  The  instrument 
performed  a  limb  scan  of  the  horizon,  with  scan  angles  varying  between  5°  and  15°  from 
the  local  horizontal  (shown  here  as  12°).  The  instrument  measured  the  Doppler  shift  of  the 
0(*D)  emission.  A  typical  volume  emission  rate  profile  for  the  is  emission  is  shown  at  top 
right,  together  with  a  weighting  function  (labelled  'path  length*)  that  indicates  the 
geometrical  path  length  through  the  atmosphere  for  the  given  viewing  geometry  shown 
below.  The  contribution  function  (top  right)  is  simply  given  by  multiplying  the  local 
volume  emission  rate  by  the  path  length  and  represents  the  contribution  to  the  given  line- 
of-sight  measurement  from  the  various  altitude  layers.  As  can  be  inferred  from  the 
sharpness  of  the  contribution  function,  the  principal  contribution  for  any  given 
measurements  comes  bom  the  tangent  point  altitude  along  the  line  of  sight 

Schematic  for  the  DE  WATS  measurement  of  wind  and  temperature  (Ret  6).  The  apenure 
of  the  dose-source  mass  spectrometer  is  scanned  systematically  in  the  vertical  and 
horizontal  planes  by  the  two  baffles  shown  in  Plate  3.  As  a  baffle  scans  in  front  of  the 
aperture,  the  signal  counts  drop  with  a  characteristic  signature  as  shown.  The  baffle 
location  at  the  time  of  maximum  signal  drop-off  is  related  to  the  wind  component  in  the 
plane  of  the  baffle  motion  and  this  wind  can  be  obtained  through  solution  of  the  vector 
triangle  shown.  This  measurement  requires  precise  knowledge  of  both  the  spacecraft  speed 
and  direction.  The  shape  of  the  signal  drop  off  is  governed  by  the  velocity  distribution  of 
the  neutral  gas  particles  and  hence  is  related  directly  to  the  kinetic  temperature. 

Horizontal  neutral  wind  vectors  at  -350  km  altitude  from  merged  FPI  and  WATS 
measurements,  a)  orbital  passes  over  the  Southern  Hemisphere  polar  cap;  b)  orbital  passes 
over  the  Northern  Hemisphere  polar  cap.  The  winds  are  given  by  the  arrows,  with  the 
arrow  length  plotted  according  to  the  scale  at  lower  right. 

Climatological  map  for  UT  bins  centered  at  07:00,  09:00,  and  ll:00hrs  UT  (top  to 
bottom).  The  left-hand  side  shows  the  satellite-  and  GBFPI-averaged  measurements  plotted 
in  geographic  polar  coordinates  (latitude  and  LT).  The  satellite-averaged  winds  are  given  by 
the  arrows  according  to  the  wind  scale  at  lower  right  The  GBFPI  wind  measurements  are 
plotted  as  standard  meteorological  symbols  with  a  barb,  100m/scc;  long  line,  SO  m/sec; 
short  line,  10  m/sec.  The  curved  line  is  the  solar  terminator.  The  left  hand  side  illustrates 
theNCAR-TGCM  predictions  for  the  midpoint  of  the  particular  UT  bin.  Model  winds  are 
plotted  according  to  the  same  scale  as  the  satellite  winds.  Figure  taken  from  Ret  15. 

Avenged  thermospheric  neutral  wind  measurements  for  a)  north  pole  By  positive;  b)  north 
pole  By  negative;  c)  south  pole  By  positive;  d)  south  pole  By  negative.  Data  collected 
between  November  1981  and  January  1982  and  between  November  1982  and  January  1983 
were  averaged  according  to  the  sign  of  the  By  component  of  the  IMF  and  plotted  in 
geomagnetic  polar  coordinates  (magnetic  latitude  and  magnetic  local  time).  Outer  circle  of 
each  polar  dial  is  at  40°  geomagnetic  latitude.  The  local  time  orientation  of  the  polar  dials 
and  vector  winds  is  given  at  top  right.  Figure  taken  from  Ref 26.. 

Averaged  thermospheric  neutral  wind  measurements  for  a)  south  pole  high  Kp;  b)  south 
pole  low  Kp;  c)  north  pole  high  Kp;  d)  north  pole  low  Kp.  Data  collected  between 
November  1981  and  January  1982  and  between  November  1982  and  January  1983  were 
averaged  according  to  the  given  range  of  Kp  and  plotted  in  geomagnetic  polar  coordinates 
(magnetic  latitude  and  magnetic  local  time).  Outer  circle  of  each  polar  dial  is  at  40° 
geomagnetic  latitude.  The  local  time  orientation  of  the  polar  dials  and  vector  winds  is 
given  at  top  right  Figure  taken  from  Ref.29„ 
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Figure  7.  Dependence  of  the  maximum  geomagnetic  polar  cap  ami-son  ward  neutral  wind  as  a 
function  of  the  geomagnetic  indices  Kp  (top)  and  >E  (bottom).  Figures  1  la)  and  c)  ux  for 
north  pole  observations  and  1  lb)  and  d)  are  for  south  pole  observations.  The  expressions 
given  represent  fits  to  the  experimental  data.  Wind  observations  over  6  months  of  the  DE-2 
mission  were  used  to  provide  these  relationships.  Figure  taken  from  Ref.  29. 

Feme  8.  Zonal  neutral  wind  observations  from  WATS  for  two  northern  hemisphere  high-latitude 
orbital  passes  are  compared  with  simultaneous  zonal  ion  drift  observations  from  IDM.  Top 
panel  for  orbit  7210.  bottom  panel  for  orbit  7212.  Data  are  plotted  at  the  same  temporal 
and  spatial  resolution  u)  provide  an  indication  of  the  integrating  effect  of  the  neutral 
atmosphere  (note  greater  spatial  structure  evident  in  ion  drift  data,  see  text).  Figure  taken 
from  Ref .27.. 

Figure  9.  a)  Simultaneously  measured  neutral  wind  vectors  from  DE  2  (orbit  1813)  and  DE-1  auroral 
image  (courtesy  L.  A.  Frank  and  J.  D.  Craven,  Universtiy  of  Iowa).  The  images  were 
obtained  using  the  SAI  instrument  viewing  at  ultraviolet  wavelengths.  The  im~g**  «re 
false-color  coded  and  are  oriented  such  that  the  direction  towards  the  sun  is  to  the  top  of  the 
figure,  dusk  to  the  left.  The  solar  terminator  is  evident,  running  roughly  horizontally 
across  each  image,  as  is  the  entire  auroral  oval  located  just  to  the  nightside  of  the 
terminator.  The  neutral  wind  vectors  are  denoted  by  the  yellow  arrows  whose  origins  are 
positioned  along  the  DE-2  orbital  track.  The  wind  scale  is  given  at  lower  right  9b,  9c,  and 
9d  are  for  DE-2  orbits  1815, 1819,  and  1823,  respectively.  Figure  taken  from  Ref33_ 

Figure  10.  Zonal  wind  average  measurements  from  WATS  at  equatorial  latitudes  in  the  altitnde  range 
300  -  500  km.  Table  shows  the  numerical  results  of  averaging  the  zonal  wind 
measurements  according  to  latitude  interval  and  altitude  range.  Figure  taken  from  Ref37„ 

Figure  11.  Comparison  of  neutral  wind  measurements  from  Dynamics  Explorer  with  the  predictions 
of  the  UCL-TGCM  for  three  dawn-dusk  crossings  of  the  spacecraft  over  the  Northern 
Hemisphere  polar  region.  The  neutral-wind  measurements  are  given  by  the  vectors  plotted 
along  the  latitudinal  track  of  the  spacecraft  (labelled  KVEC,  with  orbit  number)  and  the 
UCL-TGCM  predictions  for  two  model  runs  corresponding  to  different  ion  density  input 
models  are  given  by  the  vectors  plotted  to  the  top  and  bottom  of  the  experimental 
measurements.  Figure  taken  from  Ref.41. 

Figure  12.  Zonal  neutral  wind  component  measured  by  the  WATS  instrument  along  the  DE-2  satellite 
track  over  the  polar  caps  for  six  orbits:  a)  1639.  b)  1640.  c)  1641,  d)  1641,  e)  1644,  and  f) 
1646.  All  the  orbits,  with  the  exception  of  d)  are  for  the  Northern  Hemisphere.  Also 
shown  are  the  TGCM-predkted  neutral  temperature  pattern  (K)  and  neutral  wind  vectors 
along  the  model  z=m-2  constant-pressure  surface  (-400  km)  over  the  polar  caps  at  the  times 
of  the  DE-2  satellite  passes.  The  outer  circle  represents  a  latitude  of  either  ±373°  latitude 
depending  on  hemisphere.  Figure  taken  from  Ref.  48. 
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AFGL  LIDAR  FACILITY 
Dwight:  P.  Sipler 
Force  Cecphyuics  Laboratory 


Lidar  is  an  acronym,  Light  Detection  And  Ranging.  This  is  the  same  acronym  as  radar, 
except  that  light  is  used  instead  of  radio  waves.  Other  than  that,  the  basic  principles  are 
similar,  i.e.  a  short  pulse  of  light  is  transmitted,  hits  a  scattering  object  and  returns  to  the 
detector.  The  range  is  determined  by  the  delay  time  between  transmit  and  receive.  A  very 
simplified  description  of  an  atmospheric  sounding  lidar  is  given  by  the  equation 

i.lj  S  =  P0  K  n  r'2 

where  S  is  the  detected  signal,  Pu  is  the  transmitted  power,  K  is  a  constant  containing  system 
parameters,  n  is  the  number  of  scattering  objects  and  r  is  the  range  of  the  objects.  In  a 
radar  system  the  range  dependence  is  usually  written  as  r“*.  but  that  is  for  a  target  which  is 
small  compared  to  the  dimensions  of  the  transmitted  beam  so  that  the  illumination  of  the 
target  depends  on  r'“  and  the  light  scattered  from  the  target  into  the  detector  is  also 
proportional  to  r'~.  Since  we  are  looking  at  the  atmosphere,  the  target  completely  fills  the 
beam  and  only  the  light  scattered  into  the  detector  gives  a  range  dependence.  The  number  of 
scattering  objects  is  just  the  density  of  the  atmosphere.  The  catchall  parameter  K  includes 
such  things  as  (he  optical  system  efficiency  (both  transmit  and  receive),  extinction  of  the 
light  by  the  atmosphere  (mostly  dust  and  aerosols  in  the  lower  atmosphere)  and  the  cross 
section  of  the  atmospheric  molecules  (basically  a  constant  up  to  around  90  km). 

Several  different  scattering  processes  are  important  to  lidar  measurements  of  atmospheric 
density: 

Rayleigh  scattering 
Mle  scattering 
Raman  scattering 
Resonance  scattering. 

Rayleigh  scattering  is  the  scattering  of  light  by  small  particles  and  is  familiar  to  most 
people  as  the  process  which  causes  the  blue  color  of  the  sky.  Rayleigh  scattering  theory 
describes  the  scattering  of  electromagnetic  radiation  by  particles  which  arc  small  compared  to 
the  wavelength  of  the  radiation.  In  the  case  of  atmospheric  sounding,  the  scattering  particles 
are  the  atoms  and  molecules  which  make  up  the  atmosphere.  The  scattering  cross  section  is 
relatively  small  and  varies  as 

Mie  scattering  theory  describes  the  scattering  of  electromagnetic  radiation  by  particles 
which  are  of  the  order  of  the  wavelength  of  the  radiation  or  larger.  This  process  is 
important  in  the  lower  regions  of  the  atmosphere  (from  the  ground  to  about  25  km)  where 
dusts  and  aerosols  are  significant  components  of  the  atmosphere.  The  wavelength  dependence 
of  this  type  of  scattering  is  a  complex  function  of  particle  size  and  wavelength  but  in  general 
the  scattering  cross  section  is  larger  than  the  Rayleigh  cross  section  and  varies  as  A'1. 

Raman  scattering  is  similar  to  Rayleigh  scattering  except  that  where  Rayleigh  scattering 
is  an  elastic  process,  Raman  scattering  is  inelastic  scattering.  This  means  that  a  photon, 
when  scattering  from  an  atom  or  molecule,  leaves  some  of  its  energy  behind  in  the  atom  or 
molecule.  Since  only  discrete  energy  levels  are  available  and  these  are  different  for  different 
atmospheric  components,  the  energy  loss  will  be  characteristic  of  the  particle  doing  the 
scattering.  The  energy  loss  means  that  the  wavelength  of  the  scattered  light  will  be  longer 
(lower  energy)  than  the  light  from  the  laser.  Thus,  by  selecting  a  particular  wavelength  shift 
with  an  interference  filter,  measurements  can  be  made  on  specific  components  of  the 
atmosphere,  e  g.  N2,  C>2  °r  H2O.  The  Raman  scattering  cross  section  is  very  small,  about  10  ^ 
of  the  Rayleigh  cross  section.  This  limits  the  utility  of  this  process  to  lower  altitudes  where 
the  light  from  the  laser  is  most  intense.  However,  this  is  just  where  this  sort  of 
measurement  is  needed,  since  the  scattering  is  relatively  insensitive  to  the  presence  of  dust 
and  aerosols. 
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Resonance  scattering  is  another  process  which  depends  on  specific  energy  levels  in  a 
scattering  atom  or  molecule.  If  the  laser  can  be  tuned  to  a  wavelength  which  represents  the 
exact  energy  needed  to  raise  the  target  atom  to  an  excited  state  from  the  ground  state 
through  an  allowed  transition  then  the  cross  section  for  absorption  and  reradiation  of  the 
light  will  be  very  large.  This  large  cross  section  makes  this  process  suitable  for  the 
measurement  of  minor  components  of  the  atmosphere,  e  g.  free  atomic  sodium,  which  occurs  in 
a  thin  layer  between  approximately  85  and  105  km  altitudes.  The  sodium  density  in  this 
region  is  on  the  order  of  10^  cm'^  compared  to  the  ambient  atmospheric  density  of  around 
10^  cm'3.  However,  the  resonance  cross  section  is  about  11  orders  of  magnitude  larger  than 
the  Rayleigh  cross  section,  so  measurements  of  sodium  can  be  performed  with  moderate  laser 
power.  Other  workers  have  used  this  same  technique  to  look  at  calcium,  potassium,  lithium 
and  iron.  The  technique  could  also  be  used  on  ions  as  well  as  neutrals. 


Laboratory  Lldar  System 

AFGL  has  two  lidar  systems  for  mesospheric  sounding:  a  laboratory  system  for  testbed 
applications  and  a  mobile  version  for  field  work.  The  laboratory  system  is  located  at  AFGL 
on  the  roof  of  building  1102F.  The  system  consists  of  a  laser  transmitter,  beam  expander  and 
steering  optics,  a  telescope  detector  and  a  data  acquisition  system.  The  basic  laser  is  a 
Quantel  model  582  Nd:YAG  laser  with  doubling  and  tripling  crystals.  This  provides  us  with 
three  wavelengths:  1064  nm  (fundamental),  532  nm  (doubled)  and  355  nm  (tripled).  The  532 
and  355  nm  wavelengths  are  used  for  sounding  and  the  1064  nm  is  discarded  since  the 
Rayleigh  scattering  at  this  wavelength  is  small  and  we  have  no  good  detectors  for  1064  nm 
light.  The  use  of  two  wavelengths  enables  us  to  monitor  the  wavelength  dependence  of  the 
scattering  profiles  and  thus  to  identify  regions  in  which  Mie  scattering  from  dusts  and 
aerosols  is  significant.  We  do  not  currently  use  the  data  from  these  regions,  but  if  the 

particle  size  distribution  could  be  determined  or  accurately  modelled,  a  determination  of  the 
particle  density  could  be  extracted  from  the  two  color  returns.  The  laser  output  beam  is  1 
cm  in  diameter  and  has  a  0.6  mrad  beam  divergence.  This  beam  is  expanded  by  a  factor  of 
15,  resulting  in  a  15  cm  diameter  team  with  0.04  mrad  divergence.  The  expanded  beam  is 
directed  vertically  by  an  adjustable  45°  mirror. 

The  telescope  detector  consists  of  a  36  inch  Dall-Kirkham  telescope  coupled  to  a 

6  channel  detector  with  a  mechanical  shutter.  The  shutter  runs  at  7000  rpm  and  has  3 

apertures,  generating  350  openings  per  second.  Timing  for  the  lidar  system  is  taken  from  the 
shutter  to  ensure  synchronization.  The  laser  is  fired  before  the  shutter  opens  so  that  the 
light  scattered  in  the  troposphere  is  blocked  from  the  photomultipliers.  Since  most  of  the 
dusts  and  aerosols  in  the  atmosphere  are  located  below  the  tropopausc  and  since  r'^  is  large 
in  this  region,  the  shutter  blocks  the  major  part  of  the  returning  light.  This  is  done  to  avoid 
saturation  and  recovery  effects  in  the  photomultipliers.  After  passing  through  the  shutter,  the 
light  is  split  into  a  532  nm  beam  and  a  355  nm  beam  by  a  dichroic  beamsplitter  and  further 
split  into  three  photomultipliers  for  each  wavelength.  These  photomultipliers  have  different 
sensitivities  to  split  the  required  dynamic  range  of  the  instrument  into  smaller  ranges.  A 

second  shutter  blocks  the  highest  sensitivity  photomultiplier  for  a  time  equivalent  to  an 
additional  15  km  of  range  in  order  to  reduce  overload  from  the  high  signal  levels  from  the 
lower  altitudes.  Since  the  shutter  opens  at  350  Hz  and  the  laser  runs  at  10  Hz  several 
shutter  openings  arc  available  for  background  light  monitoring.  The  background  data  are 
separated  by  software. 

The  data  acquisition  system  is  currently  a  LeCroy  3500M  microcomputer  system  with  an 
integrated  camac  crate.  Data  from  the  6  photomultipliers  is  multiplexed  and  fed  into  two 
multichannel  scaling  modules  in  the  camac  crate.  Data  from  the  scaling  modules  is  integrated 
in  the  computer  through  software.  Since  the  read  cycle  time  for  the  camac  modules  is  1  psec 
the  maximum  sampling  rate  for  the  system  is  2  psec  per  point,  one  psec  for  each  module. 
Acquired  data  are  stored  on  floppy  disks  for  later  analysis.  We  are  currently  installing  a  new 
80386  based  PC  coupled  to  a  camac  crate  which  will  both  simplify  data  acquisition  and  enable 
some  real-time  analysis  of  the  data. 
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AFGL  is  located  under  the  tanding  pattern  for  Hanscom  field  and  is  close  to  the  flyway 
into  Logan  airport.  Interception  of  the  laser  beam  by  aircraft  poses  a  safety  hazard  due  to 
the  high  power  in  the  beam.  fn  order  to  avoid  such  problems  we  have  installed  an 
AN/APS  42  radar  system  co-aligned  with  the  laser  beam.  If  an  aircraft  is  detected  by  the 
radar  the  laser  is  automatically  shut  down  and  requires  manual  intervention  to  restart  the 
system. 

In  addition  to  the  Nd:YAG  laser,  a  Lambda  Physik  FL2002LC  narrow  band  dye  laser  is 
available  for  use  in  resonance  scattering  measurements.  This  laser  has  been  used  for 
measurements  on  the  atomic  sodium  layer  in  the  85-105  km  region.  The  ultimate  bandwidth  of 
this  laser  is  narrower  than  the  doppler  width  of  the  sodium  D2  line  from  this  region,  so  we 
will  be  able  to  use  it  for  measurements  of  the  width  of  the  resonance  line.  This  will  give  us 
a  measurement  of  the  ambient  temperature  in  this  region. 

A  second  telescope  detector  is  included  in  the  laboratory  system  for  Raman  scattering 
measurements.  This  telescope  is  a  12  inch  telescope  and  is  provided  with  two  photomultipliers 
for  detection  of  the  low  altitude  scattering  from  the  laser  beam.  One  channel  is  used  for  the 
normal  Raylcigh/Mie  scattering  signal  and  the  other  is  used  for  a  Raman  signal.  An 
interference  filter  is  used  which  has  a  passband  at  the  Raman  shifted  wavelength  but  also  has 
a  strong  band  rejection  around  the  laser  wavelength  to  avoid  contamination  of  the  weak 
Raman  signal  by  the  strong  Rayleigh/Mie  return. 


Mobile  LIdar  System 

The  mobile  system  is  similar  to  the  laboratory  system  except  that  it  is  scaled  down  a  bit 
to  get  it  into  a  small  trailer.  The  trailer  is  equipped  with  jacks  to  enable  removal  of  the 
wheels  for  air  shipment.  The  trailer  has  been  shipped  on  a  C141  to  Alaska  and  back.  The 
wheels  are  suitable  for  short  movements  of  the  trailer  system  but  not  for  extended  travel 
For  long  overland  travel  the  trailer  is  set  onto  a  flatbed  trailer. 

The  laser  in  the  mobile  system  is  a  Quantel  model  581C  Nd:YAG  laser  which  generates 
1.7  Joule  pulses  instead  of  the  2  Joule  pulses  of  the  laboratory  system.  The  telescope  is  a  30 
cm  telescope  instead  of  the  96  cm  telescope  of  the  laboratory  system.  The  telescope  is 

mounted  horizontally  on  an  optical  table  and  a  45°  mirror  is  located  in  front  of  the  telescope 

to  direct  the  field  of  view.  The  center  section  of  the  45°  mirror  is  cut  out  and  replaced  with 

a  dielectric  mirror  for  reflection  of  the  laser  beam.  This  arrangement  avoids  anv  possible 
problem  with  fluorescence  of  the  turning  nunor  iiuiu  me  ia^er  beam  since  the  light  from  the 
center  of  the  turning  mirror  is  blocked  from  the  rest  of  the  mirror  by  a  metal  mounting  ring. 
Scattered  light  from  this  section  is  blocked  by  the  secondary  of  the  telescope  so  fluorescence 
in  the  laser  turning  mirror  does  not  reach  the  detector.  The  45®  turning  mirror  can  be  rotated 
about  the  telescope  axis  to  direct  both  the  laser  beam  and  the  detector  field  of  view  toward  a 
selected  elevation  angle  from  slightly  below  the  horizon  to  20°  past  the  zenith.  Small  azimuth 
adjustments  are  made  by  moving  the  mirror  mount  about  a  vertical  axis. 

The  detector  is  equipped  with  a 
shutter  to  block  out  the  light  from  low 

altitudes.  However,  instead  of  the 
multiplexing  system  used  in  the  laboratory 
system,  the  mobil<*  shutter  is  designed  to 
extend  the  dynamic  range  of  the  detector 
by  acting  as  a  variable  aperture  in  the 
telescope.  The  primary  of  the  telescope  is 
focussed  onto  the  shutter  wheel,  which  has 
an  arrow  shaped  aperture  (see  fig.  1). 

This  shutter  opening  moves  slowly  across 
the  image  of  the  telescope  primary  and 
effectively  opens  the  aperture  slowly, 
taking  the  equivalent  of  about  20  km  of 
range  from  first  light  to  full  open.  Since 


Fig.  1:  Variable  Shutter 
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the  shutter  is  moving  across  the  image  of 
the  primary  mirror  there  is  no  time 
modulation  of  the  field  of  view  of  the 
telescope.  This  avoids  possible  timing 
changes  due  to  changes  in  the  position  of 
the  laser  return  within  the  field  of  view 
of  the  telescope.  The  shape  of  the 
transmission  function  can  be  determined  by 
looking  at  a  uniform  light  source  placed  in 
front  of  the  telescope.  An  illuminated 
reflector  can  be  used  but  in  practice,  the 
day  sky  is  the  most  uniform  source 
available  since  cloud  motions  are  averaged  in  time, 
placing  a  neutral  density  filter  at  the  primary  focal  point  of  the  telescope, 
transmission  function  is  shown  in  fig.  2. 


Fig.  2:  Shutter  Transmission  Function 

The  brightness  of  the  sky  is  reduced  by 

A  typical 


A  62  cm  detector  telescope  is  also  located  in  the  trailer  to  extend  the  altitude  range  of 
the  lidar  system.  This  telescope  is  fixed  vertically  rather  than  being  able  to  be  directed  at 
different  elevation  angles.  In  use,  the  45®  turning  mirror  of  the  laser  is  directed  vertically 
and  small  adjustments  in  the  mirror  mount  are  made  until  the  laser  beam  and  the  fixed 
telescope  field  of  view  overlap.  This  detector  includes  a  fast-opening  shutter  which  is  phase 
locked  to  the  shutter  of  the  primary  system  by  a  servo  system.  Some  jitter  is  inevitable  in 
this  system,  making  the  opening  of  this  system  somewhat  ragged.  The  initial  data  from  this 
d elector  can  be  corrected  to  some  extent  by  modelling  the  jitter  and  measuring  the  time- 
averaged  response  of  this  detector  to  a  uniform  white  light  (again  using  daylight  appears  to 
be  most  effective). 


A  Lambda  Physik  EMG201MSC  excimer  laser  is  being  installed  in  the  mobile  lidar  system 
for  use  in  Rayleigh  scattering  at  higher  altitudes.  The  excimer  will  be  used  at  351  nm  and  is 
expected  to  generate  250  mJoule  pulses  at  up  to  80  Hz.  This  should  extend  our  altitude 
measurement  range  by  about  10  km. 

The  62  cm  detector  is  being  equipped  with  an  etalon  which  will  be  used  as  a  very  narrow 
band  interference  filter.  This  will  cut  down  the  continuum  light  passing  through  the  filter 
during  daylight  operations.  We  estimate  that  the  system  will  be  able  to  measure  densities  up 
to  about  55  km  during  daylight  with  this  system. 


System  Operation 

Since  the  shutter  timing  is  essential  to  the  proper  operation  of  the  system,  a  sensor  is 
provided  which  generates  a  pulse  at  a  specified  time  before  the  shutter  opens.  The  laser 
firing  is  synchronized  with  this  signal  to  ensure  a  constant  relationship  between  the  shutter 
opening  and  the  lidar  range.  The  pulse  from  the  laser  is  detected  in  the  laboratory  system 
and  used  to  trigger  the  data  acquisition  system.  In  the  mobile  system,  a  pulse  is  sent  from 
the  control  electronics  to  the  pockels  cell  of  the  laser  to  synchronize  the  laser  firing  with 
the  shutter.  Pulses  from  the  photomultipliers  are  counted  by  the  data  acquisition  system  and 
a  sequence  of  shots  is  added  to  generate  a  lidar  profile.  The  profiles  are  stored  on  floppy 
disks  for  later  analysis. 

Tn  order  to  minimize  the  contribution  of  ambient  light  from  the  sky,  the  field  of  view  of 
the  detector  is  made  as  small  as  possible  consistent  with  the  size  of  the  laser  spot  in  the 
atmosphere.  In  order  to  align  the  laser  beam  and  the  detector  telescope,  a  rough  alignment 
can  be  done  by  observing  the  beam  through  the  telescope.  After  the  rough  alignment,  fine 
adjustments  are  made  while  taking  short  bursts  of  data.  This  allows  us  to  maximize  the 
signal. 
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Data  Analysts 

The  data  are  read  from  the  floppy  disks  onto  a  PC.  Corrections  are  made  for  the 
shutter  transmission  function  in  the  case  of  the  mobile  system.  The  multiplexed 
photomultiplier  signals  are  joined  to  make  a  smooth  profile  in  the  case  of  the  laboratory 
system.  The  resulting  profiles  are  then  multiplied  by  r^.  This  generates  a  curve  which  is 
proportional  to  the  density  of  the  atmosphere  (see  eqn.  1).  The  constant  of  proportionality 
must  be  determined  by  some  means,  since  the  extinction  in  the  lower  atmosphere  is  not 
known.  The  most  common  method  of  calibration  b  the  use  of  a  rawinsonde  profile,  which  b 
matched  at  the  corresponding  altitudes  to  the  lidar  profile.  The  use  of  Raman  scattering  to 
measure  the  low  altitude  N2  profile  b  expected  to  give  us  a  way  to  measure  the  density 
profile  closer  to  the  ground  and  thereby  enable  the  matching  of  the  profile  to  surface 
measurements. 

The  errors  in  the  profiles  are  determined  by  the  counting  statbtics  of  the  pulses  from 
the  photomultipliers.  These  pubes  are  random  and  we  acquire  as  many  pubes  as  possible,  so 
gaussian  statbtics  are  assumed.  In  order  to  determine  the  density  to  within  5%,  we  need  at 
least  400  counts.  Thb  counting  level  b  arrived  at  by  (1)  adding  lidar  shots  until  the  desired 
level  b  attained  and  (2)  smoothing  the  data  over  several  range  bins  to  increase  the  number  of 
counts  used  for  a  determination.  These  procedures  limit  the  temporal  and  spatial  resolution  of 
the  lidar.  At  present,  the  mobile  system  is  capable  of  5%  statbtics  at  a  range  of  70  k.u  with 
a  15  minute  (9000  shot)  integration  time  and  2.7  km  (9  point)  range  smoothing. 


Summary 

AFGL  has  developed  a  lidar  system  for  measurements  of  atmospheric  density  and 
temperature  to  high  altitudes.  The  system  has  been  used  in  a  campaign  of  measurements  at 
Poker  Flats,  Alaska  Current  efforts  are  under  way  to  improve  the  system  by  extending  the 
high  altitude  range  of  the  lidar  and  by  using  additional  technology  to  provide  calibration  of 
the  measurements  by  comparbon  with  ground  truth. 
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4.2  LIDAR  AND  kULKtr  M  K  A  i>  U  R  K  M  L  N  T  b  UK  ATMUbHHtKU, 

STRUCTURE  PROPERTIES 

C. R.  Phi ibr lck 

A  KUL/  L 1  U  Hanscom  APB,  flA  01720 

Workshop  on  AtmoBphenc  Density  and  Aerodynamic  Drag 
Models  tor  Air  Eorce  operations 
AKGL  20-22  October-  12a7 


ABSTRACT 

The  results  1 rom  a  mobile  lidar  have  been  used  to  derive 
density  and  temperature  proliles  between  2b  and  Bb  km  over 
central  Alaska  during  Pebruar  y  -  Apr  1 1  lOfcib.  These  results 
allow  the  temporal  changer-  ol  the  middle  atmosphere  t.o  be 
observed  on  the  scales  ol  minutes,  hours,  and  days  lor  the 
first,  time.  This  perspective  has  Jed  to  several  new  points  in 
oui  understanding  ol  atmospheric  dynamics.  The  rneusur  ements 
ol  meteorological  rockets  made  simultaneously  with  the  lidar 
measurements  have  been  studied  together  with  the  results  1  r  om 
previous  measurement  campaigns  (MAP-WINK,  TRADER,  EBC, . .  >  to 
explain  points  about  the  atmospheric  dynamics  based  on  the 
observed  structure  variations. 

1 NTRGDULl  ION 

During  the  period  between  February  and  April  lyflb  a  new  mobile 
lidar,  ill.]  NT  ( liround  -  based  Lidar  IN  vest  l  gat  ion  -  Trarispor  table  > , 
was  used  to  obtain  data  on  the  variations  which  occur  in  the 
high  latitude  region  ol  the  atmosphere  between  20  and  Ob  km. 
The  lidar  measurements  were  undertaken  to  provide  an  improved 
data  base  on  which  the  reentry  1  light  character isti cs  ol  the 
Space  Shuttle  descent  iron  a  polar  orbit  could  be  based.  The 
lidar-  cun  measure  proliles  with  improved  altitude  and  time 
resolution  compared  to  standard  rocket  techmgues.  This  is 
ampiortant  lor  better  understanding  ol  the  dynamical  processes 
ol  the  middle  atmosphere.  Data  was  obtained  on  2b  nights 
during  the  program.  Standard  meteorological  balloon  and 
rocket  payloads  were  used  to  make  measurements  and  provide  a 
oompar i eon  study  with  the  lidar  data.  More  than  a  thousand 
prelates  ol  atmospheric  density  were  obtained  with  the  lidar 
and  twenty  meteorological  rockets  were  launched. 

The  measurement  campaign  has  resulted  in  several  cone 1  us  i  ons : 
(1)  the  lidar  data  lias  provided  the  opportunity  to  observe 
small-scale  variations  and  the  background  granular  r  t  y  ol  the 
atmosphere,  (2)  planetary  waves  may  be  associated  with  more 
than  hall  o  1.  the  total  density  variation  m  the  winter  middle 
atmosphere,  cs)  density  and  temperature  variations  can  be  used 
to  study  the  source,  intensity,  and  propauation 
chdrdctt'i  ietiCB  ol  gravity  waves  in  the  high  latitude  winter 
atmosphere,  (4)  these  measurements  prov  de  a  major  step  l h 
denionsti  ution  ol  the  capabilities  ol  lidar  an  a  tool  lor 
routine  meteorological  and  atmospheric,  measur  emen  Is. 
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ALTITUDE  <KM> 


Tdbie  1.  The  principle  components  ol  the  mobile  lidar. 

Transmitter  -  Nd:YA6  laser,  10  Hz,  U-switch,  lb  watt 
fc->00  mj  fe*  532  nm,  250  m3  (•*  355  nm 
5X  beam  expander  -  divergence  0. lb  mrad 

Receiver  -  32  cm  telescope  (steerable)  photon  counting 
with  shutter  closed  below  15  km  and  variable 
aperture  up  to  40  km  -  nominal  lov  1.0  mrad 
-  62  cm  telescope  (vertical  pointed)  photon 
counting  with  mechanical  shutter  below  40  km 

Data  System  -  Lecroy  3500  multi-channel  scaler  with  2 

microsecond  range  bins  (300  meter  resolution) 

-  data  stored  on  magnetic  disk 

SaJety  System  -  Electrical  interlocks  and  radar  used  to 

automatically  disable  the  laser 

SJMMAKY  UK  RESULTS 

On  each  evening  ol  lidar  operation,  a  series  ol  3-b  proiiles 
were  obtained  using  the  ultraviolet  wave  length  at  Jbb  nm  and 
the  visible  wavelength  at  532  nm,  simultaneously.  Figure  3 
shows  an  example  ol  the  proiiles  measured  by  the  green  and 
ultraviolet  detectors  on  the  32  cm  telescope.  Above  30  km,  the 
proiiles  are  in  complete  agreement  and  below  30  km,  the  signal 
irom  the  green  detector  is  significantly  higher  than  that  lor 
the  ultraviolet  detector.  The  ultraviolet  return  is  relatively 
larger  lor  the  molecular  scattering  and  thus  less  sensitive  to 
the  particulate  scattering.  This  fact  is  due  to  the  l/\^ 
dependence  ol  the  molecular  scattering  cross  section.  By  using 
the  two  colors,  green  and  ultraviolet,  it  is  possible  to 
resolve  the  guestion  of  what  altitude  regions  can  be  properly 
analyzed  as  pure  Rayleigh  scatter  from  which  density  and 
temperature  can  be  directly  determined. 

98 
88 
78 
88 
58 
48 
38 
28 
18 
8 

Figure  3.  Comparison  ol  the 
visible  and  ultraviolet  lidar 
profiles. 
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Figure  4.  Comparison  ol  the 
seasonal  extremes  measured 
on  14  February  and  27  April . 


DENSITY  RATIO 


F  3  i]  ui  e  r)  show:.  an  oxampJ  u  o.l  the--  t  ■>:  t  j  pini-j!:  ol  the  meat  .ur  i  *ei 

COndj  tiuas;.  1  tit  ■  pi  ol  l  1  t:*r :  iii  e  sitiown  at:  a  do  ri; ;  l  L  v  i  at  n  >  t  t  <  l.  )■  t  ■ 

UiVnAV  b  mode*!  .  I  he  profile  ol  Id  Febj  um  y  j.  typical  iiI  the 

lowe*r  winter  dent  :i  1.  y  i  ri  t.  it  e  mesosphere  an  <1  tin*  hight.-t  dent.-j  t  y 
Of  till'  summer  j  s  l  f(n  i.'Ki-ntoil  by  thi’  y./  Apj  l  L  pro!  i  Jr.  I  he  high 

Lilt,  it  u  ill*  cj  t  in  1 1  : ;  f  j  h  i 1 1  e  i  •  x  h  r  Li  i  t  s  these  t  d  j  i  jt  ■  seasonal  i:hd(uji.‘i.  1 1  ill’ 

to  the  change  from  i  u  imp  J  e  l  e  il.n  kdi'i.t:  i  n  mnlt'i  to  cun  t  j  nuous 

sun  in  the  ijiiinrinT.  The  solai  contruJ  dr  lvt'i;  Jiu  ije  changes  in 

the  cir culatiori  ul  the  uttni  Kiphcj  t-  in  the  pnJdi  r  eg  ion.  Nutt' 

that  the  i  egion  between  1 H  and  .'L>  kin  exhibits  intense  nai  i  i>w 

layei  s.  Thin  signal  it;  due  to  the*  ae:  utiiil  polar 

stratospheric.  cJuuds,  ami  vnle.an.ii-.  duwt. 

The  measurements  conducted  dui  uiij  tluu  campaniii  represent  the 

first  nerioufi  attempt  to  compel!  e  .Ildar  data  and  profiles  li  . 

metporuioiiical  rockots.  An  example  o.f  the  iidcu  prof i i e 

obtained  during  the*  1  light  ol  a  dutasoude  is;  i.huwfi  in  Figure 
The  size  of  the  i  la  error  bar  is  shown  on  the  lidar  data  with 

a  spacing  which  indicates  the  length  o.l  i  fie*  sniuuttijriii  inteiv.il. 
The  example  exhibits  rather  strong  wave*  activity,  which  wm 

observed  Ireguently  i  ri  March.  Case-!.  erl  bot.n  excel  lent 

agreement  ami  signif  icant  dil  u*i  ence  be* tween  the  reicket  and 
lidar-  results  were  obtained.  Figure  t.  shows  an  example  ol  the 

planetai  y  Wave*  variations  measured  with  the  lidar  during  a 

major  stratospheric  warming  which  peaked  on  Fin  February.  The 

mean  night  profiles  provide  a  qraphic.  picture  of  the* 
development  ol  the  str atospher ie:  warming. 


V.V 

OENSITY  RATIO 


Figure  b.  Comparison  ol  the  lidar  profile  from  the  S2.  cm 
telescope  detector  (G)  with  a  standard  datascmdi*  profiled)). 
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RLTITUDE  (KM) 


Figure  6.  Sequence  ol  mean  night  profiles  ol  the  density  ratio 
to  the  US'SA  model  which  show  the  change  in  ieiiponee  to  a  major 
stratospheric  warming. 
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GRAVITY  WAVE  DENSITY  VARIATIONS  INFERRED  FROM  GROUND  BASED 
RADAR  MEASUREMENTS  OF  WINDS  IN  THE  LOWER  THERMOSPHERE. 


Robert  A.  Vincent 
Physics  Department, 
University  of  Adelaide, 
Adelaide  5001, 

Australia 


1.  INTRODUCTION 

There  are  now  a  growing  number  of  measurements  of  winds  in  the  80 
to  100  km  region  of  the  upper  atmosphere  made  with  ground  based 
radars  operating  at  frequencies  near  2  to  3  MHz  (Vincent,  1984a). 
They  have  primarily  been  used  to  study  variations  of  mean  winds  and 
tides  in  this  region  but  recently  climatologies  of  gravity  wave 
motions  have  also  started  to  appear  (e.g.  Meek  et  al.,  1985; 
Vincent  and  Fritts,  1987).  Through  these  and  other  studies  it  has 
become  possible  to  understand  better  the  fundamental  role  that 
gravity  waves  play  in  determining  the  mean  circulation  and  thermal 
state  of  the  upper  atmosphere.  While  such  radars  have  some 
advantages  over  other  ground  based  systems,  such  as  lidars,  for 
continuous  monitoring  of  the  lower  thermosphere,  it  is  not  always 
easy  to  convert  the  wave  motions  to  other  important  wave  parameters 
such  as  density.  Nevertheless,  as  discussed  below,  it  is  possible 
with  suitable  care  to  develop  a  statistical  description  of  the 
density  perturbations  and  this  is  carried  out  for  radar 
observations  made  at  a  mid-latitude  site  in  Australia. 


2.  DATA  ANALYSIS 

The  wind  measurements  were  made  with  a  spaced  antenna  radar  at 
Adelaide,  Australia  (35°S,  138°E).  Since  the  observations  were  made 
on  an  almost  continuous  basis  over  a  3  year  period  starting  in 
November  1983,  with  a  time  resolution  of  about  10  minutes,  it  has 
been  possible  to  study  the  temporal  variations  in  wave  activity 
over  almost  the  full  wave  spectrum.  Because  the  gravity  wave 
spectrum  covers  a  wide  range  of  periods  extending  from  a  few 
minutes  out  to  the  inertial  period,  which  at  Adelaide  is  21  hr,  it 
is  necessary  to  break  the  wave  spectrum  up  into  a  number  of  bands. 
For  convenience,  and  for  ease  of  comparison  with  other  studies,  the 
data  were  filtered  and  analyzed  in  three  spectral  bands,  (i)  a  low 
or  quasi-inertial  band  covering  the  8  to  24  hr  period  range,  (ii)  a 
medium-period  band  ranging  from  1  to  8  hr,  and  (iii)  a  short-period 
band  ranging  from  20  min  to  1  hr.  In  order  to  avoid  contamination 
from  the  strong  12  and  24  hr  tidal  oscillations  the  analysis  was 
carried  out  in  10  day  blocks  and  the  tidal  winds  calculated  and 
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removed  by  harmonic  analysis.  The  data  were  then  processed  to 
produce  climatologies  of  wind  and  density  perturbations. 

.J  .  ^Uli  -  JIiUU  X  LtO 

(a)  Winds 

Figure  1  shows  a  contour  plot  of  the  mean  square  amplitude,  V'  2 , 
due  to  gravity  waves  in  the  height  range  80-98  km,  where  V'^  is 

2  ~~2  ~~2 
V'  =  (u'  +  v'z) 

and  u'  and  v'  are  the  zonal  and  meridional  perturbation  velocities 
respectively.  The  wave  amplitudes  grow  slowly,  such  that  near  86  km 
the  rms  amplitude  is  about  30  ms-^,  a  value  which  agrees  well  with 
a  previous  estimate  based  on  a  smaller  amount  of  wind  data  taken  in 
the  years  1978-1981  at  Adelaide  (Vincent,  1984b).  At  96  km  the  rms 
amplitude  is  about  40  ms-*.  Typically  these  values  may  vary  by 

between  plus  or  minus  15-20  ms-^ .  At  the  lower  heights,  there  is  a 
pronounced  semiannual  variation  in  wave  activity,  with  minima  at 
the  equinoxes  whereas  at  upper  levels  there  is  a  weaker  semiannual 
variation  with  minima  occurring  near  the  solstices. 

Comparisons  of  u'  and  v'  shows  that  v'  is  larger  than  u'  and  that 
the  covariance  term,  u'v',  is  generally  non-zero.  These 
characteristics  show  that  the  wave  field  is  anisotropic,  such  that 
a  significant  fraction  of  the  wave  motions  have  a  preferred 
direction  of  orientation  (Vincent  and  Fritts,  1987).  Using  the 
'Stokes  Parameters'  approach  discussed  in  their  paper  it  is  found 
at  the  lower  frequencies,  that  up  to  50%  of  the  wave  motions  tend 
to  lie  in  the  NW/SE  quadrants  in  summer  and  to  have  a  preference 
for  NS  directions  in  winter. 


(b)  Density  Variations 

The  conversion  of  the  wind  variations  to  other  wave  quantities, 
such  as  density,  can  in  principal  be  easily  made  by  using  the  so- 
called  polarization  relations  which  relate  one  observed  quantity  to 
another  (Hines,  1960,  Gossard  and  Hooke,  1975).  The  fractional 
density  perturbations  are  given  in  terms  of  the' wind  perturbations 
by 

S p]  2  .  (u2  -  f2)  N2  ^2  (1) 

uTT?)  ? 

Here  u>  is  the  intrinsic  frequency,  f  is  the  inertial  frequency,  N 
is  the  Vai  sala-Brunt  frequency  and  g  is  the  acceleration  due  to 
gravity.  Unfortunately,  the  use  of  this  expression  is  not 
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straightforward  since  intrinsic  frequency  (i.e.  the  frequency 
relative  to  the  background  wind)  must  be  used,  rather  than  the 
frequency  observed  in  the  ground  based  frame  of  reference. 

One  way  to  estimate  w,  is  to  make  use  of  the  polarization  relation 
between  v'  and  u' ,  which  for  a  wave  propagating  horizontally  at  an 
azimuth  9,  measured  clockwise  from  east,  is 


v'  = 


( tanG  - 


i  a) 
tanQ ) 


where  a=f/w.  This  means  that  gravity  waves  are  in  general 
elliptically  polarized  such  that  the  ratio  of  the  amplitude  of  the 
wind  component  along  the  direction  of  motion  to  the  amplitude  of 
the  transverse  component  is  w/f .  A  statistical  analysis  based  on 
the  Stokes  Parameters  approach  shows  that  rather  f ortuituosly ,  most 
of  the  waves  tend  to  be  propagating  horizontally  at  right  angles  to 
the  prevailing  winds,  so  that  at  this  location  at  least,  the 
intrinsic  and  ground  based  frequencies  are  very  similar. 

Another  subtlety  that  must  be  borne  in  mind  before  (1)  can  be 
applied  is  that  a  temperature  model  has  to  be  chosen  before  the 
Vai sala-Brunt  frequencies  can  be  calculated.  In  the  present  case, 
the  recent  model  for  the  southern  hemisphere  given  by  Groves  (1987) 
was  used  to  calculate  N. 

A  contour  plot  of  the  inferred  rms  percentage  density  variations 
are  shown  in  Figure  2.  The  percentages  vary  from  values  near  5%  at 
heights  near  80  km  upto  values  near  8-10%  at  heights  near  100  km. 
This  suggests  that  at  mid-latitudes  gravity  waves  cause  the 
atmospheric  density  to  fluctuate  by  up  to  at  least  5-10%  around  its 
mean  value  in  the  lower  thermosphere.  Taking  into  account  the 
natural  variability  of  the  wind  fluctuations  suggests  that  these 
mean  fractional  densities  will  be  uncertain  by  about  2-3%.  A 
further  cause  of  uncertainty  is  the  possibility  of  systematic 
errors  in  the  temperature  model  used  to  derive  N. 


4.  DISCUSSION 


The  present  analysis  shows  some  of  the  strengths  and  weaknesses  of 
using  radar  wind  measurements  to  infer  wave  parameters  such  as 
density.  The  excellent  time  and  height  resolution  of  the  radars 
plus  their  ability  to  run  continuously  means  that  they  can  be  used 
to  derive  climatologies  of  wave  activity.  However,  the  conversion 
of  wind  data  to  relative  density  requires  information  about  the 
intrinsic  frequencies  of  the  waves,  information  which  is  not  easy 
to  obtain  with  a  ground  based  system.  A  temperature  and  density 
model  is  also  required  to  produce  absolute  density  variations. 
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FIGURE  CAPTIONS 

Fig.  1  Contours  of  mean  square  gravity  wave  amplitudes  for  the 
height  range  80-98  km. 

Fig. 2  Contours  of  percentage  gravity  wave  density  variations. 
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4.4  INCOHERENT  SCATTER  DENSITY 


W.  L.  Oliver 

M.I.T.  Haystack  Observatory 


7G  1.  Two  primary  neutral  density  measurements  are  made 
with  the  incoherent  scatter  radar  technique:  N„  from  about 
80  to  120  km  and  0  from  about  200  to  500  km  altitude.  The 
N^  density  is  determined  directly  from  the  measured  ion- 
neutral  collision  frequency;  the  O  density  is  determined 
through  the  joint  consideration  of  ionospheric  measurements 
and  a  theory  of  ionosphere-atmosphere  interaction.  Both  of 
these  measurements  are  limited  to  daytime. 

VG  2.  Several  of  the  basic  parameters  measured  by 
incoherent  scatter  radar:  ion  and  electron  temperatures,  ion 
composition,  ion-neutral  collision  frequency. 

VG  3.  The  incoherent  scatter  spectrum,  showing  the  effects 
of  increasingly  strong  ion-neutral  collisional  effects. 

These  collisional  effects  die  out  above  about  120  km 
altitude  as  the  neutral  density  decreases.  Below  about  90 
km  the  spectrum  is  so  strongly  dominated  by  collisions  that 
it  becomes  Gaussian  in  shape;  then  it  becomes  impossible  to 
distinguish  between  temperature  and  collision  effects  and 
only  a  collision-frequency/temperature  ratio  is 
determinable.  The  degree  of  collisional  effect  is  dependent 
upon  the  ratio  of  the  radar  operating  frequency  and  the  ion- 
neutral  collision  frequency,  such  that  the  region  of  such 
density  determination  will  be  translated  several  kilometers 
upward  or  downward  according  to  the  radar  operating 
frequency. 

VG  4.  Example  of  the  neutral  density  profile  measured  at 
Arecibo,  assuming  the  CIRA  1972  temperature  profile  and  ion 
masses  as  indicated.  Negative  ion  effects  contaminate  the 
results  at  the  lower  altitudes. 

VG  5.  The  collection  of  ion-neutral  collision  frequency 
results  for  measurements  covering  a  full  day  at  Arecibo. 

Much  of  the  data  spread  results  from  real  geophysical 
variability. 

VG  6.  Example  of  daytime  measurement  of  neutral  temperature 
and  ion-neutral  collision  frequency  throughout  a  day  at 
Millstone  Hill. 

VG  7.  Daily  average  neutral  temperature  and  collision 
frequency  profiles  for  three  different  days  as  measured  at 
Millstone  Hill. 

VG  8.  The  first  incoherent  scatter  density  measurements, 
from  20  years  ago. 
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VG  9.  Daily  averages  of  neutral  density  made  at  different 
altitudes  above  Saint  Santin  during  experiments  spanning  8 
years . 

VG  10.  Mean  density  profiles  for  the  4  seasons,  computed 
from  the  data  in  VG  9 . 

VG  11.  Seasonal  amplitudes  and  phases  of  the  density, 
computed  from  the  data  in  VG  9. 

VG  12.  Solar  dependence  of  the  density,  computed  from  the 
data  in  VG  9 . 

VG  13.  Comparison  of  the  seasonal  variation  of  neutral 
density  at  100  km  altitude  as  measured  at  Saint  Santin  with 
that  of  the  Jacchia  1971  model. 

VG  14.  Comparison  of  the  radar  density  (measured  at  100  km 
and  extrapolated  to  400  km  altitude  using  the  radar-measured 
temperature  profile)  with  sevexal  model  densities. 

VG  15.  Upper  atmospheric  temperature  profiles,  illustrating 
the  method  by  which  the  neutral  temperature  and  density  are 
deduced  from  the  ionospheric  temperatures.  The  T  profile 
represents  the  neutral  temperature  computed  from  ?he  T  and 
T  profiles  (measured)  if  we  choose  the  correct  neutraf 
density  profile  for  use  in  the  ion-neutral  heat  transfer 
equations.  The  Tn'  and  T  ' '  profiles  represent  the  neutral 
temperatures  thatone  would  compute  if  he  were  to  choose  a 
neutral  density  exactly  twice  or  half,  respectively,  of  the 
correct  density.  As  we  know  from  theory  that  the  neutral 
temperature  profile  should  be  asymptotic  at  high  altitudes, 
we  then  have  a  means  of  discerning  both  the  neutral 
temperature  and  neutral  density  in  the  upper  thermosphere  by 
adjusting  our  neutral  density  assumption  until  the  deduced 
neutral  temperature  profile  becomes  asymptotic. 

VG  16.  Comparison  of  radar-measured  densities  at  400  km 
altitude  (dots)  with  a  simultaneous  satellite  measurement 
(x),  and  two  models,  showing  the  excellent  agreement  during 
the  middle  of  the  day  but  discrepancies  normally  seen  at 
other  times  of  day. 

VG  17.  Seasonal  variation  of  midday  density  at  400  km 
altitude  as  measured  at  Saint  Santin  compared  with  two 
models . 

VG  18.  The  local-time/latitude  behavior  of  exospheric 
temperature  and  neutral  density  at  400  km  altitude  as 
measured  from  Millstone  Hill  during  the  7-day  ETS  campaign 
of  17-24  September  1984. 
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V  G  1 


Two  Primary  Density  Measurements 


80  -  120  km  altitude 
0  200  -  500  km  altitude 


ion-neutral  collision  frequency 

ionosphere  measurements 
plus  theory  of 
ionosphere/atmosphere 
interaction 


Temperature  Height  Profile  Measurement 
Extrapolates  density  measurements  to  other  altitudes 


Daytime  Measurements  Onl^ 
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The  power  spectrum  of  the  scattered  signal  for  no  magnetic  field  and  various  values 
(indicated  on  the  curves'  of  the  normalized  ion  collision  frequency  T(.  \l>,  is  taken  to  be 
<r  t/10  Oi  is  the  normalized  Doppler  shift  d'T:ni  1  by  (i,  —  (w/k)  (.m</2KT)' 3.  The  full  curve 
is  svmmctric  about  0 1  =  0. 


ALTITUDE  (km) 
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NEUTRAL  NUMBER  DENSITY  (xlOl9m‘3) 
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VG  5 


80  85  90  95 


ALTITUDE  (km) 

Collision  frequency  (»/,„),  approximate  \p,  values,  and  in¬ 
ferred  neutral  number  densities  are  plotted  on  the  ordinate  (logarith¬ 
mic  scale)  versus  altitude  on  the  abscissa  (linear  scale).  The  solid  linen 
the  best  linear  approximation  to  the  96  data  points  and  can  be  com¬ 
pared  with  the  model  atmosphere  shown  as  the  dashed  line.  The 
extremes  of  temperature  in  the  model  atmosphere  over  the  altitude 
range  shown  is  from  1 95°K  at  80  km  to  I  83°K  at  90  km  and  to  193#K 
at  97  km  with  a  mean  of  I90°K  over  the  entire  range.  The  observations 
indicate  an  atmosphere  which  is  both  hotter  and  denser  than  ik 
model  atmosphere. 
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COLLISION  FREQUENCY  (sec' 


Daily  oscillations  of  tho  middatitudo  thormosphoro 


Collision  frequency  measurements  on  14  April  1971  at  Millstone  Hill 
converted  to  neutral  number  density. 


NEUTRAL  NUMBER  DENSITY  (cm'3) 


ALTITUDE  <Kf1>  ALTITUDE  ( KH ) 


'115 
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N2  CONCENTRATIONS  (01-3) 


M2  CONCENTRAT I ONS  <0*1-3  > 


The  seasonal  fit  of  observed  molecular  nitrogen  concentrations  leads  to 
characteristic  values  which  are  shown  for  F^q  ^  =  125  as  a  function  of 
altitude.  The  error  bars  are  computed  from  the  uncertainties  on  the 
coefficients  of  equation  (2).  The  full  lines  correspond  to  mixing 
distributions  with  the  temperature  profiles  of  Fig.  7  and  the  dashed 
lines  represent  Jacchia’s  (1977)  values. 
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VG  17 
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VG  18 
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Introduction 


The  incoherent-scatter  technique  is  unique  in  its  broad  capabilities  for  upper- 
atmospheric  and  ionospheric  research.  From  the  electromagnetic  radiation  backscattered 
from  the  ionosphere,  an  extensive  set  of  parameters  descriptive  of  the  physical,  dynamical, 
and  chemical  state  of  the  upper  atmosphere  and  ionosphere  can  be  determined.  This  tech 
nique  allows  an  extended  altitude  region  to  be  probed  simultaneously.  Depending  on  the 
observing  mode,  measurements  over  a  wide  range  of  latitudes  and  longitudes  can  be  made 
continuously  for  extended  intervals  (up  to  several  days)  with  time  resolution  ranging  from 
seconds  to  minutes. 

The  return  signal  spectrum  and  auto-correlation  function  (ACF)  obtained  from  the 
E  region  allow  several  quantities  to  be  determined,  including  the  electron  density,  the 
electron  and  ion  temperatures,  the  ion  velocity,  and  the  ion-neutral  collision  frequency. 
Figure  1  shows  an  example  of  some  of  these  measured  quantities  on  a  day  of  considerable 
geomagnetic  activity  in  the  auroral  zone  at  Chatanika,  Alaska  (65° N).  These  quantities, 
in  conjunction  with  atmospheric  and  ionospheric  models,  can  be  used  to  determine  other 
parameters  important  at  E- region  heights.  These  include  ionospheric  conductivities,  en¬ 
ergy  deposition  by  auroral  particles,  Joule  heating,  electrojet  currents,  and  neutral  wind 
velocity,  temperature,  and  density. 

Of  particular  interest  here  is  the  calculation  of  neutral  density  and  winds  in  the  lower 
thermosphere.  Both  of  these  calculations  depend  critically  on  the  ion-neutral  collision 
frequency,  i/tn,  which  is  non-negligible  at  these  heights.  The  ion-neutral  collision  frequency 
can  be  measured  directly  from  the  changes  in  shape  of  the  ACF  that  result  from  collisions 
between  ions  and  neutrals.  Alternatively,  it  can  be  estimated  by  using  atmospheric  and 
ionospheric  models.  Determination  of  the  neutral  winds  at  a  given  altitude  additionally 
requires  knowledge  of  the  ion  velocity  at  that  height  and  the  electric  field  vector  imposed 
from  above,  both  of  which  axe  measured  directly  by  the  radar. 
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In  the  next  section,  the  calculation  of  Uin  and  the  total  neutral  number  density  U 
from  incoherent  scatter  radar  data  will  be  discussed.  The  following  section  addresses  the 
calculation  of  the  neutral  wind  velocity  u  and  its  dependence  on  i/,n.  The  importance  of 
lower  thermospheric  tides  and  geomagnetic  activity  effects  at  high  latitude  is  also  discussed. 
Finally,  the  potential  for  improvements  on  the  determination  of  >/  with  this  technique  is 
discussed. 


Ion-neutral  collision  frequency  and  neutral  density 

Determination  of  i /tn  from  incoherent-scatter  measurements  proceeds  from  an  analysis 
of  the  scattered  signal  ACF.  As  collisions  of  ions  with  neutrals  become  more  frequent  with 
decreasing  altitude,  the  characteristic  two-humped  form  of  the  collisionless  incoherent- 
scatter  spectrum  changes  into  a  Lorentzian  profile  [Dougherty  and  Farley,  1963).  Over  the 
altitude  range  in  which  this  shape  transition  occurs,  corresponding  roughly  to  the  lower 
thermosphere,  Tj,re,  and  Vin  can  theoretically  all  be  calculated  by  a  least-squares  fitting 
procedure,  although  usually  only  two  parameters  can  be  fitted  simultaneously  with  suc¬ 
cess  [Lathuillere  et  al.,  1983],  Below  those  altitudes,  the  ACF  shape  at  upper- mesospheric 
heights  is  a  function  of  the  ratio  of  the  to  t/«n;  additional  assumptions  are  therefore 
required  to  obtain  Vin  from  the  measurement  [Tepley  et  al.,  1981;  F1&  et  al.,  1985].  Ad¬ 
ditionally,  the  ACF  broadens  to  the  point  that  i /,-«  cannot  be  determined  efficiently  from 
a  single-pulse  or  multipulse  measurement.  Instead,  a  pulse-to-pulse  technique  must  be 
employed  [Mathews,  1976;  Mathews,  1978;  Tepley  and  Mathews,  1978;  Tepley  et  al.,  1981; 
Kofman  et  al.,  1984],  As  a  result,  the  method  of  data  acquisition  and  analysis,  as  well 
as  the  requisite  assumptions  differ  in  the  upper  mesosphere  and  the  lower  thermosphere. 
Our  interest  here  is  in  the  lower  thermosphere,  where  ACFs  formed  from  single-pulse  or 
multipulse  experiments  can  be  analyzed  to  obtain  uin.  This  technique  has  been  used  by 
Wand  and  Perkins  [l968j,  Salah  [1974],  Salah  et  al.  [1975],  Wand  [1976],  Schlegel  et  al. 
[1980],  Wickwar  et  al.  [1981],  Lathuillere  et  al.  [1983],  Kofman  and  Lathuillere  [1985], 
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F1A  et  al.  [1985],  Kofman  et  al.  [1986],  and  fluuskonen  [1986]  at  Arecibo,  Millstone  Hill, 
Chatanika,  and  EISCAT. 

One  difficulty  encountered  in  determining  from  incoherent-scatter  measurements 
arises  in  the  theoretical  specification  of  collision  process  in  incoherent-scatter  theory 
[Dougherty  and  Farley,  1963;  Tanenbaum,  1968;  Waldteufel,  1970a, b;  Hagfors  and  Brock- 
elman,  1971].  Depending  on  the  treatment  of  collisions  in  the  derivation  of  the  incoherent- 
scatter  spectrum,  slight  differences  result  in  the  dependence  of  the  ACF  on  uin.  Hagfors 
and  Brockelman  [1971]  suggested  that,  as  a  result,  U{n  obtained  using  the  Dougherty  and 
Fariey  jl963j  theory  would  be  slightly  underestimated,  and  T,  slightly  overestimated.  How¬ 
ever,  in  a  comparison  of  the  two  theories,  Salah  et  al.  [1975]  found  that  in  the  107  to  119 
km  height  range  at  Millstone  Hill,  the  collision  frequencies  determined  with  the  Dougherty 
and  Farley  model  exceeded  those  of  Hagfors  and  Brockelman  at  the  lowest  levels  by  up  to 
30%,  attaining  better  agreement  as  altitude  increased.  Ganguly  and  Coco  [1987]  compared 
the  theoretical  results  of  Dougherty  and  Farley  [1963]  and  Waldteufel  [1970a, b]  with  one 
another  as  well  as  with  experimental  data  in  the  collision  dominated  regime,  and  gave 
strong  support  to  the  formulation  of  Dougherty  and  Farley.  There  has  been  no  satisfac¬ 
tory  resolution  of  these  problems  to  date,  theoretically  or  experimentally.  Except  for  the 
studies  by  Salah  et  al.  [1975]  and  Ganguly  and  Coco  [1987],  most  work  in  this  area  has 
proceeded  from  the  theory  of  Dougherty  and  Farley. 

Another  difficulty  arises  in  determining  the  relationship  of  uxn  to  the  neutral  den¬ 
sity.  The  form  of  this  dependence  is  itself  dependent  on  the  dominant  mode  of  collisional 
interaction  between  ions  and  neutrals.  It  is  generally  assumed  at  these  heights  and  tem¬ 
peratures  that  the  induced  dipole  interaction  is  the  most  important  collisional  mechanism 
[Banks,  1966].  However,  resonant  charge  exchange  has  also  been  suggested  as  a  competing 
process. 

In  the  height  range  of  interest  here  (90-150  km),  the  dominant  ionic  species  are  N0+ 
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and  0^"  (of  comparable  concentrations),  while  N2  and  O2  are  the  dominant  neutral  species. 
As  altitude  increases,  the  concentration  of  atomic  oxygen  ( [O])  increases  until  it  exceeds  [02j 
above  120  km,  while  0+  does  not  become  an  important  constituent  relative  to  0^  and  N0+ 
until  about  150  km.  Resonant  charge  exchange  between  0^  and  O2  becomes  important 
only  at  temperatures  above  «500  K,  while  for  0+  and  0  that  transition  temperature  is 
«200  K  [Banks,  1966).  Thus  resonant  charge  exchange  can  compete  with  induced  dipole 
interactions  at  the  upper  limits  of  the  altitude  range  considered  here,  where  temperatures 
approach  the  500  K  level  for  0^" —  O2  collisions.  However,  the  importance  of  resonant  charge 
exchange  for  0+ —  0  collisions  is  reduced  because  both  0+  and  0  are  minor  constituents. 
Thus,  for  normal  E-region  conditions,  where  temperatures  range  from  200  to  400  K  as 
altitude  increases  from  90  to  about  125  km,  and  where  [0+]-c[N0+]  or  [0^],  resonant 
charge  exchange  can  be  neglected. 

Under  such  conditions,  the  model  of  polarization  potentials  is  adequate  to  describe 
the  collisional  process.  As  was  shown  by  Lathuillere  et  al.  [1983],  the  total  ion-neutral 
collision  frequency  at  E-region  heights  can  be  expressed  as 

E  2  59  x  10“ 11 

ujn  ^  — - (5.07%jV2  +  4.98%02  +  2.88%0)  U,  (l) 

uU.J 

j  =  NO*tO* 

where  utn  is  the  collision  frequency  of  the  jth  ion  with  neutral  species  n,  M  is  the  total 
neutral  number  density,  and  %n  represents  the  relative  proportion  of  the  neutral  species 
n  to  the  neutral  gases  as  a  whole.  In  the  derivation  of  this  formula,  a  single  ion  with  a 
mass  intermediate  between  N0+  and  0^"  is  assumed,  and  collisions  of  0+  with  the  neutral 
constituents  are  neglected.  This  expression  shows  that  if  a  model  of  neutral  composition  is 
available,  U  can  be  determined  from  the  ion-neutral  collision  frequency  determined  with 
the  incoherent-scatter  technique. 

It  is  difficult  to  assess  the  variability  of  the  neutral  composition  with  time,  since  few 
measurements  are  available  in  the  literature.  The  results  of  Offerman  et  al.  (1981j  indicate 
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variations  in  [N2]  and  [O2]  of  about  5%,  while  atomic  oxygen  has  been  observed  to  vary 
by  factors  of  4  to  5.  However,  over  the  entire  height  range  in  question,  N2  dominates 
over  atomic  oxygen  ranging  from  a  factor  of  100  at  90  km  to  a  factor  of  2  at  150  km. 
Thus,  the  effect  of  potentially  variable  atomic  oxygen  concentration  on  deduction  of  the 
neutral  density  is  reduced  because  the  contributions  from  N2  and  O2  collisions  dominate 
the  calculation  over  most  of  the  altitude  range. 

As  was  mentioned  briefly  above,  usually  the  data  quality  do  not  allow  Tt,Te,  and 
1 /{n  to  be  simultaneously  determined  from  the  ACF.  Most  studies  to  date  have  therefore 
assumed  thermal  equilibrium  between  ions  and  electrons,  fitting  instead  for  T,  =  Te  and 
Uin  alone.  While  the  assumption  of  thermal  equilibrium  is  valid  at  low  to  mid  latitudes, 
several  physical  processes  which  occur  in  the  auroral  zones  can  cause  selective  heating  of  the 
electrons  or  the  ions  [I  athuillere  et  al,  1983].  As  a  result,  periods  of  high  electric  fields  must 
be  excluded  from  determination  of  z/,n  when  working  with  high-latitude  measurements. 

Figure  2  shows  an  example  of  i/m  calculated  from  Chatanika  radar  data  obtained  in 
March  and  November  of  1978,  in  comparison  with  values  calculated  using  the  Jacchia-71 
model  [Lathuillere  et  al.,  1983].  This  example  of  good  agreement  between  experimental 
i/in  profiles  and  empirical  model  profiles  is  typical  of  the  results  obtained  over  the  upper- 
mesosphere  and  lower-thermosphere  region  [Wand  and  Perkins,  1968;  Salah,  1974;  Wand, 
1976;  Tepley  and  Mathews,  1978;  Schlegel  et  al.,  1980;  Tepley  et  al.,  1981;  Kofman  and 
Lathuillere,  1985;  F14  et  al.,  1985;  Kofman  et  al.,  1986;  Huuskonen  et  al.,  1986;  Gan¬ 
guly  and  Coco,  1987].  These  studies  provide  empirical  justification  for  the  approach  of 
Dougherty  and  Farley  (1963),  as  well  as  the  assumption  that  induced-dipole  interactions 
are  of  dominant  importance  in  this  height  range. 
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E-region  neutral  winds 


The  ion  velocity  along  the  line-of-sight  of  the  radar  is  obtained  from  ‘he  Doppler 
shift  of  the  return  signal  spectrum.  Through  the  use  of  different  pulse  lengths  for  E- 
and  F-region  measurements,  altitude  resolution  appropriate  to  both  of  these  regions  can 
be  obtained.  Combination  of  three  successive  line-of-sight  measurements  with  different 
pointing  directions  allows  vector  ion  velocities  with  time  resolution  of  the  order  of  ten 
minutes  or  less  to  be  obtained.  The  E-region  ion  velocity,  along  with  the  electric  field, 
enter  into  the  calculation  of  the  neutral  wind. 

The  height-dependent  neutral  wind  velocity,  u(z),  is  deduced  from  the  steady  state 
ion  equation  of  motion  appropriate  to  E-region  heights.  Perpendicular  to  B,  this  equation 
reduces  to 

«x(*)  =  v±{z)  ~  (~^  (£l  +  x  &)  »  (2) 

where  E±  is  the  ionospheric  electric  field  derived  from  F-region  ion  velocity  measurements, 
v±(z)  is  the  height-dependent  E-region  ion  velocity,  and  fl,  is  the  ion  gyrofrequency.  Par¬ 
allel  to  B,  the  ion  equation  of  motion  reduces  to 

U|((«)  =  t>||(*)  -Vambi(z),  (3) 

where  vomt»t  is  the  ambipolar  diffusion  velocity.  At  E-region  heights,  vambt{z)  is  usually 
negligible,  as  are  vertical  velocities  relative  to  the  horizontal  neutral  wind  components. 

Results  from  calculations  of  the  horizontal  neutral  wind  usually  show  the  presence  of 
large  amplitude  low-frequency  motions  [Rino  et  al.,  1977;  Wand,  1983a, b;  Johnson  et  al., 
1987].  The  semidiurnal  fluctuation  is  particularly  important  at  high  latitudes.  Figure  3 
summarizes  the  tidal  results  for  summer  months  at  Chatanika  [Johnson  et  al,  1987].  The 
four  different  polar  plots  show  the  variation  with  local  time  of  the  diurnal  (see  panel  a)  and 
semidiurnal  (panel  b)  tidal  components  obtained  from  a  least-squares  fit  to  the  average  of 
the  data  at  115  km  altitude.  Also  shown  are  the  resultant  net  tidal  wind  field  (panel  c), 
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and  the  total  average  flow  field  in  which  the  mean  flow  at  115  km  has  been  added  (panel 
d).  Wind  velocities  of  75  to  100  m/s  are  common  at  this  altitude,  and  the  semidiurnal 
component  is  clearly  dominant. 

Equation  (2)  shows  the  dependence  of  the  neutral  wind  on  the  local  ion  velocity  and 
electric  field.  During  intervals  of  geomagnetic  activity,  enhanced  electric  fields  drive  the 
ions  into  motion  and  produce  variations  in  the  neutral  winds  as  well.  The  study  by  Johnson 
et  al.  [1987]  showed  that  neutral  winds  responded  in  a  consistent  and  systematic  fashion 
above  sslOO  km  at  Chatanika  during  summer  months.  This  response  was  particularly 
strong  in  the  local  morning  sector,  where  zonal  flow  enhancements  of  several  hundred  m/s 
were  found  on  occasion.  Figure  4  shows  the  zonal  geographic  neutral  wind  at  115  km 
obtained  for  June  11,  1980  (open  circles)  as  well  as  the  average  zonal  wind  for  Kp  <  2 
conditions  (thin  line)  as  a  function  of  universal  time.  At  Chatanika,  local  dawn  occurs  at 
16  UT.  Throughout  the  local  morning  sector,  the  zonal  wind  was  strongly  enhanced  over 
the  Kp  <  2  average  values,  up  to  «300  m/s  at  20  UT. 

Equation  (2)  also  shows  that  the  neutral  wind  velocity  is  dependent  on  t/in.  The 
smaller  i/,n  becomes,  the  more  the  ion  velocity  diverges  from  the  neutral  wind.  The  inverse 
form  of  this  dependence  causes  uncertainties  in  and  v  to  result  in  proportionately 

larger  uncertainties  in  ujl(«)  at  greater  heights.  Figure  5  shows  how  uncertainty  in  i/tn  can 
affect  the  resultant  neutral  wind  obtained  through  Equation  (2).  Three  different  neutral 
density  profiles  were  calculated,  corresponding  to  75%,  100%,  and  125%  of  the  MSIS-83 
value.  The  panel  at  the  left  compares  the  profiles  of  the  zonal  neutral  wind  determined  us¬ 
ing  these  density  profiles  at  a  time  when  the  second  term  on  the  right-hand  side  of  equation 
(2)  dominates  the  calculation.  On  the  right-hand  side  of  the  figure,  the  zonal  ion  velocity 
dominates  the  calculation.  Thus,  depending  on  the  relative  magnitudes  of  the  ion  velocity 
components  and  the  electric  field,  the  uncertainty  produced  as  a  result  of  uncertainty  in 
i/,n  can  range  from  negligible  to  the  order  to  the  neutral  velocity  itself  (assuming  a  25% 
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uncertainty).  Obviously,  in  order  to  determine  the  neutral  wind  accurately,  the  ion-neutral 
collision  frequency  must  be  known  as  accurately  as  possible. 


Potential  for  Improvements 

At  Chatanika  [Schlegel  et  al.,  1980;  Lathuillere  et  al.,  1983]  and  at  EISCAT  [F1A  et  ah, 
1985;  Kofman  and  Lathuillere,  1985;  Kofman  et  ah,  1986;  Huuskonen  et  ah,  1986],  the  ion- 
neutral  collision  frequency  could  not  be  determined  reliably  from  fits  to  the  experimental 
data  at  altitudes  above  about  110  km.  Above  this  height,  the  influence  of  collisions  on 
the  autocorrelation  function  (ACF)  of  the  return  signal  was  too  small  to  significantly 
distort  the  shape  of  the  ACF  at  the  radar  operating  frequencies  (1290  MHz  and  933  MHz, 
respectively) . 

The  parameter  which  characterizes  this  distortion  is  the  normalized  ion  collision  fre¬ 
quency,  =  (^nA/47r)(m,/2«;T,)1/2,  where  A  is  the  radar  operating  wavelength,  and 
Ti  are  the  ion  mass  and  temperature,  and  k  is  Boltzmann’s  constant.  As  (or  i/tn) 
decreases,  the  ACF  of  the  return  signal  gradually  develops  the  normal  oscillatory  shape 
characteristic  of  the  collisionless  case  and  the  damping  influence  of  U{n  can  eventually  no 
longer  be  detected. 

Figure  6  [after  Dougherty  and  Farley,  1963]  shows  curves  of  corresponding  to 
operating  frequencies  of  50,  430,  and  1290  MHz  versus  altitude.  The  vertical  line  at 
ss  0.04,  which  intersects  the  1290  MHz  line  at  110  km,  indicates  the  minimum  value  of 
'J',  from  which  t/*n  could  be  determined  reliably  from  the  experimental  ACFs  at  Chatanika. 
At  lower  frequencies,  this  same  value  of  is  attained  at  higher  altitudes  («120  km  at 
430  MHz  and  wl40  km  at  50  MHz).  This  indicates  that  with  the  same  system  sensitivity 
available  at  Chatanika,  the  influence  of  i/in  on  the  ACFs  could  be  detected  to  higher 
altitudes  at  lower  frequencies.  Indeed,  Salah  [1974]  and  Wand  and  Perkins  [1968]  report 
measurements  of  i/{n  up  to  120  km  at  Millstone  Hill  (440  MHz)  and  at  Arecibo  (430  MHz), 


4-50 


respectively.  With  increased  system  sensitivity,  obtained  through  increased  transmitter 
power  and  the  use  of  low-noise  receivers  these  upper  altitude  limits  could  probably  be 
raised  by  about  half  an  atmospheric  scale  height  (several  km). 


Whether  these  techniques  can  be  used  to  deduce  uxn  (and  thereby  neutral  densities)  at 
altitudes  above  125  km,  as  is  suggested  by  the  50  MHz  curve  on  the  figure,  depends  on  the 
influence  of  other  factors.  As  was  noted  by  Dougherty  and  Farley  [1963],  other  scattering 
processes  may  affect  the  measurement  as  the  operating  frequency  decreases.  Another  factor 
is  degradation  in  signal-to-noise  ratio  at  low  frequencies  because  of  the  increase  in  noise 
temperature  from  the  sky.  Finally,  the  interpretation  of  uxn  as  a  measure  of  neutral  density 
needs  to  be  reexamined  at  greater  heights  due  to  the  increased  importance  of  0+  as  well  as 
resonant  charge-exchange  interactions  and  the  increased  importance  of  potentially  variable 
atomic  oxygen  density.  A  study  of  these  factors  and  their  effect  on  the  measurement  and 
interpretation  of  win  is  required  in  order  to  determine  the  potential  of  this  technique  at 
these  greater  altitudes. 
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Figure  Captions 


Fig.  1  —  Contour  plots  of  E  region  electron  density  (cm-3,  upper  left),  ion  temperature 
(K,  upper  right),  electron  temperature  (K,  lower  left),  and  ion  velocity  (m/s,  lower 
right)  versus  altitude  and  time  on  June  11,  1980.  Minimum  and  maximum  values  are 
shown  in  the  lower  left  and  right  comers  of  each  panel,  respectively. 

Fig.  2  —  Profiles  of  mean  collision  frequency  for  March  and  November,  1978  Chatanika 
radar  measurements  (solid  lines  with  error  bare)  compared  with  Jacchia-71  model 
values  (dashed  lines)  (after  Lathuillere  et  al.  (1983]). 

Fig.  3  —  Diurnal  (a)  and  semidiurnal  (b)  wind  vectors  determined  from  a  least-squares  fit 
to  the  average  of  the  summer  Chatanika  data  at  115  km  in  a  polar  plot  format.  The 
sum  of  these  two  tidal  components  is  shown  in  panel  (c),  and  panel  (d)  includes  the 
mean  flow  as  well  [Johnson  et  al.,  1987j. 

Fig.  4  —  Zonal  neutral  wind  at  115  km  in  geomagnetic  coordinates  on  June  11,  1980  (open 
circles)  compared  with  the  average  for  Kp<  2  conditions  (dot-dashed  curve).  Local 
midnight  at  Chatanika  occurs  at  10  UT  (Johnson  et  al.,  1987], 

Fig.  5  —  Altitude  profiles  of  the  zonal  neutral  wind  on  June  11,  1980  at  two  different 
universal  times.  The  three  curves  on  each  panel  show  the  neutral  wind  obtained 
using  the  MSIS-83  density  profiles  multiplied  by  0.75  (dashed),  1.0  (solid),  and  1.25 
(dotted). 

Fig.  6  —  Normalized  ion  collision  frequency,  =  (i/mA/47r)(m,/2KT1,) as  a  function 
of  altitude  at  sunspot  maximum.  Curves  for  three  different  operating  frequencies  are 
shown.  The  vertical  line  indicates  the  minimum  value  of  from  which  i/tn  could  be 


deduced  with  the  Chatanika  radar  (1290  MHz)  (after  Dougherty  and  Farley,  [1963)). 
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1.0  INTRODUCTION 

Orbital  reentry  of  large  maneuverable  vehicles  and  the  potential  to  modify 
orbital  parameters  using  the  aerobreaking  technique  have  increased  the  need  to 
understand  in  detail  the  variability  of  atmospheric  density  in  the  mesosphere  and 
lower  thermosphere.  Acceleration  data  1  obtained  during  shuttle  reentry,  shown  in  I 
figure  1.,  illustrate  the  fact  that  the  atmosphere  exhibits  variation  which  often 
exceed  expectation  and  which  raise  questions  regarding  the  safe  return  of  heavy 
Shuttle  payloads  from  polar  orbit. 

The  variability  of  mass  density  in  the  mesosphere  has  been  apparent  for  many 
years,  being  seen  in  sounding  rocket  data,  backs,  utter  radar  profiles,  and  airglow 
optical  emissions.  The  systematic  variations  tha'  have  been  determined  from  these 
data  sets  are  the  obvious  seasonal,  latitudinal,  and  local  time  dependences  that  are 
prescribed  in  the  major  models  in  use  today.  These  models,  however,  are  based  on 
limited  data  and  thus,  are  not  useful  in  making  predictions  or  forecasts  of  the  effects 
of  planetary  waves  or  gravity  wave  breaking.  The  wave  processes  will  never  be 
able  to  be  accurately  predicted,  but  like  the  weather  in  the  troposphere  we  may 
ultimately  be  able  to  base  a  forecast  on  operational  information.  A  primary  factor 
holding  back  the  improvement  of  these  models  and  the  development  of  a  forecasting 
technique  is  the  lack  of  a  sufficiently  large  and  accurate  morphological  data  base. 

This  global  data  base  must  rely  on  a  sampling  technique  that  will  provide  a 
morphology  of  the  atmospheric  density  in  the  critical  60  to  120  km.  region,  with 
the  consistency  required  to  quantify  the  higher  frequency  variability.  The  only 
practical  method  of  providing  global  information  is  to  use  a  remote  sensing 
technique  from  earth  orbit,  given  that  it  has  the  required  accuracy  and  spatial 
sensitivity.  Such  a  technique  could  be  tested  using  first  sub-orbital  evaluation, 
followed  by  shuttle  validation,  and  finally  deployed  on  a  satellite  in  earth  orbit  to 
provide  the  required  global  information. 
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Several  possible  remote  sensing  techniques  have  been  proposed  which  have  the 
potential  to  measure  the  atmospheric  density  in  this  important  atmospheric  region. 
These  include  measuring  the  brightness  of  sunlight  that  is  rayleigh  scattered  by  the 
molecular  nitrogen  and  oxygen  2,  measuring  the  absorption  of  starlight  during  stellar 
occultation  3,  and  a  new  technique  where  the  fluorescence  of  molecular  oxygen  is 
measured  photometrically.  This  new  method  has  evolved  from  the  studies  of  the 
University  of  Michigan's  Space  Physics  Research  Laboratory  which  is  providing  the 
High  Resolution  Doppler  Imager  for  the  Upper  Atmosphere  Research  Satellite  4.  In 
order  to  evaluate  the  potential  for  a  spacebome  remote  measurement  of  density  in  the 
60  to  120km.  region  it  is  necessary  to  compare  the  relative  merits  of  these  related 
techniques.  Each  of  these  methods  of  sensing  the  density  has  inherent  advantages  as 
well  as  weaknesses. 

2.0  TECHNICAL  DESCRIPTION  AND  APPROACH 

The  University  of  Michigan  Space  Physics  Research  Laboratory  is  in  the 
process  of  evaluating  several  means  of  obtaining  the  density  morphology  in  the 
Mesosphere  and  lower  Thermosphere  by  remote  sensing  techniques.  At  this  time 
three  candidate  methods  have  been  identified,  these  being: 

1.  Resonance  fluorescence  of  molecular  oxygen 

2.  Absorption  of  starlight  during  stellar  occultation 

3.  Rayleigh  scattering  of  sunlight . 

The  order  that  these  have  been  listed  is  based  on  the  expectation  of  their 
relative  usefulness  in  accomplishing  the  objective,  however,  the  present  study  may 
on  completion  lead  to  a  re-ordering  of  these  methods,  or  addition  to  uic  list  uf 
potential  techniques. 

2.1  Resonance  fluorescence  of  molecular  oxygen 

The  daytime  Atmospheric  "A"  Bands  of  molecular  oxygen  provide  one  of  the 
brightest  spectral  emission  features  that  are  seen  in  the  visible  and  near  infra-red 
region  of  the  spectrum  of  light  emitted  from  the  earths  atmosphere.  Yet  this  band 
system  has  recieved  very  little  attention  since  the  early  investigations  of  Wallace 
and  Hunten  5,  and  the  more  recent  observations  obtained  from  the  Dynamics 
Explorer  Satellite  by  Skinner  and  Hays 6 .  This  band  system  results  from 
competition  between  processes  which  excite  molecular  oxygen  to  the  b'£g+  state 
and  those  that  depopulate  the  state,  where: 

A)  Excitation  processes 

1.  Resonance  absorption  (7650  A) 
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02(X3!)  +  hv  ->  02(b!E) 

2.  Production  of  O(^D)  and  quenching  (  Ultra-violet) 

02(X3!)  +  hv  ->  O('D)  +  0(3P)  N 

Production 

03  +  hv-4  0(>D)  t-  02  J 

0(!D)  +  0:(X3X)  ->  02(b'Z)  +  0(3P)  }  Quenching 

3.  Chemical  recombination 

0(3P)  +  0(3P)  +  M  02*  +  M 
02*  +  0(3P)  -4  02(b‘Z)  +  0(3P) 

B)  Loss  processes 

1.  Radiation 

02(b:I)  — >  02(X3X)  +  hv 

2.  Quenching 

02(b‘I)  +  M— >  02(X3X)  +  heat 


The  various  sources  of  production  of  molecular  oxygen  in  the  singlet  sigma 
state  are  shown  in  figure  2,  where  resonance  excitation  and  photodissociation  of  02 
are  combined.  When  losses  are  considered  with  these  production  processes  in  a 
simple  photochemical  equilibrium  calculation,  we  are  able  to  determine  the 
contribution  each  makes  to  the  volume  emission  rate  as  a  function  of  altitude. 

Figure  3  shows  the  volume  emission  rate  due  to  each  of  these  sources  of  excited 
molecular  oxygen.  Note  that  the  dominant  sources  of  excited  molecular  oxygen  in 
the  60  to  140  km.  range  are  all  directly  proportional  to  the  ground  state  oxygen 
density,  and  thus  are  directly  related  to  the  total  density  in  that  region  of  the 
atmosphere  where  the  mixing  ratio  is  well  known.  Below  60  km.  the  dissociation  of 
ozone  must  be  taken  into  account.  We  note  that  in  this  region  the  source  of  emission 
can  be  obtained  from  an  examination  of  the  intensity  of  the  (1-0  )  Atmospheric  band 
at  8650A  which  is  optically  thin. 


4-67 


The  limb  brightness  observed  in  the  (0-0)  "A"  band  at  7650 A  is  determined 
by  emission  in  the  various  rotational  lines  of  the  band,  attenuated  through  the 
molecular  oxygen  between  the  observer  and  the  point  of  the  emission.  Two  typical 
band  distributions  are  shown  in  figure  4  ,  where  the  distribution  of  emission  is  seen 
for  200  and  1000  0  K  .  Since  these  lines  are  self  absorbed,  the  column  density  of  the 
absorber  (here  molecular  oxygen)  can  be  determined  from  low  resolution  band 
spectra  of  the  light  scattered,  not  only  in  the  optically  thin  region  at  high  altitude, 
but  in  the  optically  thick  region  below.  This  is  achieved  by  an  inversion  process, 
illustrated  symbolically  in  figure  5,  where  resonance  absorption  is  indicated  to  be 
the  principal  source  of  excited  molecular  oxygen.  Detailed  spectra  have  been 
calculated  as  a  part  of  the  High  Resolution  Doppler  Imager  studies,  conducted 
under  the  Upper  Atmosphere  Research  Satellite  program  at  Michigan  since  the  late 
1970's.  Typical  rotational  line  brightness  for  various  line  strengths  is  shown  in 
Figure  6,  where  a  line  with  strength  of  20  cm’/atmos.km.  has  a  brightness  of  nearly 
300  kilorayleighs  when  viewed  at  a  tangent  ray  height  of  80  km.  This  system  is 
extremely  bright  and  offers  one  of  the  best  possible  sources  of  information  about 
the  density  structure  in  the  upper  mesosphere  and  lower  thermosphere. 

2.2  Stellar  Occultation  and  Absorption 

Stellar  occultation  is  used  to  determine  the  density  of  absorbing  species  in  the 
atmosphere  by  a  technique  that  is  similar  to  classical  absorption  spectroscopy.  The 
Sun,  Moon,  or  a  star  are  used  as  a  source  of  light  and  the  Satellite  photometer  as  a 
detector,  with  the  atmosphere  between  acting  as  the  absorption  cell.  As  the  satellite 
moves  in  its  orbit,  the  source  untimately  is  occulted  by  the  Earth.  During  the 
occultation  process,  the  intervening  atmosphere  absorbs  progressively  more  of  the 
light  due  to  the  strong  photo-absorption  features  in  the  atmospheric  gases.  The 
occultation  intensity  data  can  be  related  to  the  tangential  column  number  density  of 
the  absorbing  species  if  the  absorption  cross-sections  are  known.  This  technique  has 
been  investigated  in  detail,  with  the  early  theoretical  studies  of  Hays  and  Roble  3 
describing  the  technique,  errors  were  examined  by  Roble  and  Hays  7  ,  and  the 
recent  work  has  been  reviewed  by  Atreya  8 .  Figure  7  illustrates  the  geometry  and 
methodology  used  in  stellar  occultation. 

Molecular  oxygen  is  the  obvious  candidate  for  observation,  the  ultra-violet 
absorption  cross-section  for  this  species,  shown  in  figure  8,  has  features  that  allow 
the  density  to  be  observed  throughout  the  region  of  interest.  Typical  absoption 
profiles  are  shown  in  figure  9  taken  from  reference  10  for  a  few  selected  spectral 
regions.  We  note  that  observations  near  Lyman  alpha  line  of  hydrogen  allow 
density  to  be  observed  down  to  about  90  km.,  this  is  illustrated  in  the  observations 
reported  by  Atreya  8  and  shown  in  figure  10.  Lower  altitudes  would  be  observable 
using  bands  in  the  Schumann  Runge  Bands  near  1900  A.  This  region  was  observed 
on  the  OAO-A2  spacecraft ,  but  were  not  analyzed  due  to  the  complexity  of  the 
band  system.  Since  that  time  detailed  studies  of  these  bands  have  been  carried  out 
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and  it  would  be  reasonable  to  expect  that  densities  could  be  extracted  in  the  critical 
60  to  100  km.  region. 


2.3  Ravleish  Scattering  of  Sunlight 

The  major  gases  of  the  atmosphere  scatter  light  due  tc  the  processes  of 
Rayleigh  scattering  in  which  non-resonance  polarization  of  the  molecules 
results  in  scattering  of  incomming  light.  During  the  daytime  period  the  blue 
sky  results  from  the  preferential  scattering  at  shorter  wavelengths.  Thus, 
like  the  resonance  scattering  technique  Rayleigh  scattering  can  be  used  to 
determine  the  density  of  the  atmosphere  at  high  altitudes.  This  method  of 
observation  is  illustrated  in  figure  1 1  where  the  observation  and  inversion 
processes  are  indicated  schematically.  This  technique  has  been  used  on  the 
Solar  Mesosphere  Explorer  Satellite  to  determine  the  density  and 
temperature  in  the  upper  stratosphere  2.  Questions  regarding  this  method  of 
measuring  the  density  arise  due  to  the  existance  of  aerosols  and  Polar 
Mesospheric  Clouds9.  These  particles  and  clouds  cause  scattering  that  is 
difficult  to  distinguish  from  density  variations  in  the  basic  atmosphere.  This 
problem  may  be  resolved  using  a  multiplicity  of  wavelengths. 
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ABSTRACT 


Measurement  of  emission  from  the  earth's  atmosphere  in  the 

far  infrared  due  to  transitions  between  fine  structure  levels  of 

3  3  3 

the  ground  state  (  P)  of  atomic  oxygen  at  63pm  (  P^  — >  P2)  and 

147pm  ( 3Pq  — >  3p^)  is  proposed.  These  magnetic-dipole  allowed 

trcinsitions  with  long  radiative  lifetimes  ( -  3.2  hrs  for  the  63pm 

transition  and  -  16.3  hrs  for  the  147pm  transition)  are  assumed  to 

be  in  equilibrium  with  the  local  translational  temperature.  A 

one-dimensional  onion-peel  inversion  of  the  limb  emissions  at 

63pm  and  147pm  from  a  model  atmosphere  is  shown  to  yield 

reasonable  results  for  both  the  temperature  and  atomic  oxygen 

density  in  the  90-250  km  altitude  range. 


! 
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I.  INTRODUCTION 


Atomic  oxygen  is  an  important  constituent  of  earth's  atmos- 

phere  above  90  km  and  the  most  abundant  constituent  above  about 

200  km  [U.S.  Standard  Atmosphere,  1976).  It  has  long  been  known 

to  play  an  important  role  in  the  heat  balance  and  chemistry  of 

the  upper  atmosphere.  Bates  [1951]  suggested  that  the  fine  struc- 

3  3 

ture  transition  0(  P^  )  — >  0(  P 2 )  at  63/m  plays  an  important  part 
in  cooling  the  upper  thermosphere.  The  5.3pm  emission  from  vibra- 
tionally  excited  nitric  oxide,  which  in  the  quiescent  atmosphere 
is  produced  solely  [Caledonia  and  Kennealy  (1982)]  by  collisions 
of  nitric  oxide  in  the  ground  state  with  atomic  oxygen,  was  shown 
by  Kockarts  [1980],  Gordiets  et  al.  [1982],  Stair  et  al.  [1985] 
and  Zachor  et  al.  [1985]  to  be  the  dominant  cooling  mechanism  in 
the  lower  thermosphere.  Around  100  km,  15pm  emission  from  the 
bending  mode  of  C02<  which  again  is  excited  primarily  by  the 
collisions  of  oxygen  atoms  with  ground  state  carbon  dioxide 
[Sharma  et  al.  (1981)],  plays  a  dominant  role  in  cooling  the 
atmosphere  [Dickinson  (1984)).  At  still  lower  altitudes  [Sharma 
(1985)],  atomic  oxygen  plays  an  important  role  in  the  emission 
for  the  hydroxyl  airglow  layer  and  chemiluminescence  from  03  via 
the  reactions 


H  +  03  — >  OH*  +  02 
H02  +  0  — >  OH*  +  02 
O  +  02  +  M  — >  03*  +  M 
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W t  thus  see  that  atomic  oxygen,  directly  or  indirectly,  plays  a 
vital  role  in  the  heat  balance  of  the  upper  atmosphere. 

Offermann,  et  al.  (1981]  have  reviewed  the  methods  used,  all 
ir.  situ,  for  the  measurement  of  atomic  oxygen  densities  in  the 
upper  atmosphere.  Sharp  11985)  has  commented  on  results  of  these 
measurements:  "The  data  show  a  variation  in  the  concentrations  by 
over  a  factor  of  10  at  120  km,  a  factor  of  40  for  the  peak  at 
about  95  km  and  much  larger  factors  below  this.  Peak  concentra¬ 
tions  range  from  6  x  10^  cm~^  to  2  x  10^  cm  Within  any  one 
technique  there  are  factors  of  5  variations  over  altitude.  Some 
of  this  variation  can  be  due  to  time  of  day,  season  and  dynamics. 
Hc-vever,  the  principal  source  of  the  variation  is  probably  the 
differing  measurement  techniques.” 

While  several  methods  for  measuring  the  density  of  atomic 
ox/gen  in  the  thermosphere  have  been  undertaken  and  critically 
re/iewed,  we  have  not  come  across  corresponding  experiments  or 
critical  reviews  for  the  measurement  of  thermospheric  tempera¬ 
ture.  Zachor  et  al.  (1985)  tried  obtaining  temperature  as  a 
fuiction  of  altitude  from  the  envelope  of  the  limb  emission 
observed  from  vibration-rotation  transitions  of  nitric  oxide  at 
5.3  pm  observed  in  the  110-150  km  tangent  height  range.  The 
results  obtained  were  not  satisfactory.  This  may  be  due,  in  large 
pact,  to  the  rather  weak  dependence  of  the  vibration-rotation 
baid  envelope  on  temperature.  Although  the  problem  may  perhaps  be 
re nedied  by  measuring  emission  from  carefully  selected  lines  of 
thi  vibration-rotation  band,  it  has  not  been  tried  thus  far.  In 
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sedition,  the  nitric  oxide  emission  decreases  rapidly  [Stair  et 
al .  (1985)]  above  150  km  due  to  decreasing  density,  and,  there¬ 
fore,  above  this  tangent  height  other  methods  for  determining  the 
temperature  have  to  be  developed. 

In  this  work  we  describe  a  method  of  measuring  the  oxygen 
atom  density  and  temperature  by  remote  sensing .  The  method  con¬ 
sists  of  simultaneous  measurement  of  intensity  in  the  63  pm 
3P1  — >3P2  transition  and  the  147pm  3pq— >3PX  transition  (Figure 
1)  between  the  ground  state  fine  structure  levels.  We  assume  that 
these  fine  structure  levels  are  in  thermodynamic  equilibrium  with 
the  local  translational  temperature  (LTE),  and  that  the  atmos¬ 
phere  is  one-dimensional.  Radiance  in  the  63pm  and  147  pm  lines 
is  calculated  using  a  model  atmosphere  and  a  modified  form  of  the 
radiation  transport  code  of  Sharma  et  al.  [1983J.  An  onion-peel 
inversion  is  then  carried  out  to  show  that  it  is  possible  to 
recover  the  density  and  temperature.  Although  the  present  inver¬ 
sion  was  performed  on  the  "data"  which  had  no  noise,  work  in 
progress  adds  noise  to  the  radiance  "data"  before  carrying  out 
tho  inversion.  The  principal  aim  of  this  work  is  to  show  that  it 
is  possible  to  develop  reliable  procedures  for  the  recovery  of 
ox/gen  atom  density  and  temperatures  as  a  function  of  altitude 
from  the  limb  radiance  data. 

The  assumption  that  the  fine  structure  levels  are  in  ther¬ 
modynamic  equilibrium  with  the  translational  temperatures  (LTE) 
neods  some  elaboration.  Grossmann  and  Offermann  [1978]  concluded 
from  their  measurements  at  63pm  that  this  assumption  is  not  a 
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valid  one  and  that  the  population  of  the  ^P^  level  is  less  than 
thcit  given  by  the  LTE  model.  Iwagami  and  ugawa  (1980]  argue  that 
the  Grossmann-Of fermann  results  are  consistent  with  LTE  formula¬ 
tion  if  radiative  transfer  is  properly  treated.  Our  computation 
of  the  63  radiance  using  the  model  atmosphere  kindly  supplied 
by  Dr.  D.  Offerman  agrees  with  Grossmann-Of fe rmann  results. 
Grossman  et  al  (1985)  have  re-examined  the  situation  and  have 
suggested  that  if  a  value  of  the  Einstein  Coefficient 
a12“4 . 5x10“5s_1  for  the  63 pm  transition  is  used  (nearly  one-half 
the  calculated  value),  the  computed  radiances  are  in  reasonable 
agreement  for  one  (1981)  of  the  two  sets  of  rocket  experiments. 

We  are  not  quite  sure  of  uncertainties  involved  in  the  model 
atmosphere.  For  example,  an  increase  of  about  25%  in  atomic 
oxygen  density  would  bring  the  results  obtained  from  this  model 
atmosphere  in  agreement  with  the  emission  observed  by  the  1981 
rocket  experiment.  There  are  serious  reasons  not  to  either  alter 
the  value  of  the  radiative  lifetime  of  the  excited  state  or 

to  discard  the  LTE  hypothesis.  In  an  elaborate  study,  Fischer  and 
Saha  (1983]  have  calculated  the  radiative  lifetime  of  the  63pm 
transition  to  be  3.192  hrs  (11,491  seconds)  (Fig.  1).  Allison  and 
Burke  [1969]  estimate  the  cross-section  for  transition  between 
the  fine  structure  levels  during  collisions  between  two  oxygen 
atoms  to  be  between  0.01  and  1  times  the  gas  kinetic  cross- 
section.  At  an  altitude  of  250  km  such  collisions  occur  about 
or.ce  every  second  or  about  11,000  collisions  per  radiative 
lifetime,  enough  to  populate  the  P^  state  in  accordance  with 
LTE. 
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Lin  et  al.  (1987)  have  conducted  a  baloon  experiment  which 
measures  63pm  emission  at  elevation  angles  of  +30°  and  +1°.  From 
this  experiment  these  authors  assuming  LTE  and  using  calculated 
value  of  Einstein  Coefficient  A^»  have  derived  peak  atomic 
oxygen  density  of  3.6  (+  1.9)x  10**  cm**^  at  an  altitude  of  about 
110  km.  We  would  like  to  reiterate  that  we  are  not  foreclosing 
the  possibility  that  either  the  LTE  hypothesis  or  the  values  of 
the  Einstein  Coefficients  may  have  to  be  modified  in  light  of 
further  observations.  However,  the  best  evidence  available  at 
this  time  favors  these  hypotheses. 

it  has  been  pointed  out  earlier  that  the  method  proposed 
here  is  a  passive  remote  sensing  technique.  The  advantages  of  the 
method  are  best  exploited  by  a  space-shuttle  or  a  satellite 
experiment  which,  in  principle,  can  yield  not  only  diurnal, 
seasonal,  latitudinal  and  variations  due  to  geomagnetic  distur¬ 
bances  but  also  time  constants  for  reaching  new  steady  states. 
There  is  no  reason  why  the  method  outlined  here  should  not  be 
used  to  study  other  planets.  For  example,  in  both  Venus  and  Mars, 
the  atomic  oxygen  is  believed  to  be  an  important  constituent  of 
the  atmosphere  above  100  km  (Chamberlain  (1978));  the  method 
outlined  here  may  be  used  to  obtain  quantitative  results. 

II.  THEORY  AND  METHODOLOGY 

A  derivation  of  the  radiative  transfer  equation  has  been 
given  by  Sharma  et  al.  (1983],  Sharma  and  Zachor  (1983)  and 
Zachor  and  Sharma  (1985).  We  present  here  a  brief  summary  of  the 
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formalism.  In  one-dimensional  models  of  the  atmosphere,  the 
atmospheric  parameters,  e.g.,  temperature,  density,  etc.,  depend 
orly  on  the  altitude  z.  To  calculate  the  emission  and  transmis¬ 
sion  of  radiation  we  divide  the  line-of-sight  (LOS)  into  segments 
of  thickness  (Az)  such  that  the  atmospheric  parameters  over  each 
segment  may  be  assumed  to  be  constant. 

We  express  the  radiance  in  terms  of  n^,  the  number  of 
photons  (of  frequency  v  era-*)  per  wavenumber  per  unit  area  per 
steradian  which  are  observed  at  a  point  z  along  the  line-of-sight 
per  unit  time.  The  change  in  n^  in  a  path  length  ds  is  described 
by 

dnyds  -  {-n^(  hvQ/c)  [  n-j^Bt  l->u)-nuB(u->l )  ]+Aulnu/4  ji)  f  (  v-vQ  )  (1) 

whare  h  and  c  are  Planck's  constant  and  the  speed  of  light;  n^ 
an i  nu  are  the  number  densities  of  the  lower  and  upper  radiating 
states;  B(l->u)  and  B(u->1)  are  Einstein  coefficients  for  absorp- 
ti  an  and  induced  emission;  is  the  Einstein  coefficient  for 

spontaneous  emission,  and  f  is  a  normalized  lineshape  function 
su  :h  that 


j  f ( v-\q ) dv  -  1 

The  Einstein  coefficients  are  related  by 
g1B(l->u)  -  guB(u->l) 


(2) 


(3) 


4-88 


and 


Au]>  -  8  nhVgB  ( u->l )  (4) 

where  g^  and  gu  are  the  statistical  weights  of  the  lower  and 
upper  radiating  states.  Through  use  of  Eqs.  (3)  and  (4),  Eq.  (1) 
can  be  cast  in  the  form 

dnv/ds  -  (hv0/c)B(l->u)n1(-nv(l-r)  +  2cv^ y ) f ( v-vQ )  (5) 

where  y  is  defined  as 

Y  "  SlWl  (6) 

The  factor  (1-y)  corrects  the  simple  absorption  term  (which  is 
proportional  to  n^n^)  for  stimulated  emission.  The  second  term  in 
brackets  corresponds  to  spontaneous  emission. 

We  now  define  the  optical  depth,  tv,  by 

dr^/ds  -  (hv0/c)B(l->u)f(v-v0)n1(l-y)  (7) 

Suostituting  Eq.  (7)  in  Eq.  (5)  we  get  the  result 

dnv/dxv  -  [-nv  +  2cVqy/(1-y))  (8) 

This  equation  has  a  very  simple  physical  meaning.  The  rate  of 
change  of  the  number  of  photons  with  absorption  depth  decreases 


4-89 


as  the  numbers  of  photons  incident  on  the  segment  and  increases 
with  the  source  function.  Under  conditions  of  LTE  the  source  term 
is  seen  equal  to  the  black-body  function  when  we  recall  that  y 
( Eq .  5)  under  these  conditions  equals 

Y  -  expl -(  Eu-E1  )/kT]  -  exp  [ -C2n>q/T  ]  (9) 

where  Eu  and  E^  are  the  energy  levels  of  the  upper  and  lower 
states;  k  is  the  Boltzmann  constant;  C2  [■  1.4308  (K/cm-^)]  is 
the  second  radiation  constant,  and  T  is  the  temperature. 

This  equation  can  easily  be  integrated  across  a  segment  of 
constant  temperature  and  density.  The  solution  is 

n\.(zi  +  l)  "  nv<  )exp(-Axvl )  +  2cvQ(l-exp(-ATvi  )  ]  Yi/{  l-ri  )  (10) 

where  the  segment-dependence  of  r  has  been  made  explicit  with  the 
subscript  i.  In  Eq .  (10),  the  first  terra  represents  the  absorption 
within  the  i th  segment  of  radiance  incident  upon  it  at  z^,  and 
the  second  term  represents  the  contribution  from  the  ith  segment 
alone;  Ar^  the  change  in  optical  depth  across  a  segment  can  be 
written  as 


ATv  “  (Aulgu/8nc\.Qg1)n1(l-Y)f (v-v0)As 
where  As  is  the  length  of  the  segment. 


(11) 
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In  this  study  the  Einstein  coefficients  for  ground-state 
atomic  oxygen  along  with  assumed  density,  temperature,  and  pres¬ 
sure  atmospheric  profiles  were  input  data  for  a  modified  form  of 
the  NLTE  computer  code  developed  at  the  Air  Force  Geophysics 
Laboratory  (Sharma  et  al .  [1983),  Wintersteiner  and  Sharma 
[1985]).  The  modifications  were  accomplished  such  that  the 
program  calculates  the  radiance  for  the  0  P2  <-  0  transition 
at  63/>m  and  the  <-  O^Pq  transition  at  147pm  in  the  earthlimb 

(or  Zenith)  geometry  for  any  range  of  tangent  heights  (or  al¬ 
titudes)  within  the  range  of  the  assumed  atmospheric  profile. 

The  atmospheric  density  and  temperature  profiles  used  in 
this  study  are  plotted  in  Figure  2.  This  density  profile  was 
selected  by  comparison  with  previous  estimates  of  O-atom  den¬ 
sities,  and  represents  a  reasonable  guess  at  the  true  density. 

Two  additional  density  profiles  corresponding  to  atmospheres  two 
times  and  five  times  as  rich  in  atomic  oxygen  as  that  given  in 
Figure  2  were  also  studied  in  order  to  test  the  accuracy  and 
reliability  of  the  numerical  inversion  technique. 

The  integrated  line  radiances  calculated  from  Eq.  2  are 

plotted  versus  tangent  height  in  Figure  3  for  both  transitions 

appropriate  for  the  atmospheric  profile  given  in  Figure  2.  In 

these  calculations  it  was  assumed  that  the  energy  levels  were  in 

local  thermodynamic  equilibria.  A  voigt  lineshape  with  a  col- 

lisional  line  width  equal  to  0.1  cm”1  per  amagat  was  assumed  in 

-1/2 

all  of  the.  work  reported  here.  A  weak  T  '  dependence  of  col- 
lisional  linewidth  upon  temperature  [Penner  (1959)]  was  ignored. 
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For  the  purpose  of  this  initial  study,  as  was  pointed  out  ear¬ 
lier,  it  is  assumed  that  this  data  corresponds  to  "exact"  total 
integrated  radiance  data  for  the  two  transitions.  In  the  next 
section,  we  describe  a  numerical  procedure  to  "invert"  this  data 
for  the  purpose  of  recovering  the  assumed  density  and  temperature 
profiles.  Consideration  of  errors  due  to  finite  signal  to  noise 
and  other  experimental  factors  will  be  presented  and  discussed  in 
a  later  paper. 

III.  Numerical  Inversion  of  the  Integrated  Radiance 

In  order  to  start  the  inversion  procedure,  it  was  assumed 
that,  at  the  highest  tangent  heights  considered,  the  transitions 
are  optically  "thin".  That  is,  no  significant  self-absorption 
takes  place.  In  this  case,  the  integrated  radiance  for  the  tran¬ 
sition  i <  —  j ,  Njj,  depends  only  on  the  excited  state  (03Pj)  den¬ 
sity  profile,  which  is  directly  related  to  the  ground  state 
density  through  the  temperature  under  the  assumed  conditions  of 
thermodynamic  equilibrium.  Formally,  we  may  write 

Nij(V  “  <l/2n)J  e.  j(z}W(z,HT)dz,  (12) 

where  ht  is  one  tangent  altitude  and  e^j(z)  is  the  total  volume 
emission  rate  at  altitude  z. 

j ( z  )  -  n*( z)Ai j  ( 13 ) 

In  Eq.  (12),  W(z,HT)  is  a  geometrical  weighting  factor  that 
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converts  distance  dl  along  limb  line-of-sight  to  altitude  dz,  Wdz 
-  dl , 

W(z,Ht)  -  { R+z ) [ ( R+z ) 2  -  (R+HT)2 J"1/2  (14) 

where  r  denotes  the  radius  of  the  earth. 

The  high  altitude  density  profile  is  obtained  by  first 
assuming  that  the  density  decreases  exponentially 

n j M z )  -  n j * ( HT) exp( -0z )  (15) 

Using  Eqs.  (13)  and  (IS)  in  Eq.  (12)  we  obtain  the  result 

Nij(HT)  -  (l/2n)n*(HT)Al;JJ  exp(-0z)W(z,HT)dz  (16) 


The  integral  in  Eq.  (16)  can  be  evaluated  analytically  with 
th-2  result 

Nlj(HT)  -  Ai  .nJ(HT)exp[0(R+HT)  ](R+HT)K10(R+HT)  ]  (17) 

is  the  modified  Bessel  function.  The  two  unknowns  in  this 
equation  [n^  and  0]  can  be  evaluated  by  considering  the  in¬ 
tegrated  radiance  at  the  two  highest  tangent  heights  for  which 
radiance  data  is  available.  An  approximate  excited  state  density 
profile  for  the  highest  altitude  regions  can  then  be  constructed 
from  Eq.  (15).  We  have  made  the  additional  assumption  that  the 
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temperature  in  this  region  varies  in  a  linear  fashion,  und  can  be 
approximated  from  known  temperature  profiles.  These  approxima¬ 
tions  appear  to  be  reasonable,  and,  in  fact,  as  long  as  the 
highest  measurement  is  far  removed  from  the  atmospheric  region  of 
interest  these  "starting'  values  have  little  effect  on  the 
results  at  lower  tangent  heights. 

Inversion  of  the  radiance  through  the  lower  tangent  heights 
if.  accomplished  by  solving  for  the  minimum  in  the  function 

UT,n,HT)  -  [N21(HTj  -  nP1(Ht)]2  +  (N;0(Ht)  -  nPq(Ht)]2  (10) 

n2i(ht)2  nio(ht)2 

nPj(Ht)  denotes  the  predicted  radiance  at  tangent  height,  .  The 
tvo  terras  in  this  equation  have  been  separately  normalized  in 
Older  for  them  to  be  weighted  equally.  In  the  wor,.  described 
here,  we  have  also  multiplied  this  function  by  a  "penalty"  func¬ 
tion  in  order  to  avoid  mathematical  diversions  involving  large 
temperature  changes.  Such  diversions  in  temperature  are  not 
ur common  in  problems  of  this  type  in  which  the  variabl  s  T  and  n 
aie  somewhat  correlated.  The  resultant  function  is  given  below 

f(T,n,HT)  -  UT,n,HT)exp(0.01|AT|  )  (19) 

where  AT  -  T  -  T  .  T  denotes  the  temperature  of  the  previous 

r  r 

slab,  and  the  bar  denotes  that  this  function  has  been  normalized. 
II  some  constraint  is  not  placed  upon  the  temperature  variable, 
then  relatively  large  oscillations  in  temperature  may  occur  in 
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certain  regions  of  function  space.  These  ten>oorary  oscillations 
around  the  true  solution  are  due  to  a  very  flat  function  in  whicn 
relatively  large  changes  in  temperature  can  be  offset  by  opposing 
all  changes  in  density. 

The  inversion  technique  that  we  have  chosen  to  use  in  this 
initial  study  is  a  so-called  "onion-peel"  method.  In  this  ap¬ 
proach,  one  determines  the  best  value  for  the  density  and  tem¬ 
perature  for  the  atmospheric  slab  of  thickness  1  km  at  the  cur¬ 
rent  tangent  height.  Following  the  generation  of  a  high  atmos¬ 
phere  profile,  as  described  above,  we  begin  the  solution  of  eq. 
(19),  by  choosing  as  our  initial  guess,  the  oxygen  atom  density 
and  temperature  of  the  previous  slab.  The  limb  radiance  at  the 
current  tangent  height  depends,  of  course,  on  the  density  and 
temperature  of  all  slabs  of  higher  altitude.  This  requires  the 
inversion  to  be  performed  with  high  accuracy  at  each  tangent 
height,  so  that  the  cumulative  errors  do  not  result  in  making 
inversion  at  lower  tangent  heights  impossible. 

In  order  to  minimize  the  number  of  computations,  for  each 
tangent  height  all  of  the  frequency  independent  parts  of  the 
radiative  transfer  equation  are  calculated  only  once.  Even  with 
this  effort,  because  the  path  length  through  each  slab  changes 
for  each  tangent  height,  the  calculated  total  radiance  for  the 
two  transitions  is  computationally  expensive,  especially  for  low 
tangent  heights  where  the  number  of  slabs  is  large.  It  is  impor¬ 
tant,  therefore,  to  select  a  numerical  procedure  that  will  lead 
to  rapid  convergence.  A  common  approach  to  rapid  solutions  of 
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problems  of  this  complexity  is  to  use  so-called  conjugate  methods 
in  which  the  good  attributes  of  more  than  one  method  are  com¬ 
bined.  For  this  reason,  we  have  chosen  to  solve  for  the  minimum 
in  Eq .  (19)  through  a  combined  "steepest  descent"  and  Newton- 
Raphson  numerical  procedure.  Efficiency  of  the  solution  of  Eq. 
(19)  is  improved  due  to  the  insensitivity  of  the  steepest  descent 
method  to  the  initial  guess,  and  the  rapid  (quadratic)  conver¬ 
gence  of  the  Newton-Raphson  method.  This  numerical  technique  is 
commonly  referred  to  as  a  "modified  Newton  method"  (Bazarra 
(19  79)  ]  . 

Formally,  we  write  the  solution  vector  for  the  (k+l)th  slab, 
xk.i  ■  (nk+i'Tk.i>  as  £ollows 

xk.l  *  xk  -*  M  +  Xdl»gtH»  I-1  VJ|X  (201 

whe re  X  is  a  parameter  selected  to  make  the  matrix  in  brackets 
negative  definite,  and  H  is  the  Hessian  matrix  of  second  deriva¬ 
tives. 

(H)ij  a  (32J/3oiOj)  {«i  -  n;  -  T)  (21) 

Additional  savings  in  computation  is  achieved  by  approximat¬ 
ing  the  Hessian  matrix  of  second  derivatives  (which  is  required 
in  the  Newton-Raphson  procedure)  by  the  vector  product  of  the 
gradient  ( vj )  and  its  transpose.  If  the  matrix  in  eq.  (20)  is  not 
invertible-,  then  the  e igenvalue  and  eigenvectors  of  the  Jacobian 
matrix  of  derivatives  are  evaluated.  This  process  essentially 
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transforms  the  dependent  set  (n,  T)  into  an  approximately  inde- 
penlent  set  of  pseudovariables  which  are  derived  from  the  eigen¬ 
vectors.  The  eigenvalues  of  the  Jacobian  are  used  for  scaling 
purposes.  This  numerical  method  is  a  variation  of  the  method 
developed  by  Marquardt  [1970].  In  the  work  described  here,  this 
conjoined  method  is  found  to  converge  for  initial  guesses  far  from 
the  true  value,  while  it  converges  at  least  supralinearly  for 
convex  regions  "close"  to  the  true  solution  [Harlow  and  Riehl  (to 
be  publ i shed ) ] . 

One  cycle  of  the  numerical  procedure  described  above  yields 
a  next  guess  for  the  solution  vector  x^.  For  each  new  estimate  of 
T^  and  n^,  boundary  conditions  are  checked,  and,  if  exceeded,  the 
result  is  reflected  back  into  a  preset  allowed  region.  The 
process  is  repeated  until  UT,n,HT)  <  1  x  IQ-6. 

IV.  Results  and  Conclusions 

The  results  from  the  inversion  procedure  described  above  are 
plotted  in  Figures  4,  5,  and  6  for  the  three  different  atmos¬ 
pheric  profiles.  The  solid  lines  in  these  figures  are  the  assumed 
temperature  and  density  profiles.  In  Figure  4,  it  can  be  seen 
that  the  inversion  was  started  at  a  tangent  height  of  200  km, 
i.e.  below  the  highest  altitude  used  in  the  calculation  of  the 
radiances  (230  km),  in  order  to  simulate  the  kind  of  data  that 
would  be  available  from  an  actual  experimental  measurement.  In 
Figures  5  and  6,  the  numerical  inversion  was  begun  at  the  highest 
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value  of  the  assumed  atmospheric  profile.  (In  these  cases  the 
profile  was  extended  out  to  300  km.)  As  can  be  seen  in  these 
figures,  the  errors  associated  with  the  assumptions  for  the  high 
altitude  regions  are  not  large,  and  are  significantly  reduced 
after  a  small  number  of  inversion  steps. 

As  is  evident  from  the  results  in  Figures  4-6,  the  inversion 
is  quite  successful  for  both  temperature  and  density  until  the 
very  "optically  thick"  regions  of  the  atmosphere  are  reached. 
Below  about  90  km  (  95  km  for  the  high  density  atmospheres)  self¬ 
absorption  is  so  strong  that  the  integrated  radiances  for  the  two 
transitions  are  almost  constant.  This  can  be  seen,  for  example, 
in  the  plots  of  radiance  versus  tangent  height  in  Figure  3, 
Continuation  of  the  inversion  through  these  regions  leads  to 
large  oscillations  in  both  temperature  and  density. 

All  three  results  show  one  or  more  small  discontinuities  in 
the  inversion  results  especially  in  the  recovered  temperature 
values.  In  Figure  4,  for  example,  these  occur  at  approximately 
120  and  130  km.  These  discontinuities,  which  are  not  unusual  for 
nonlinear  optimization  problems  of  this  complexity,  are  due  to 
the  fact  that  multiple  solution  pathways  exist  at  these  par¬ 
ticular  tangent  heights.  These  oscillations  are  significantly 
reduced  through  use  of  the  "penalty"  function  as  described  above. 

Extensions  of  the  work  reported  here  are  underway.  In  par¬ 
ticular,  more  advanced  numerical  techniques  are  under  study,  and 
procedures  to  invert  radiance  data  in  which  noise  and  finite 
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spectral  resolution  have  been  added  are  under  development.  It  is 
felt  that  measurements  of  emission  at  63pm  and  147pm  have  enor¬ 
mous  potential  in  remote  sensing  of  atomic  oxygen  and  tempera¬ 
ture. 
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Figure  Captions 


1.  Schematic  energy  level  diagram  for  the  ground  state  (^P)  of 
atomic  oxygen. 


2.  Assumed  atmospheric  temperature  and  density  profile  for  atomic 
oxygen . 

2 

3.  Calculated  radiance  (watts/cm  ster.)  versus  tangent  height 
(km)  for  the  63  pm  and  147  pm  transition  of  atomic  oxygen 
using  the  atmospheric  profile  given  in  Figure  2. 


4.  Oxygen  atom  density  and  temperature  results  obtained  from 
numerical  inversion  for  the  atmospheric  profile  given  in 
Figure  2.  The  solid  lines  are  the  assumed  density  and 
temperature  profiles. 


5.  Oxygen  atom  density  and  temperature  results  obtained  from 
numerical  inversion  for  an  atmospheric  profile  in  which  the 
atomic  oxygen  density  is  two  times  that  given  in  Figure  2. 
The  solid  lines  are  the  assumed  density  and  temperature 
profiles  . 


6.  Oxygen  atom  density  and  temperature  results  obtained  from 
numerical  inversion  for  an  atmospheric  profile  in  which  the 
atomic  oxygen  density  is  five  times  that  given  in  Figure  2. 
The  solid  lines  are  the  assumed  density  and  temperature 
profiles . 
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ATMOSPHERIC  DENSITY  MEASUREMENTS 
WITH  THE  REMOTE  ATMOSPHERIC  AND 
IONOSPHERIC  DETECTION  SYSTEM  (RAIDS) 
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INTRODUCTION 

The  scientific  objectives  of  the  Remote  Atmospheric  and 
Ionospheric  Detection  System  (RAIDS)  include  the  measurement  of  the 
dayside  and  nightside  F-region  electron  density  profile,  neutral 
composition  and  temperature,  energetic  electron  flux  distribution 
and  several  minor  species  abundances.  The  dayside  electron  density 
profile  will  be  determined  from  measurement  of  the  multiply 
scattered  0+  dayglow  emission  in  the  extreme  ultraviolet  (EUV) 
primarily  at  834  8.  The  nightside  electron  density  and  temperature 
will  be  determined  from  analysis  of  0+  and  recombination  in  the  far 
ultraviolet  (FUV)  1300-1700  8  spectral  region  and  at  7774  8  and 
from  C>2  recombination  at  6300  8.  The  neutral  composition  will  be 
determined  from  N2  and  emissions  and  02  absorption  m  the  1300-1700 
8  band,  and  neutral  0  atom  recombination  in  the  2500-3200  8  O2 
Herzberg  I  bands.  Energetic  electron  information  will  be  obtained 
from  the  N2  airglow  in  the  Lyman - B i r ge - H op f i e 1 d  (LBH)  bands  (  1  3  0  0  - 
1700  8)  and  the  N2  second  positive  and  first  negative  systems  in 
the  near  ultraviolet.  Density  distributions  of  several  minor 
species  such  as  H,  He,  N,  NO,  Mg,  Mg+,  and  Na  will  be  obtained  from 
measurements  of  their  emission  lines  and  bands.  The  sodium 
measurements  will  also  provide  verification  of  the  satellite 
attitude.  A  more  detailed  description  of  the  RAIDS  remote  sensing 
techniques  is  given  below. 

The  RAIDS  experiment  will  be  flown  on  a  N0AA  TIROS  satellite 
which  occupies  a  sun- synchronous  polar  orbit  at  830  km  at  a  local 
time  of  0730.  RAIDS  is  currently  selected  to  fly  on  the  TIROS  J 
satellite  which  is  scheduled  for  launch  in  September  1991  although 
there  is  some  possibility  that  RAIDS  will  fly  on  an  earlier  TIROS 
scheduled  to  launch  15  months  sooner.  The  actual  launch  date 
depends  on  the  health  of  one  or  more  of  the  operational  satellites 
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in  the  orbiting  TIROS  constellation. 

lse  RAIDS  instrument  complement  includes  seven  limb-scanning 
and  one  limb  imaging  sensors  viewing  in  the  orbital  plane  of  the 
satellite.  The  limb- scanning  sensors  are  mounted  to  a  mechanical 
scan  mechanism  wnich  scans  the  instrument  f i e 1 d s - o f - v i e w  through 
tangent  altitudes  75-750  km.  The  limb-imager  continuously  views 
altitudes  75-300  km  tangent  altitude.  The  scan  platform  takes 
about  97  seconds  to  complete  a  limb  scan  in  which  time  the 
satellite  has  moved  about  5  degrees  in  latitude. 


REMOTE  SENSING  OF  NEUTRAL  TEMPERATURE: 


Neutral  atmospheric  temperatures  will  primarily  be  determined 
with  the  RAIDS  experiment  by  measurement  of  topside  scale  heights 
of  emission  profiles  of  a  number  of  airglow  emitters  including  N2, 
0,  and  NO.  The  primary  candidates  for  these  temperature 
measurements  are  those  optically  thin  emissions  produced  by 
photoe lec t ron  impact  excitation.  At  the  higher  altitudes  the 
photoelectron  flux  is  essentially  independent  of  altitude  so  the 
scale  heights  depend  on  the  density  of  the  emitting  species. 
Similarly,  scale  heights  of  optically  thin  emission  produced  by 
resonant  or  fluorescent  scattering  of  sunlight  depend  only  on  the 
density  of  the  emitting  species  as  long  as  the  chemical  lifetimes 
are  long . 

Additional  techniques  for  temperature  sensing  with  RAIDS  will 
include  measurements  of  the  rotational  line  intensity  distribution 
in  molecular  bands  such  as  the  02  Atmospheric  bands.  Even  with 
limited  spectral  resolution  it  will  be  possible  to  sample  the 
intensity  distribution  over  several  parts  of  the  band.  The  results 
will  be  compared  to  synthetic  spectra  to  provide  the  temperature. 
This  technique  will  b *  most  useful  during  the  daytime  at  low 
altitudes  . 


ATOMIC  OXYGEN  DENSITY 


DAYSIDE : 

The  most  promising  technique  so  far  for  obtaining  the  0  density 
involves  measurement  of  the  ratio  of  emission  lines  at  1304  R  and 
1641  R. .  Meier  and  Conway  (Geophys.  Res.  Lett.,  12,  601,  1985)  have 

shown  that  this  ratio  is  a  linear  function  of  the  0  density  and 
should  provide  the  0  density  above  120  km.  Additional  techniques 
for  0  determination  involve  modelling  of  photoelectron  impact 
excitation  of  0  emissions  (eg.  1356  ft)  and  0+  emissions  (eg.  617 

ft)  . 

Atomic  oxygen  densities  can  be  inferred  indirectly  by 
elimination  of  the  effects  of  quenching  of  airglow  emissions.  For 


1  10 


example,  the  upper  excited  state  of  N2  Ve r ga r d - Ka p 1  an  (VK)  bands 
are  susceptible  to  quenching  and  observations  of  these  emissions 
could  yield  the  0  density  in  the  120-200  km  altitude  region. 

The  delta  bands  of  nitric  oxide  are  produced  by  the  two-body 
association  of  N  and  0  atoms.  If  the  NO  delta  bands  can  be  reliably 
observed  on  the  dayside  above  the  rayleigh  scattering  background 
these  emissions  would  yield  the  0  density  down  to  100  km  or  lower. 

NIGHTSIDE  : 

The  principal  techniques  for  measuring  0  at  night  using  the 
RAIDS  observatory  will  involve  the  measurements  of  nightglow 
emissions  produced  by  the  release  of  chemical  energy  associated 
with  0  atom  recombination.  The  primary  candidates  are  the  02 
Herzberg  I  bands  (0  +  0  +  M  ->  02  +  M),  the  NO  delta  bands 
(N  +  0  ->  NO),  and  the  N02  continuum  (NO  +  0  ->  N02).  All  three 
band  systems  have  been  observed  in  the  terrestrial  nightglow.  Since 
these  emissions  mechanisms  involve  chemical  reactions,  some 
modelling  will  be  required  to  extract  the  0  density.  Ultimately, 
the  success  of  these  methods  will  depend  on  the  quality  of  the 
signal  obtained  from  these  nightglow  emissions. 


MOLECULAR  NITROGEN  DENSITY 


Pho t o e 1 e c t r on  impact  excitation  of  N2  produces  observable 
airglow  emissions  in  a  number  of  band  systems  including  the  Birge- 
Hopfield  (BH),  LBH,  VK,  Second  Positive,  and  First  Negative. 
Observations  of  these  emissions,  when  combined  with  a  photoelectron 
flux  model  as  discussed  above,  can  be  used  to  infer  the  N2  density 
on  the  dayside.  Emissions  from  N+  at  1085  S  and  N  at  1493  S  are 
produced  by  solar  photodissociation  and  ionization  of  N2  and 
observations  of  these  line  emissions  can  be  used  to  infer  the  N2 
dens ity . 


MOLECULAR  OXYGEN  DENSITY 


DAYSIDE : 

The  RAIDS  experiment  will  employ  a  number  of  techniques  to 
determine  the  02  density  on  the  dayside.  The  near  infrared 
spectrometer  (NIR)  on  RAIDS  will  measure  the  02  Atmospheric  bands 
which  are  produced  by  solar  scattering  by  02.  Since  the  02 
S chumann - Rung e  bands  are  strong  absorbers  in  the  FUV,  comparisons 
of  measured  N2  LBH  band  intensity  ratios  to  synthetic  spectra  will 
yield  the  02  density  above  120  km.  Similarly,  the  Schumann-Runge 
bands  will  attenuate  starlight  and  the  spectra  of  stars  ppssing 
through  the  RAIDS  instrument  line  of  sight  may  provide  a  means  for 
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02  density  measurement. 

On  the  dayside,  the  NO  gamma  bands  are  produced  by  resonant- 
fluorescent  scattering  of  sunlight.  Above  140  km,  the  NO  molecule 
is  in  chemical  equilibrium  between  source  and  loss  reactions  given 
by  : 


N(  AS)  +  02  ->  NO  +  0 
NO  +  N(4S)  ->  N2  +  0 

resulting  in  an  equilibrium  NO  density  given  by: 

[NO]  =  k  {  T  )  C  0  2  ]  , 

where  k(T)  is  the  temperature  dependent  ratio  of  source  and  sink 
reaction  rates.  If  the  temperature  is  known,  measurements  of  the 
NO  gamma  bands  above  140  km  can  be  used  to  infer  the  02  density. 


NIGHTSIDE  : 

As  mentioned  in  the  introduction,  0+  ions  undergo  charge 
exchange  with  02  molecules  forming  02  ions  which  can  then  recombine 
with  electrons  and  dissociate  producing  atomic  oxygen  airglow 
(principally  at  6300  R)  .  The  6300  R  emission  rate  is  proportional 
to  the  product  of  the  electron  and  02  density  and  can  be  used  to 
infer  the  02  density  at  night.  Absorption  of  starlight  in  the 
Schumann-Runge  bands  will  also  provide  a  means  for  02  remote 
sens ing  at  night . 


CONCLUSION 


The  eight  optical  instruments  that  make  up  the  RAIDS  orbiting 
observatory  will  provide  a  valuable  database  of  airglow  emission 
profiles  which  will  be  used  to  test  techniques  for  remote  sensing 
of  the  neutral  atmosphere  and  ionosphere.  All  of  the  airglow 
emissions  mentioned  above  have  been  observed  before  by  a  number  of 
investigators  and  all  of  the  remote  sensing  techniques  discussed 
here  have  been  attempted  before  -  at  least  tentatively.  RAIDS  will 
provide  a  unique  opportunity  to  test  and  intercompare  these 
techniques  with  a  global  dataset  of  airglow  observations. 
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4.9  THE  ATMOSPHERIC  DENSITY  SPECIFICATION  (ADS)  SATELLITE 
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David  C.  Kayser 
Aerospace  Corporation 


and 

Timothy  L.  Killeen 
University  of  Michigan 


1.  INTRODUCTION 

The  statistical  accuracy  of  empirical  models  has  remained  at  about 
15%  (one  sigma)  during  the  past  15  years.  Model  improvements  will 
require  parameterizations  of  upper  atmosphere  dynamics  and 
appropriate  representations  of  actual  physical  energy  inputs  (to 
replace  present  "proxy"  indicators) .  There  is  no  alternative  to 
a  satellite  experiment  to  gather  the  inputs  critical  for 
understanding  the  dynamic  processes  that  affect  the  variability  of 
the  neutral  atmosphere.  The  proposed  ADS  satellite  mission  will 
generate  a  coordinated  and  comprehensive  global  data  set  of  lower 
thermospheric  state  variables  needed  to  improve  the  performance  of 
operational  density  models. 


2 .  OBJECTIVE 

The  purpose  of  this  proposed  low  altitude  satellite  experiment  is 
to  obtain  coordinated  measurements  of  aerodynamic  drag,  density, 
neutral  composition,  winds  and  temperature  to: 

a)  develop  new,  more  accurate  operational  satellite  drag 
models  by  measuring  atmospheric  dynamic  processes  and  relating  them 
to  high  latitude  energy  inputs 
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b)  meet  AF  requirements  for  10%  absolute  density 

specification  by  advancing  state-of -the-art  sensing  techniques 
to  measure  neutral  density  within  5%. 


3 .  INSTRUMENTATION 

The  instrumentation  consists  of  the  AFGL  Satellite  Electrostatic 
Triaxial  Accelerometer  (SETA)  and  Quadrupole  Ion/Neutral  Mass 
Spectrometer  (QINMS) ,  the  University  of  Michigan  Fabry-Perot 
Interferometer  (FPI)  and  the  Aerospace  Corporation  Composition  and 
Density  Sensor  (CADS) .  This  coordinated  payload  provides  the 
inputs  critical  for  understanding  the  dynamic  processes  that  affect 
the  variability  of  the  neutral  atmosphere.  A  brief  description  of 
each  experiment  is  given  below: 

a)  Satellite  Electrostatic  Triaxial  Accelerometer  (SETA) 

This  instrument  measures  satellite  orbital  accelerations  at 

micro-g  levels.  It  is  an  outgrowth  of  the  single  axis  accelerometer 
systems  flown  by  AFGL  on  the  three  Atmosphere  Explorer  satellites 
and  on  six  Air  Force  satellites.  The  triaxial  version  has  been 
successfully  flown  on  Air  Force  satellites  in  1979,  1982  and  1983, 
and  on  the  S85-1  spacecraft  in  1984.  The  accelerometer  sensing 
element  consists  of  an  electrostatically  suspended  proof  mass. 
Input  accelerations  cause  this  proof  mass  to  be  displaced  slightly 
from  its  null  or  zero-acceleration  position.  Movement  of  the  order 
of  a  millionth  of  a  centimeter  causes  electrostatic  forces  equal 
and  opposite  to  the  applied  acceleration  to  re-center  the  proof 
mass  along  three  orthogonal  axes.  The  output  is  a  digital  pulse 
rate  directly  proportional  to  the  input  acceleration.  High 
accuracy  atmospheric  density  and  neutral  winds  are  obtained  from 
the  along-track  and  cross-track  accelerations  respectively. 

b)  Quadrupole  Ion-Neutral  Mass  Spectrometer  (QINMS) 

A  major  feature  of  this  experiment  is  its  capability  to  obtain 
the  first  satellite  high  accuracy  total  mass  density  measurements. 
The  objective  is  to  measure  the  density  to  within  five  percent. 
The  QINMS  sensor  is  closely  related  to  the  instrument  successfully 
employed  on  STS-4  for  ion  and  neutral  spacecraft  contamination 
studies  and  to  three  other  instruments  awaiting  space  shuttle 
flights.  It  consists  of  a  self-contained  vacuum  system  containing 
a  switchable  electron  impact  ion  source,  a  total  density  sampling 
grid,  a  retarding  potential  analyzer,  a  quadrupole  mass  selection 
filter,  a  high  sensitivity  ion  detector/electron  multiplier 
combination,  and  a  wide  dynamic  range  logarithmic  electrometer 
amplifier.  After  launch,  a  commandable  cap  opening  mechanism 
exposes  the  instrument  to  the  ambient  atmosphere  to  initiate  the 
data  collection  process. 

This  versatile  package  is  capable  of  providing  measurements 
of  neutral  density,  neutral  composition,  ion  density,  and  ion 
composition  at  a  very  fast  sample  rate  (80  meter  spatial  resolution 
composition  measurements  and  8  meter  spatial  resolution  density 
measurements) .  Gases  with  molecular  weights  between  1  and  64 
atomic  mass  units  can  be  measured  over  an  altitude  range  of  140  km 
to  more  than  500  km.  Neutral  measurements  are  made  possible  by 
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electron  impact  ionization  of  the  incident  neutral  gas  and 
subsequent  analyses  of  the  resulting  ions.  An  accommodation  sphere 
surrounding  the  ion  source  is  used  to  thermalize  the  energetic 
incoming  gas  molecules  to  a  known  temperature  and  molecular  state 
and  thereby  allow  highly  accurate  measurements  of  neutral  mass 
density.  A  movable  flag  within  the  accommodation  sphere  allows 
unperturbed  measurements  of  the  reactive  ambient  gasses.  The 
instrument  is  highly  programmable;  it  permits  measurement  of  peaks 
associated  with  an  ambient  gas  (e.g.  0,  02,  N2)  at  a  100  times  per 
second  rate.  A  programmable  and  commandable  ROM  allows  tailoring 
the  sampling  program  to  retarding  potential  analyzer  scan,  ion  mode 
mass  scan,  or  neutral  mode  mass  scan  while  simultaneously  sampling 
total  ion  density  or  neutral  pressure. 

c)  Composition  and  Density  Sensor  (CADS) 

The  CADS  design  incorporates  flight-proven  concepts  to  measure 
in-situ  densities  of  neutral  thermospheric  reactive  species  such 
as,  atomic  oxygen  and  nitrogen,  nonreactive  gases  such  as  ,  He,  N2 
and  Ar,  the  neutral  gas  temperature  and  the  neutral  wind  component 
along  the  satellite  velocity  vector.  Composition  of  the  major 
species  can  be  measured  from  perigee  to  above  4  00  km.  This 
measurement  is  accomplished  through  the  electron-impact  ionization 
of  the  incoming  neutral  gas  followed  by  mass  filtering  and  ion 
detection.  Temperature  and  wind  are  derived  from  a  special 
technique  involving  a  retarding  potential  analysis  of  ion  energy. 

The  CADS  sensor  consists  of  a  quadrupole  mass  filter, 
electron-impact  ionization  ion  source  with  two  independent  electron 
guns,  a  low  sensitivity  electrometer  ion  detector,  and  a  high 
sensitivity  spiraltron  ion  detector — all  enclosed  in  a  vacuum 
housing.  Once  on  orbit  a  cutter  mechanism  opens  the  sensor  to  the 
ambient  atmosphere.  Several  measurement  sequences  can  be  invoked 
through  the  command  system.  Among  these  are:  scan  from  3  to  48  AMU, 
dwell  on  N2  for  composition  and  energy  analysis,  visit  a  selected 
set  of  important  masses.  The  chosen  sequence  is  repeated  every 
few  seconds,  thereby  providing  high  along-the-track  resolution. 

d)  Fabry-Perot  Interferometer  (FPI) 

The  FPI  is  a  remote  sensing  instrument  designed  to  measure  the 
meridional  neutral  wind  and  the  temperature  of  the  upper 
atmosphere.  Optical  measurements  of  thermospheric  emission 
features  provide  winds  to  5  m/sec  from  Doppler  line  shifts  and 
temperature  to  10  K  from  the  line  broadening.  The  measurements  are 
made  with  an  extremely  stable,  high-resolution  Fabry-Perot  etalon. 
A  wavelength  analysis  is  performed  on  the  detected  light  by 
spatially  scanning  the  interference  fringe  plane  with  a  multi¬ 
channel  array  detector.  A  commandable  horizon  scan  mirror  permits 
viewing  the  earth's  atmosphere  at  various  tangent  heights, 
providing  a  cross-sectional  view  of  the  thermodynamic  and  dynamic 
state  of  the  thermosphere  in  16  steps,  over  a  200  km  vertical 
region. 

The  etalon  is  the  central  component  of  the  instrument.  It 
performs  the  high  resolution  spectral  analysis  by  multiple  path 
interference  of  the  transmitted  light.  This  etalon  consists  of 
two  fused  silica  substrates  upon  each  of  which  is  deposited  a  semi- 
reflecting  dielectric  multi-layer.  These  mirrors  are  flat  to  1/100 
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wavelength.  Other  major  components  of  the  FPI  are  the  main 
scanning  telescope  and  baffle,  filter  wheel,  atomic  emission  line 
calibration  source,  etalon  objective  telescope  lenses,  and  the 
image  plane  detector  that  converts  the  intensity  of  light  falling 
on  the  image  plane  into  a  set  of  discrete  electron  pulses. 

A  more  detailed  description  of  this  experiment  is  given  in 
this  section’s  paper  "Fabry-Perot  Interferometer  for  ADS"  by 
Timothy  L.  Killeen. 


4 .  ORBITAL  REQUIREMENTS 

A  low  altitude,  high  inclination  mission  is  required  to  make 
progress  in  modeling  the  atmospheric  region  below  250  km.  A 
perigee  as  low  as  140  km  is  desired  with  a  relatively  low  apogee 
(preferably  about  500  km)  to  permit  data  collection  over  a  wide 
latitude  region  on  each  pass.  Nominal  lifetime  desired  is  fifteen 
months  and  the  minimum  is  seven  months.  The  low  orbit,  long 
lifetime  requirements  are  mutually  exclusive  without  an  on-board 
propulsion  system.  It  is  recognized  that  tradeoffs  in  perigee  and 
apogee  may  be  necessary.  A  three-axis  stabilized  satellite  is 
preferable  to  maximize  duration  and  accuracy  of  the  sensor 
measurements . 


5.  SUMMARY 

The  ADS  satellite  experiment  complement  consists  of  the  following 
flight-  proven  sensors:  a)  Satellite  Electrostatic  Triaxial 
Accelerometer  (SETA)  to  measure  in-situ  drag,  density  and  cross¬ 
track  winds,  b)  Quadrupole  Ion/ Neutral  Mass  Spectrometer  (QINMS) 
to  measure  in-situ  absolute  density  within  5%  and  neutral  and  ion 
composition,  c)  Composition  and  Density  Sensor  (CADS)  to  measure 
in-situ  along-track  winds  and  neutral  composition  and  d)  Fabry- 
Perot  Interferometer  (FPI)  to  remotely  measure  along-track  winds 
and  neutral  temperature.  A  low  perigee  is  necessary  to  provide  data 
in  the  altitude  regions  of  interest  while  the  high  inclination 
permits  global  coverage.  This  well-designed  and  focused  mission 
to  determine  thermospheric  state  variables  with  high  absolute 
accuracy  and  appropriate  geophysical  coverage  is  critically  needed 
to  significantly  improve  our  density  specification  capability. 


4  -  116 


4.  10 


ACCELEROMETER  MEASUREMENTS  OF  DENSITY  AND  WINDS  FOR  ADS 

FRANK  A.  MARCOS 

AIR  FORCE  GEOPHYSICS  LABORATORY 

AND 

ANDREW  J.  MAZZELLA,  JR. 

RDP,  INC. 


1.  INTRODUCTION 

The  Satellite  Electrostatic  Accelerometer  (SETA)  experiment  is 
designed  to  measure  density  and  winds  on  the  low-perigee  altitude 
ADS  satellite.  This  instrument  has  been  flown  successfully  on  Air 
Force  satellites  in  1979,  1982,  1983  and  1984.  It  has  provided  our 
only  extensive  lower  thermosphere  neutral  density  data  base  during 
Solar  Cycle  21.  The  ADS  mission  extends  these  measurements  into 
the  present  solar  cycle. 


2.  EXPERIMENT  DESCRIPTION 

The  accelerometer  consists  of  a  proof  mass  in  the  shape  of  a 
flanged  hollow  cylinder,  a  temperature  monitor  and  supporting  power 
and  signal  conditioning  electronics.  A  complete  description  of  the 
instrument  operating  principles  is  given  by  Marcos  and  Swift 
(1982).  The  one  gram  proof  mass  is  constrained  in  all  three  axes 
by  means  of  electrostatic  forces  generated  by  electrodes  (Fig  1) . 
The  proof  mass  cylinder  is  0.025  cm  thick  and  its  flange  is  0.0375 
cm  thick.  The  spacer  has  a  thickness  of  0.056  cm.  The  gap  between 
the  flange  and  each  Z-axis  carrier  is  0.089  cm.  Between  the  X-  and 
Y-  axes  carrier  electrodes  and  the  inside  of  the  proof  mass 
cylinder  the  gap  is  0.033  cm.  The  proof  mass  is  suspended 
electrostatically  along  its  X  and  Y  (radial)  axes  by  the  eight 
electrodes  on  the  cylindrical  carrier  extending  through  the  proof 
mass.  In  the  Z  (longitudinal)  axis  its  position  is  controlled  by 
the  three  annular  electrodes  on  each  side  of  the  proof  mass 
housing.  The  relative  motion  between  the  proof  mass  and  these 
electrodes  is  sensed  in  three  orthogonal  directions,  and  restoring 
forces  are  generated  in  three  separate  feedback  constrainment  loops 
to  keep  the  proof  mass  centered.  Figs  2  and  3  show  the  force 
rebalance  loops  for  the  cross-track  and  along-track  axes 
respectively.  Each  restoring  voltage  is  proportional  to  the  input 
acceleration  in  the  corresponding  direction.  Typical  accelerations 
for  atmospheric  drag  data  are  at  the  micro-g  level. 
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3.  DATA  ANALYSIS 


The  analysis  method  for  extracting  atmospheric  density  from  the  acceleration 
measurements  has  been  described  by  Marcos  and  Swift  (1982).  Fig.  4  shows  a  sample 
of  3-axis  raw  acceleration  data  obtained  forgone  orbit.  The  acceleration  data,  in  units  of 
micro-g's,  are  plotted  vs  universal  time.  Numerical  filtering  techniques  are  used  to  remove 
high-frequency  acceleration  components  which  are  due  to  vehicle  dynamics.  Fig.  5  shows 
the  result  of  applying  a  filter  passing  only  frequencies  less  than  0.01  Hz  to  the  data  of 
Fig  4.  An  example  of  the  detailed  density  structure  obtained  from  measurements  on 
these  low  altitude  satellites  is  given  in  Fig.  6.  These  data  were  obtained  during  one  orbit 
on  14  July  82,  a  period  of  very  enhanced  geomagnetic  activity.  Coverage  is  from  34 
degrees  north  latitude  on  the  nightside  (local  time  of  2200  hr),  up  to  83  degrees  north 
latitude  and  down  to  57  degrees  south  latitude  on  the  dayside  (local  time  0930  hrs).  This 
coverage  is  typical  of  the  density  data  obtained  on  each  orbit. 

To  illustrate  the  derivation  of  winds,  wind  measurements  during  an  active  geomagnetic 
period,  when  Kp  was  7-,  are  described  (Marcos  and  Forbes,  1985).  Fig.  7  shows  raw 
acceleration  data  from  the  X  axis  plotted  against  universal  time,  geographic  and 
geomagnetic  latitude,  local  time  and  altitude.  Atmospheric  rotation  velocity  is  included  in 
the  data.  The  result  of  applying  the  numerical  filtering  is  shown  in  Fig  0.  Neutral 
horizontal  winds  derived  from  these  data  are  shown  in  Fig  9.  Unsmoothed  data  are 
shown  as  a  light  line.  To  account  for  very  low  frequency  small-amplitude  oscillations  in 
the  raw  data,  the  heavy  line  was  added  by  Marcos  and  Forbes.  These  results 
demonstrated  that  information  concerning  the  general  positions  and  magnitudes  of  the 
nighttime  and  morning  cells  characteristic  of  the  polar  thermosphere  wind  system  can  be 
derived  from  the  SETA  data. 


4.  EXPERIMENT  OBJECTIVES 

Specific  goals  of  the  SETA  experiment  on  ADS  are  to  measure  atmospheric  density  from 
perigee  to  about  300  km  and  to  measure  cross-track  winds  from  perigee  to  about  200 
km.  The  density  data  will  be  used  as  follows: 

a.  To  extend  the  density  measurements  made  by  this  instrument  during  the 
previous  solar  cycle  into  solar  cycle  22. 

b.  To  provide  a  cross-check  on  the  new  absolute  density  mass  spectrometer 
measurements.  Due  to  its  extensive  flight  history,  the  errors  in  extracting  density  from  the 
accelerometer  are  well  understood.  Comparison  with  the  mass  spectrometer  density 
values  can  provide  confidence  in  this  new  technique.  Also,  the  orbit-to-orbit  variations 
will  be  accurate  to  within  a  few  percent,  except  where  in-track  winds  are  significant 
(greater  than  about  200  meters  per  second). 

c.  To  provide  a  cross-check  on  the  total  mass  densities  derived  from  summing 
up  mass  spectrometer  individual  constituent  densities  obtained  in  the  neutral  composition 
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mode.  This  technique  is  less  accurate  than  densities  obtained  in  the  absolute  density 
mode. 

d.  To  attempt  to  derive  the  spacecraft  drag  coefficient  as  a  function  of  altitude  and 
composition.  With  measured  drag  acceleration,  absolute  density  and  winds,  and 
knowing  the  area-to-mass  ratio,  the  only  unknown  in  the  drag  equation  is  the  drag 
coefficient.  This  approach  assumes  that  the  absolute  density  accuracies  are 
independently  determined.  Thus  ADS  provides  a  unique  opportunity  to  determine  this 
important  spacecraft  parameter. 


5.  SUMMARY 

The  SETA  experiment  can  measure  atmospheric  neutral  density  and  winds.  For  the  1979 
flight,  drag  accelerations  (proportional  to  density)  were  sensed  to  apogee,  about  280  km. 
Winds  were  determined  to  a  height  of  about  200  km.  The  altitude  limits  of  useful  ADS 
data  will  depend  on  the  satellite's  area-to-mass  ratio,  vehicle  dynamic  "noise",  perigee 
altitude  and  level  of  solar  activity. 
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Figure  4.  Haw  Data  for  Test  Flight 
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Figure  b.  Haw  Data  for  One  Orbit.  Day  7  90114 
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Table  5.  rhe  stagnation  ratio  Measured  for  various 
spectrometers. 
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Verify  reproducibility  and  time  stability 


Summary 
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In  Situ  Measurements  Of  The  Neutral  Wind  And 
Their  Application  To  The  ADS  Mission 

(D.  C.  Kayser,  The  Aerospace  Corporation) 


The  upper  atmosphere  composition  spectrometer  (UACS)  was  flown  on  the 
Department  of  Defense  Space  Test  Program  mission  S85-1.  UACS  was  designed  to 
measure  neutral  composition,  temperature,  and  the  in-track  wind  component  in  the  lower 
thermosphere.  In  the  high  inclination  orbit  of  S85-1,  UACS  measured  the  meridional 
component  of  wind  at  all  but  the  highest  latitudes.  A  general  discussion  of  UACS  can  be 
found  in  Kayser,  et  al.  [1986].  A  later  paper  [Kayser,  1988]  is  devoted  to  a  description 
of  the  specific  instrumental  technique  which  is  used  to  derive  the  in-track  component  of 
the  neutral  wind. 

As  part  of  the  ADS  mission,  UACS  has  been  renamed  CADS  (for  composition 
and  density  sensor)  and  will  contribute  to  the  coordinated  measurement  of  density  and 
composition.  Its  special  contribution  will  be  in  the  area  of  temperature  and  wind 
measurements.  Later  in  this  discussion  it  will  be  shown  that  measurements  of  the  in¬ 
track  wind  component  must  be  used  in  the  analysis  of  density  and  composition  data  if  the 
goal  of  5%  accuracy  is  to  be  reached.  Also,  in-track  wind  in  combination  with  cross¬ 
track  wind  measured  by  the  accelerometer  prescribes  the  horizontal  neutral  wind  field. 
These  data  can  be  compared  with  the  output  of  dynamical  models,  thus  contributing  to 
improvements  in  the  formulation  of  such  models. 

Briefly,  in-track  neutral  wind  is  measured  by  CADS  as  an  energy  offset  of  the 
retarding  potential  analysis  (RPA)  curve.  Figure  1  is  an  example  of  such  a  curve,  and 
shows  the  normalized  detector  signal  (data  are  circles)  for  molecular  nitrogen  versus 
retarding  potential  voltage.  The  position  in  volts  of  the  50%  point  on  the  curve 
corresponds  approximately  to  the  average  energy  of  a  molecule  as  seen  in  the  satellite  (~8 
km  s1)  reference  frame.  An  in-track  wind  component,  if  it  is  present,  will  shift  the 
RPA  curve  to  higher  (head  wind)  or  lower  (tail  wind)  energies.  In  the  figure  a  best-fit 
curve  is  shown  along  with  two  curves  which  are  displaced  0.1  V  either  side  of  the  best- 
fit  voltage.  These  curves  are  calculated  from  the  neutral  gas  distribution  function  which 
is  convoluted  with  the  instrument  response  function. 
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Table  1  shows  the  theoretical  change  in  position  for  the  50%  point  on  the  RPA 
curve  for  a  range  of  head-wind  velocity,  W.  With  ~44  m/s  corresponding  to  an  easily 
detected  0.1  volt,  winds  of  this  magnitude  or  less  can  be  readily  measured.  Sensitivity 
varies  by  gas  specie  (Table  2),  due  to  the  mass  dependence  of  kinetic  energy.  In  terms  of 
millivolts  energy  shift  per  m/s  of  wind,  heavier  gases  give  greater  sensitivity. 

Abundance  of  each  specie  is  also  important  however,  as  the  sensor  output  signal  is 
proportional  ambient  density.  Combined,  these  considerations  make  molecular  nitrogen 
the  best  gas  specie  to  use  in  the  lower  thermosphere,  with  atomic  oxygen  a  possible 
candidate  for  higher  altitudes. 

The  wind  measurements  on  S85-1  were  obtained  from  analysis  of  the  nitrogen 
signal.  The  experiment  obtained  a  small  but  significant  set  of  wind  data  for  2  days  in 
July  1984  at  altitudes  between  200  and  235  km,  at  local  solar  times  of  1030  and  2230 
hours.  The  measured  winds  compare  favorably  with  published  results  from  ground-based 
incoherent  scatter  radar  and  Fabry-Perot  interferometer  experiments  and  with  satellite- 
based  cross-track  measurements  of  the  meridional  wind:  In  Figure  2  the  new  data  have 
been  averaged  into  2-degree  latitude  bins  and  plotted  using  (+)  symbols.  Also  plotted  are 
combined  data  sets  of  the  nighttime  meridional  wind  for  the  period  corresponding  to 
summer  in  the  northern  hemisphere.  Letters,  and  letters  with  vertical  bars,  are  previously 
published  data  from  the  AE-E  satellite  and  a  variety  of  ground  stations.  As  much  as 
possible  the  data  shown  are  associated  with  solar  minimum  conditions.  Winds  to  the 
north  are  positive.  The  same  format  is  used  in  Figure  3  except  that  the  combined  data 
sets  are  for  the  daytime  meridional  wind.  S85-1  measurements  are  plotted  using  dashes 
and  circles.  Overall,  these  favorable  comparisons  provide  assurance  that  the  RPA 
technique  can  be  used  effectively  to  measure  neutral  wind  under  typical  thermospheric 
conditions. 

Some  additional  comments  can  be  made  regarding  the  accuracy  of  in-track  wind 
measurements:  Table  3  summarizes  the  error  levels  realized  on  the  S85-1  mission 
employing  the  RPA  technique.  Using  only  4%  of  the  available  signal,  UACS  provided 
wind  data  with  estimated  point  errors  of  16-21  m/s  in  the  200-225  km  altitude  range.  If 
those  errors  are  adjusted  to  obtain  consistency  with  along-the-track  scatter,  values  of  20- 
28  m/s  are  calculated.  Using  typical  density  (hence  signal)  changes  between  225  and  300 
km,  one  can  estimate  an  adjusted  error  of  130  m/s  at  300  km.  Had  100  %  of  the 
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nitrogen  signal  been  used,  projected  errors  of  4,  6,  and  26  m/s  are  obtained  for  altitudes 
of  200,  225,  and  300  km  respectively.  Thus,  a  sensor  optimized  for  the  task  could  in 
principle  produce  highly  competitive  wind  measurements. 

The  extent  to  which  density  and  composition  measurements  are  affected  by  wind 
depends  on  the  specific  measurement  technique  employed.  Table  4  summarizes  the  wind 
dependence  which  is  implicit  in  each  of  the  3  common  techniques:  Drag  techniques  sense 
vehicle  drag  (A)  which  is  a  function  of  the  B-factor,  the  mass  density,  and  the  square  of 
the  sum  of  satellite  velocity  plus  in-track  wind.  Thus,  the  fractional  error  in  density  due 
to  wind  is  equal  to  2  times  the  ratio  VV/V.  A  mass  spectrometer  run  in  normal  mode 
produces  a  signal  S  which  is  a  product  of  a  calibration  factor  C,  the  ambient  composition 
N,  and  the  sum  of  satellite  velocity  plus  in-track  wind.  The  fractional  error  in 
composition  due  to  wind  is  therefore  W/V.  A  mass  spectrometer  run  in  RPA  mode 
produces  a  signal  which  to  first  order  is  independent  of  velocity  and  wind  and  is 
proportional  to  ambient  composition.  A  measurement  using  this  technique  is,  however, 
difficult  to  calibrate  and  gives  lower  sensitivity  for  the  determination  of  composition. 

Table  5  gives  the  measurement  error  in  percent  for  drag  and  normal-mode  mass  - 
spectrometer  measurements  in  the  presence  of  in-track  winds  from  10  to  500  m/s. 
Roughly  speaking,  if  an  accuracy  of  5%  or  better  is  required  from  drag  measurements 
then  winds  over  50  m/s  need  to  be  known  and  corrected  for  in  the  data  analysis.  The 
equivalent  level  for  mass-spectrometers  is  100  m/s.  Thus,  it  is  anticipated  that  CADS  can 
easily  provide  wind  data  of  sufficient  accuracy  to  remove  this  source  of  error  from 
concurrent  measurements  of  composition  and  density.  In  the  auroral  zones,  where  winds 
of  over  500  m/s  are  routinely  experienced,  this  capability  is  especially  important. 
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Figure  1.  A  retarding  potential  analysis  curve. 
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Figure  2.  S85-1  wind  data  compared  to  other  techniques. 
Nighttime  meridional  winds. 
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Figure  3.  S85-1  wind  data  compared  to  other  techniques 
Daytime  meridional  winds. 


Table  1.  Retarding  potential  voltage  shift  for  a  range  of 
of  typical  neutral  wind  values. 
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Table  2 . 


RPA  sensitivity  using  various  gases. 
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Table  3.  In-track  wind  component  measurement  errors  (given 
in  m/s) . 
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Table  4 


Density  and  composition  measurement  error  due  to 
an  in-track  wind. 


(V-  SATELLITE  VELOCITY,  W-  IN-TRACK  WIND) 

DENSITY  FROM  DRAG: 

A  =  1/2  B  p  (V+w)2 
d p/p  ~  2  w/v 

COMPOSITION  FROM  MASS  SPECTROMETER  (NORMAL  MODE): 
S  =  C  N  (V+W) 
dN/N  -  W/V 

COMPOSITION  FROM  MASS  SPECTROMETER  (RPA  MODE): 

S  =  C  N 


dN/N  ~  0 


Table  5.  Density  and  composition  measurement  error  as  a 
function  of  wind  speed. 
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4.13  FABRY-PEROT  INTERFEROMETER  FOR  ADS 
T.  L.  Killeen 

Space  Physics  Research  Laboratory 
University  of  Michigan 


The  Fabry-Pcrol  Interferometer,  FPI,  measures  the  intensity  and  spectral  location  of  visible 
atmospheric  emissions  in  altitude  regions  centered  ai  -98,  185  and  225  km.  From  these  measurements  the 
primary  FPI  state  variables,  winds  and  temperatures,  can  be  calculated  to  better  than  ±  5  m/s.  Secondary 
measurements  are  meiasiable  atom  densities  to  ±  10%.  Density  profiles  can  be  obtained  by  using  measured 
temperature  profiles  and  the  hydrostatic  approximation,  together  with  an  in-silu  observation  of  density  to 
normalize  the  profile.  Since  FPI  is  a  remote  sensing  instrument,  measurements  can  be  made  continuously, 
however  data  will  degrade  at  the  higher  altitudes  of  an  elliptical  orbit.  The  FPI  was  flown  successfully  on 
Dynamics  Explorer  2  (Hays  et  ah,  1984)-  Since  the  FPI  is  an  instrument  that  has  not  been  flown 
previously  on  an  Air  Force  spacecraft  arid  since  it  provides  critical  wind  and  temperature  information  needed 
for  the  modeling  efforts,  we  discuss  ihe  insimmcnt  in  greater  detail  than  the  in  situ  sensors  described  above. 
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FPI  Instrument  Description 

The  instrument  is  composed  of  three  basic  elements:  the  interferometer,  an  electronics  package  ;md  the 
telescope  sub-assembly.  A  schematic  of  the  optical  path  of  the  instrument  designed  to  meet  the 
requirements  discussed  above  is  shown  in  figure  1  and  the  telescopes  are  shown  in  figure  2.  The  basic 
simplicity  of  the  instrument  is  evident  here.  Also,  shown  is  the  IPD  Photometer  which  will  be  discussed  in 
the  following  section. 

Incoming  light  passes  through  both  fore  and  aft  telescopes.  The  telescopes  tire  oriented  at  -45°  with 
respect  to  the  satellite  velocity  vector.  The  two  look  directions  provide  two  vector  wind  components  and, 
since  vertical  winds  arc  normally  small  compared  to  horizontal  winds,  this  permits  the  measurement  of  the 
vector  winds.  A  simple  scan  mirror  in  each  look  direction  chooses  the  altitude.  A  fore/aft  selector  chooses 
which  telescope  signal  will  be  passed  through  a  simple  relay  optics  system  to  the  interferometer.  The 
telescopes  are  copies  of  die  DE-FPI  design. 

Two  line-to-circle  convertors  provide  entry  to  the  FPI/IPDPH.  The  line  input  matches  the  need  to 
provide  a  narrow,  horizontal  sample  of  the  horizon.  In  this  way,  good  vertical  resolution  can  be  obtained 
while  retaining  good  sensitivity  using  a  broader  horizontal  l ic'd  of  view.  The  line  shape  is  then  converted  to 
the  circular  input  required  by  both  the  F'PI  and  tile  IPDPH.  Additionally  the  line  inputs  can  be  located  close 
enough  together  such  that  they  can  share  telescopes  and  operate  simultaneously,  sampling  slightly  different 
altitude  regimes. 

The  light  from  the  linc-io-circle  convertor  is  then  partially  collimated  to  pass  through  a  dielectric 
iniorfcictice  filter  mounted  in  a  filter  wheel.  The  interferometer  light  beam  is  then  expanded  and  further 
collimated  to  meet  theoretical  requirements  of  the  fixed-gap,  high-rcsoluuon  etalon.  The  etalon  objective 
telescope  then  images  the  light,  which  now  contains  spectral  information  in  an  axially  concentric  pattern, 
onto  the  image  plane  detector.  Figure  3  shows  the  basic  interferometer  output  in  the  form  of  a  high 
resolution  12  channel  spectrogram  which  contains  the  wind,  temperature  and  brightness  information. 

The  optical  throughput  or  etendue  is  held  constant  throughout  the  optical  system  to  maximize 
sensitivity  and  minimize  size.  The  etendue  is  controlled  by  the  spectral  resolution  and  instrument  size 
constraints.  Folded  optics  arc  used  to  minimize  size  as  well  as  to  avoid  mechanical  and  thermal  distortions. 

Fubry-Peroi  Interferometer 

The  heart  of  the  FPI  is  the  Fabry-Peroi  interferometer  that  performs  the  high  resolution  spectral  analysis 
of  light  emitted  from  die  earth's  atmosphere.  The  interferometer  is  a  single,  fixed  gap  etalon  system  widi  a 
folded  optics  objective  telescope  and  a  twelve-channel  image  plane  detector.  The  interferometer  section  is 
made  up  enurcly  of  flight-spare  hardware  and  optics  from  the  Dynamics  Explorer  program. 

The  etalon  and  etalon  mount  are  the  products  of  the  development  effort  leading  to  the  successful  DE-FPI 
instrument.  The  cialon  has  three  zerodur  spacers  glued  to  two  fused  silica  end  plates  on  whose  inner  surfaces 
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arc  deposited  ihe  semi-refiecting  dielectric  coatings.  The  etalon  is  mounted  in  a  fully  kinematic  mount 
designed  to  provide  isolation  from  stresses  due  to  temperature  gradients,  etc.,  while  maintaining  optical 
alignment  and  instrumental  ruggedness.  (Killeen  et  a).,  l'J82). 
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IPD  PHOTOMETER  CHANNEL 


FABRY-PEROT  INTERFEROMETER 
AND  IPD  PHOTOMETER 


Figure  I.  Optical  schematic  of  ADS  Fabry-Peroi  interferometer 


Figure  2.  Baffle  schematic  for  ADS  Fabry-Perot  interferometer 
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INTERFER0GRAM8  OBTAINED  WITH  THE  HE/NE 
LASER  AND  DE  ETALON  AND  DETECTOR.  A 
PRESSURE  STEP  WAS  TAKEN  TO  SIMULATE  A 


INTERFEROMETER  SPECTROGRAM  OUTPUT 


Figure  3.  ADS  Fabry-Perot  interferometer  spectrogram  output.  The  peak  detected  by  the  multichannel 
IPD  contains  the  Doppler  information  necessary  to  calculate  winds,  temperature  and  volume  emission 
rates.  The  effect  of  a  100  m/sec  wind  on  the  spectrum  is  illustrated. 
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The  image  plane  detector  (IPD)  is  similarly  a  product  of  the  DE  development  program  and  is  a  12 
channel  concentric  ring  anode  device  incorporating  a  micro-channel  plate  electron  multiplication  stage  and 
an  S-20  photocathode.  Each  anode  ring  is  single  photon  counting  and  measures  one  spectral  element.  The 
LPD  effectively  "scans"  the  etulon  in  angle  without  the  need  to  mechanically  change  any  etalon  parameter.  It 
also  detects  light  in  all  spectral  elements  simultaneously,  thereby  enhancing  instrument  sensitivity  over 
conventionally  scanned  interferometers.The  IPD  is  identical  to  that  used  successfully  on  DE-FPI  (Killeen  et 
al.,  1983). 

Scanning  telescope  optical  design. 

The  telescope  is  shown  schematically  in  figure  2.  Two  look  directions  at  45  0  with  respect  to  the 
velocity  vector  are  provided  to  obtain  full  vector  wind  measurements.  The  basic  telescope  design  is  die  same 
as  DE-FPI  consisting  of  a  protective  cover,  light  shielding  baffle,  and  horizon  scan  mirror.  For  this 
mission,  an  optical  switch  selects  fore  and  aft  views. 

Electronics 

The  electronics  block  diagram  is  shown  in  figure  4.The  instrument  is  functionally  divided  into  four 
sections:  telescope,  detector,  instrument  electronics  and  satellite  interface.  The  telescope  drive  consists  of 
two  stepper  motors  for  horizon  scan  and  another  stepper  motor  for  fore  and  aft  selection.  Scan  will  be 
controlled  by  a  state  sequencer  simdar  to  that  flown  on  DE-FPI.  The  detector  electronics  are  a  basic  copy  of 
those  flown  on  DE-FPI  and  consist  of  a  HV  supply  for  the  IPD,  12  channels  of  analog  and  pulse  detection 
of  the  IPD  output  signals,  calibration  sources  and  filter  wheel  drive  and  position.  The  IPDPH  detector 
electronics  are  identical  to  die  FPI  detector  electronics.  The  instrument  electronics  sub-system  distributes 
supply  voltages,  commands  and  uming  signals  to  the  other  instrument  blocks.  Instrument  monitors  such  as 
temperature  measurements,  high  voltage  monitors,  scan  position  decoder  outputs  etc.,  enter  this  section  for 
transmission  to  the  telemetry  stream.  The  satellite  interface  section  is  designed  to  provide  simple  interfaces 
to  facilitate  testing,  integrauon  and  flight  operations  and  to  minimize  interconnects.  The  electronics  design 
is  almost  enurcly  a  carryover  from  the  DE-FPI,  AE-VA.E  systems  with  minor  updates  from  the 
HRDI/UARS  currently  in  final  assembly  and  lest. 

FPI  Sensitivity 

The  Low  Altitude  Doppler  interferometer  instrument  described  in  the  previous  sections  has  the  inherent 
high  throughput  of  die  Fabry-Perot  interferometer  and  has  the  spectral  selectivity  to  measure  winds  and 
temperatures  in  the  earth's  atmosphere.  The  basic  instrument  sensiuvity  is  tabulated  in  Table  1. 
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FPI  Error  Analysis 

Statistical  Errors 

Statistical  errors  in  the  measurements  of  winds,  temperatures  and  emission  line  brightnesses  are 
dependent  on  the  integration  lime  period  and  the  brightness  of  the  given  emission.  We  have  modified  the 
full  DE-FPI  simulator  computer  program  to  calculate  statistical  measurement  precisions  *br  the  ADS-FP1. 
These  calculations  involved  the  simulation  of  spectrograms  using  a  measured  instrument  function  and  the 
analysis  of  the  spectrograms  using  the  DE-FPI  flight  analysis  routines.  The  accuracies  are  depicted  in  figure 
5.  !t  ran  be  seen  that  lor  typical  daytime  brightnesses  for  the  5571 A  emission,  accuracies  of  5m/scc  and 
10  0  K  are  obtained  with  integration  periods  of  1  second.  For  typical  nighttime  brightnesses,  these 
accuracies  are  degraded  to  20  m/sec  and  50  °  respectively.  These  errors  are  sufficiently  small  to  enable  the 
generation  of  vector  altitude  profiles  as  required  for  the  proposed  investigations. 
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TYPICAL  'ONE  SHOT'  INTEGRATION 
ES  FOR  SPECIFIED  ACCURACY 


BRIGHTNESS  (RAYLEIGHS) 

irs  on  derived  winrD  and  temperatures  for  the  ADS  Fabry-Perot  interferometer. 


Systematic  Errors 

Possible  systematic  errors  can  occur  due  to  instrument  spectral  drift.  Experience  with  exhaustive 
laboratory  calibrations  and  spaceflight  for  the  DE-FPl  have  indicated  that  a  combination  of  continuous 
etalon  temperature  measurements  to  0.01°  and  periodic  in-flight  spectral  calibrations  are  sufficient  to  provide 
knowledge  of  instrument  drifts  to  ~lm/sec.  Two  further  types  of  systematic  error  arc  cf  imnortance.  These 
are  due  to  uncertainties  in  the  attitude  and  velocity  of  the  instrument  with  respect  to  the  earth  and  due  to 
calibration  uncertainties  in  the  instrument  discussed  below. 

The  spacecraft  has  an  inherent  orbital  velocity  of  about  7  km/scc  which  appears  as  a  large  DC 
background  shift  on  any  spectral  signal  attached  to  the  earth.  In  our  worst  case,  looking  at  45°  with  respect 
to  the  spacecraft  one  finds  that  an  altitude  uncertainty  of  0.05°  leads  to  a  5m/sec  error  in  the  deduced  wind. 
This  is  a  very  important  point  and  clearly  the  wind  measurements  become  degraded  for  large  uncertainties  in 
the  attitude.  The  accuracy  of  the  post-flight  knowledge  of  attitude  is  therefore  a  matter  of  importance  to  the 
wind  measurement  and  is  of  greater  concern  than  errors  in  the  determination  of  the  satellite  velocity  itself. 
Several  considerations,  however,  mitigate  the  negative  impact  of  large  attitude  errors  on  the  proposed 
science. 

1)  It  is  straightforward  to  use  the  Rayleigh  scattered  light  signal  from  the  lower  atmosphere  during  the 
day  and  the  5577A  emission  during  the  night  to  provide  for  a  rapid  and  accurate  in-flight  determination  of 
the  instrument  zenith  angle  using  simple  discrimination  electronics  to  optimize  the  field  of  view. 

2)  The  volume  emission  rate  altitude  profiles  and  the  kinetic  and  rotational  altitude  profiles  arc 
insensitive  to  errors  in  attitude  and  much  of  the  science  is  therefore  unaffected. 

3)  Many  of  the  interesting  science  problems  related  to  the  dynamics  of  the  mcsosphcre/lowcr 
thermosphere  involve  "AC"  rather  than  "DC'wind  systems  (e.g.  gravity  waves,  eddies  etc.)  and  therefore 
relative  measurements  of  winds,  which  arc  largely  unaffected  by  attitude  errors,  arc  sufficient  for  many 
purposes. 

4)  The  wind  measurements  themselves  (for  example  in  the  upper  thermosphere  at  low  latitudes  where  the 
wind  field  is  known  to  be  smoothly  varying)  will  provide  an  in-flight  indication  of  the  magnitude  of  the 
attitude  drift  for  later  comparison. 

IPDPH 

The  Image  Plane  Detector  Photometer,  IPDPH,  is  shown  in  figure  5.1.  The  IPDPH  is  a  integral  part  of 
the  ADS-FPI  package  which  shares  the  telescope,  filter  wheel  and  many  of  its  electronics  with  the  FPI. 
Also,  the  IPD  and  its  electronics  are  direct  copies  of  the  DE-FPI  1PD.  Additional  elements  are  a  second  line 
to  circle  convertor  and  three  simple  lenscs.The  IPDPH  complements  the  measurement  package  by  providing 
for  accurate  rotational  temperature  measurements  from  90-120  km  using  the  02  (0-0)  atmospheric  band 
emission  and  direct  O2  density  measurements  from  60  to  120  km  using  the  method  of  oxygen  fluorescence. 
Rotational  temperature  measurements  are  obtained  to  very  high  precision  by  detecting  the  spectral  changes 
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ill  the  shape  of  the  02  atmospheric  band  (Figure  6).  Figure  7  illustrates  the  accuracies  with  which 
rotational  temperatures  can  be  measured  as  a  function  of  altitude  in  the  daytime  thermosphere  by  using  a 
combination  of  the  observations  from  die  FP1  and  IPDPH  channels  of  the  ADS  remote  sensing  instrument. 
As  can  be  seen,  rapid  measurements  of  the  temperature  profile  to  accuracies  of  -1°  can  be  obtained  in  the 
altitude  region  90  -  120  km.  The  technique  is  limited  below  90  km  by  the  increasing  optical  depth  of  the 
(0-0)  emission  and  above  120  km  by  the  fall  off  in  CH  density.  In  this  region,  however,  the  hydrostatic  law 
can  be  used  to  derive  On  density  profiles. 
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RELATIVE  EMISSION  LINE  STRENGTHS 


7590  7600  10  20  30  40  50  60  70  80  90  7700  10  20  30  40  50  60  7770 


WAVELENGTH  (A) 


Temperature  Dependence  of  the  02  (0-0)  atmospheric 

molecular  emission  band  as  seen  from  orbit  (P-branch) 

Figure  6  Changes  in  spectral  signature  ot  the  02  (0-0)  atmospheric  band  due  to  changes  in  rotational 
temperature.  The  rotational  temperature  can  be  measured  by  ustng  two  relatively  low  sjvclt.iJ  resolution 
mea-.uiemcnts  with  the  FPI  and  LDPPH  ''h.u-,nels. 
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ALTITUDE  (km) 


INTEGRATION  TIME  (SEC) 


10  100 


NUMBER  OF  SAMPLES 

Figure  7  Number  of  0.22  seconds  integration  limes  r  •  ,ary  to  achieve  a  1°  rotnu'onal  temperature  error 
in  the  daytime  thermosphere  as  a  function  of  altitude  using  the  ADS-FP1. 


ADS-FPI  Sample  time  and  integration  period  for  a 
statistical  error  of  one  degree  on  a  measurement  of 
rotational  temperature  from  the  C>2  (0-0)  atmospheric 
band.  Below  ~70km,  the  optical  depth  of  the  emission  limits 
the  measurement  (although  the  technique  can  be  somewhat 

‘ tms  such 


uie  iiicdbui  triiitriii  \diuiuuyu  uitr  Leuuuque  ucjii  bumewna 

extended  into  the  optically  deep  region  using  algoritms  sue 
as  those  developed  for  the  UARS/HRDI  program,  Hays,  Privat 
commun i cat  ion ) 
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SURFACE  RANGE,  1000  km 


REENTRY  GROUND  TRACKS  BETWEEN  160  AND  60  km 
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LOCAL  SOLAR  TIME  COVERAGE 
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Match  measurements 


CONTROL  SURFACE  DEFLECTIONS 
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COMPARISON  OF  ACCELEROMETRY  AND  PRESSURE 
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ANNUAL  VARIABILITY  IN  THE  AMPLITUDE  OF  THE 
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COMPARISON  OF  AVERAGE  MEASURED  DENSITY 
NORMALIZED  TO  VARIOUS  ATMOSPHERE  MODELS 
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Shuttle  Accelerometers  Provide  Unique  Density  Measurements 
•  Horizontal  slice 

•Altitude  region  not  frequently  accessed 


